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0 . 
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SATURATION 

CURVE VALUES 

FOR 

STATOR 


132. 

18. 

1. 

40. 

2. 

66. 

5. 

76. 

8. 

85. 

14.5 

102 

101. 

114. 

300. 

132 

. 1000. 




SATURATION 

CURVE VALUES 

FOR POLE 



140. 

0 . 

1.6 

10 . 

2.3 

20. 

3.3 

30. 

4.2 

60. 

7.3 

70. 

9.2 

80. 

12.5 

85. 

15. 

90. 

• 

o 

CM 

100. 

40. 

108. 

100. 

112 

160. 

126. 

500. 

140. 

1000. 




SATURATION 

CURVE VALUES 

FOR YOKE 



140. 

0 . 

1.6 

10 . 

2.3 

20. 

3.3 

30. 

4.2 

60. 

7.3 

o 

• 

9.2 

• 

o 

CO 

12.5 

85. 

15. 

90. 

20. 

100. 

40. 

108. 

100. 

112 

160. 

126. 

500. 

140. 

1000. 

♦ 
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COMPUTER PROCEDURE FOR 


SINGLE OR 2 COIL 

OUTSIDE COIL LUNDELL DESIGN CALCULATIONS 

1. Clear core (no switch control). 

2. Insert output Form #1 into typewriter, set margin for correct 
output• 

3. Load pass #1 followed by input #1 (output printed plus punched cards). 

il. Reset and load pass #2 followed by output from pass #1 (output 

punched cards). 

3. Reset and load pass #3 followed by output from pass #2 (output 
punched cards). 

6. Reset and load pass #U followed by output from pass #3(output 
punched cards). 

7. Reset and load pass #3 followed by saturation curve values* 
and output from pass #U (output punched cards). 

3. Reset and load pass #6 followed by output from pass #5 (output 
printed plus punch cards). 

9. Reset and load pass #7 followed by output from pass #6 (output 
punched cards). 

10. Reset and load pass #8 followed by saturation curve values* and 
output from pass #7 (output punched cards). 

11. Reset and load pass #9 followed by output from pass #8 (output 
printed plus punched cards if no load saturation curve required). 

12. If there is punch card output from pass #9 a no load saturation curve 
is required. Insert output Form #2 into typewriter and reset margin. 
Load pass #10 followed by saturation curve values and output from 
pass #9 (output printed). 

* Saturation curves are loaded in order shown on Input Form #1. 
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ALL INPUT PARAMETERS ARE IN FORMAT F7.0 (FIG. 1) 


1 . 10 . 100 . .001 .1 .01 
0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 D 0 0 0 U 0 0 0 oil 

I 7 3 4 5 8 1 a IIOIIUI3mi«infc!iaii»aNHl»nH»30J!32U3«35»3?3|M4a4f 

m n m m m ii m i n|t 111 m|i m m|m i m|i 


22222222 
33333333 
44444444 
55555555 
66666666 
17777777 
8 8 8 8 8 8 8 8 


9 9 9 9 9 9 9|9 

i j ? < s«?r 


2222222 

3 3 3 3 3 3 3 

4 4 4 4 4 4 4 

5 5 5 5 5 5 5 

6 6 6 6 6 (6 

7 7 7 7 / 7 7 

8 8 8 8 8 8 8 


9 9 9 9 9 9) 

9 10 II 12 n li 

son 


2 2 2 2 2 2|2 2 2 2 2 2 2|2 2 2 2 2 2 2|2 2 2 2 2 2 2|2 
3 3 3 3 3 3|3 3 3 3 3 3 3|3 3 3 3 3 3 3|3 3 3 3 3 3 3|3 
44444 4|4 4444444444444444444 4|4 
55555 5|5 5555555555555555555 5|5 
66666 6|6 6666666666666666(66 H]s 
77777 7|7 77777 7 I 7 777777777777 7| 


9)9 99999 |9 99999 9|9 99993 

75 76 17 78 19 20 2*22 23 24 25 2t 27 2*79 30 31 12 33 > : 


0 0 0 0 0 ofo 

<4 45 <e 47 41 

111 111ll 


10 . 


0 0 0 0 0 0 0 
4$0 51 52 53 58 55 Mfl 
111111; 


8 3 8 8 8 *8 8 8 8 8 8 8 8 8 8 8 8 8 B 8 8 8 8 8 8 8 8 8 8 8 S 8 8 8 8 8 U 8 8fea 8 8 8 8 8 8 8 8 8 8 8 3fe fi 8 8 8 S fl 8 8 8 


3399999999999999999999H 

' 3 K 37 Jt 39 80 47 4! 13 44 45 46 17 48 49 50 51 5? 51 54 55 556/ 


FIG. 1 


1.0 


1000. 


0 0 0 0 0 
nsts9toei» 
11 HI 111 


qo 0 0 0 0 0 olo 

qw e * si s> s> ra|; 
1111111 


333333333333333333333 

444444444444444444444 


2 2 2 2 2 2|2 2 2 2 2 2 2|2 2 2 2 2 2 2|2 2 2 2 2 2 2|2 
3 3 3 3 3 3 3 
4444444 

555555555555 5|5 55555 5|5 55555 5|5 
6(666 6|E 66666 (|6 66666 6|6 66666 6|C 
7 7 7 7 7 7 7|7 7 7 7 7 7 7p 7 7 7 7 7 7|7 7 7 7 7 7 7b 


10 . 


9 9 9 9 9 99 1 

* M w *> £«’ Sb 66 b 66 69 70tl 72 73 74 75 76 H Tfl 79 SQ 


000000000 

|71 7? T3 14 75 n 77 78 79 60 

1111111111 


222222222 

333333333 

444444444 

555555555 

658666666 


9 9 9 S 9 9|9 9 9 9 9 9 9 9 9 9 


ALL SATURATION CURVE VALUES ARE IN FORMAT FI0.0 (FIG. 2) 
(ALL SATURATION CURVES MUST HAVE 5 CARDS) 


100 . 

0 0 0 0 0 0 0 0 o o|o 

> I J 4 I I / | • 

11111111111 

22222222222 

33333333333 

44444444444 

5S55S5S555S 

6666 G 666666 

77777777777 

008(01(8(81 


999999999 9|S 

I 2 3 ♦ S I 7 8 • 18 
m SOI 


10 . 


1. 


oooooooo op oooooooo o|o 

I2 0MIS1II7W19 2 ofct 22 23 24 M26272IM 


11111111111111111111 

22222222222222222222 

33333333333333333333 

44444444444444444444 

55555555555555555555 

E 66 (( 6 (( 6 ( 6 ((( 666 (t( 

77777777777777777777 


100 . 

oooooouoob 

3^31 I2I3M39B17n» Stic 
1111111111 


222222222222222222222222722222 
333333333333333333333333333333 

444444444444444444444444444444 

555555555555555555555555555555 
( 6 ( 6 ( 66666666666(6666666666666 
777777777777777777777777777777 

S8l8(((8((l(8((l88IIS0g|(tlt(ISI(l008SI((l(88(888S 

99999999999999999909 

iii 12 u h « « n is is n n 22 n ?« 23 « v n »30 ji 


10 . 
oooooooo olo 


42 41 44 49 IS 4) 41 IS sail 

I I I I 1 I 1 I 111 


9999999991999999999 9|9 99999999 3|9 

M 33 34 35 X 37 38 39 4W1 « *3 44 45 41 47 48 49 54^51 52 S3 54 55 K 57 51 51 Wjtl 

FIG. 2 


.01 


0 0 0 0 0 0 0 o o|o 

II SJ 53 54 55 34 II SI 19 

1111111111 


0000000000000000000 
M|S1 S7l3S4SSSS67£8SS 7071 7ir3 74 B7S 77 7t7SIO 
11111111111111111111 


2222222222222222222 

3333333333333333333 

4444444444444444444 
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7777777777777777777 
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SINGLE ) 

OR ) OUTSIDE COIL LUNDELL 
TWO COIL) 


Calculation Electrical Fortran 

Number Symbol Symbol 


(128) 

A 

A, a 

A 

(46) 

a c 

AC 

(153) 

a cf 

AS 

(70) 

A g2 

A2 

(124) 

A y2 

AY2 

(79) 

■p 

AP 

(20) 

B 

B, b 

BK 

(22) 

b o 

BO 

(94) 

B c> 

BC1 

(76) 

b pl 

BP1 

(76) 

b p2 

BP2 

(95) 

B g , 

BG1 

(122) 

B g 2 

BG2 

(224) 

B g2FL 

BG2L * 

(58) 

b t 


(57a) 

bt 1/3 

SM 

(91) 

B t , 

BT1 

(57) 

b tm 

TM 
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Calculation 

Number 

Electrical 

Symbol 

Fortran 

Symbol 

(126a) 

B y 

by 

(229b) 

B yL 

BYL 

(228) 

B y 2 L 

BV2L 

(15) 

W 

BV 

(125) 

B y 2 

BY2 



c 

(32) 

C 

C 

(71) 

Cl 

Cl 

(74) 

C M 

CM 

(73) 

Cp 

CP 

(75) 

C q 

CQ 

(72) 

C W 

CW 



_d_ 

(12) 

D 

DU 

(ID 

d 

PI * 

(35) 

d b 

DB 

( 11 a) 

d r 

DR 

(78) 

d g2 

DQ2 

(78) 

d y 2 

DY2. 
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Calculation 

Number 

Electrical 

Symbol 

Fortran 

Symbol 

* 


e 

(3) 

E 

EE 

(55) 

ef TOP 

ET 

(56) 

efbot 

EB 

(238) 

e ffl 

EPFL 

(127b) 

e fnl 

EPNL 

(4) 

e ph 

EP 




(5a) 

f 

F 

(98) 

F c> 

FC 

(236) 

f fl 

FFL 

(96) 

Fg, 

FG 

(123) 

F g2 

FG2 

(225) 

F g2FL 

FG2L 

(127) 

f nl 

FNL 

(104) 

F P 

FP 

(126b) 

F 

y 

FY 

(180) 

F sc 

FSC 

(97) 

f t, 

FT 

(183) 

F & W 

WF 

(126) 

F y2 

FY2 
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Calculation 

Number 

Electrical Fortran 

Symbol Symbol 

(229c) 

F yL 

FYL 

(229) 

F y2L 

FY2L 



Ajl 

(59) 

g 

GC 

(59a) 

S2 

GZ 

(69) 

Se 

GE 



H, h 

(24) 

h c 

HC 

(38) 

t 

h ST 

SD 

(39) 

h ST 

SH 




(237) 

X FFL 

AIFL 

(127a) 

X FNL 

AINL 

(8) 

!PH 

PI 

(182) 

i 2 r f 

FE L 

(241) 

i 2 r f 

FCUL 

(194) 

i 2 r 

PS 

(245) 

i 2 r l 

SCUL 
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Calculation 

Number 

Electrical 

Symbol 

Fortran 

Symbol 


A 

_k^ 

(19) 

k 

WL 

(9a) 

K c 

CK 

(43) 

Kd 

DF 

(63) 

K e 

EK 

(16) 

Ki 

SF 

(44) 

Kp 

CF 

(67) 

Ks 

CC 

(42) 

KsK 

FS 

(2) 

k va 

VA 

(61) 

K x 

FF 


_L, 

1 

(13) 

X 

CL 

(80) 

h 

A LI 

(80) 

*2 

AL2 

* 

(82) 

^3 

AL3 

(48) 

l E 

EL 

(36) 

^e2 

CE 

(161) 

l f 

SI 

(78) 

3g2 

AU52. 
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Calculation 

Number 

Electrical Fortran 

Symbol Symbol 

(76) 

'p 

ALP 

(17) 

h 

SS 

(49) 

it 

HM 

(147) 

& tf 

FE 

(78) 

l y 

ALY 



M, m 

(5) 

m 

PN 



N^n. 

(146a) 

N co 

COILS 

(146) 

n f 

PT 

(45) 

n e 

EC 

(30) 

n s 

SC 

(34) 

n st 

SN 

(34a) 

n st 

SN1 

(14) 

n v 

HV 



Ajl 

(6) 

P 

4 

PX 

(9) 

PF 

PF 

(80) 

P 1 

PI 

(81) 

p 2 

P2 

(82) 

P 3 

P3 
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Calculation 

Number 


Electrical 

Symbol 


Fortran 

Symbol 


(83) 

p 4 

P4 

(84) 

p 5 

P5 

(86) 

p 7 

P7 



q 

(23) 

Q 

QQ 

(25) 

q 

QN 



R^ 

(154) 

% (cold) 

FK 

(155) 

% (hot) 

FR 

(7) 

RPM 

RPM 

(53) 

R SPH(cold) RG 

(54) 

r RPH (hot) RP 



S, s 

(181) 

SCR 

SCR 

(127c) 

S F 

CD 

(47) 

S S 

S 




(177) 

T a 

TA 

(178) 

T d 

T5 

(176) 

T do 

TC 
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Calculation 

Number 

(78) 

(78) 

(145) 

(185) 

(186) 
(243) 
(242) 
(184) 

(129) 

(131) 

(132) 

(133) 

(167) 

(168) 

(166) 

( 160 ) 


Electrical Fortran 

Symbol Symbol 


t y 

TV 

ty2 

TYZ 

V^JL 

v r 

VR 

W, w 

W C 

WQ 

W NPL 

WN 

W PFL 

WNL 

W TFL 

WTFL 

W TNL 

WT 

J 

X 

X 

XR 

X ad 

XD 

X aq 

XQ 

Xd 

XA 

Xd 

xs 

x d 

XX 

Xdu 

XU 

Xp 

XF 
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aiculation 

Number 

Electrical 

Symbol 

Fortran 

Symbol 

(150) 

Xf °c 

T2 

(130) 

x * 

XL 

(169) 

f ! 

XY 

(134) 

X 

q 

XB 

(50) 

x s °c 

TI 

(170) 

*2 

XN 

(172) 

*o 

xo 


JL 

JL 

(31) 

y 

YY 


0_ 


(108) 

CM 

tO 

o 

PGY 

(221) 

0g2L 

PqZL 

(100) 

01 

PL 

(118) 

0 l 5 

PL5 

(99) 

0 n 

PL7 

(226) 

$ 5L 

PL5L 

(207) 

0 7L 

PL7L 

(92) 

0 P 

FQ 


MA-18 


Calculation 

Number 

Electrical Fortran 

Symbol Symbol 

(213) 

0PL 

PQL 

(102) 

0pT 

PLT 

(88) 

0T 

TG 

(227) 

0 y2L 

PG2L 



_r 

(41) 

T 

1 P 

TP 

(26) 

r s 

TS 

(40) 

^SK 

SK 

(27) 

1/3 

TT 



A 

(64) 


EW 

(62) 

Ai 

PC 

(64a) 

Az 

yO ALZ 

(151) 


RR 

♦ 

(152) 

A (hot) 


(51) 

A 

RS 


(77) PE 
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C PASS 1 OUTSIDE COIL LUNDELL 

1 FORMAT(El 1«5,E11,5,E11.5,E11.5,E11«5,E11.5) 

2 F0RMAT(F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0) 

3 FORMAT(9X F12.5.2X F12.5) 

7 READ2,VA,EE,EP,PN,F,PX,RPM,PI,PF,CK 
READ2,POL,DI,DU,CL,HV,BV,SF,WL,BK,ZZ 
READ2,B0,B1,B2,B3,BS,HO,HX,HY,HZ,HS 
READ2,HT,HW,QQ,W,RF,SC,YY,C,DW,SN 
READ2,SN1,DW1,0B,CE,SH,SD,PBA,SK,T1,RS 
READ2,ALG2,DG2,GC,G2,C1,CW,CP,EL,CM,CQ 
READ2,PE,BP1,BP2,TP1,TP2,ALP,DR,WR,D1,P1 
READ2,P2,P3,P4,P5,P7,P8,DY2,TY2,TY,ALY 
READ2.AL1,COILS,PT,FE,RD,RT,T2,RR,SNL,WF 
READ1, ALZ 
SS-SF*(CL-HV*BV) 

HC-(DU-DI-2.0*HS)*0.5 

QN-QQ/(PX*PN) 

TS-3.142*DI/QQ 
IF(ZZ-4.0)9,10,9 

9 TT«(0.667*HS+DI)*3.142/QQ 
oa TQ U 

10 TT-3.1416*(DI+2.*H0+1.32*BS)/QQ 

11 IF(ZZ-1.0)12,12,13 

12 BO-BS 

CC-(5.*GC+BS)*TS/((5.*GC+BS)*TS-BS*BS) 

GO TO 14 

13 QC«(4.44*GC+0.75*B0)*TS 
CC-QC/(QC-BO*BO) 

14 CS-YY/(PN*QN) 
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TP«3.1416*DI/PX 
IF(SK)18,18,19 

18 FS-1.0 
GO TO 20 

19 FS«SIN(1.571*SK/TP)*TP/(1.571*SK) 

20 IF(PBA-60.)21,21,22 

21 D-1.0 

GO TO 95 

22 D-2.0 
95 l-QN 

U-1 

IF(QN-U)23,23,24 

24 U-PX*PN 
XX-U 
N-U 

DO 25 K»1,N 

Z-U/XX 

l-Z 

Z1-I 

IF(Z-Z1)26,26,25 

26 ZY-QQ/XX 
1-ZY 
Z1-I 

IF(ZY-Z1)27,27,25 

25 XX-XX-1. 

23 ZY-QH 

27 DF-SIN(.5236*D)/(ZY*D*SIN(.5236/ZY)) 
CF-SIN(YY*1.571/(PN*QN)) 
EC«QQ*SC*CF*FS/C 
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GE*=CC*GC 

IF(C1)29,28,29 

28 C1»0.649*L0G(PE)+1.359 

29 IF(CW)30,30,31 

30 CW=0.7 07*EE*C1*DF/(EP*PN) 

31 IF(CP)32,32,33 

32 CP=PE*(L0G(GC/TP)*.0378+1.191) 

33 IF(EL)34,34,42 

34 IF(RF)35,35,41 

35 IF(PX-2.0)36,36,37 

36 U=1.3 

GO TO 40 

37 IF(PX-4.0)38,38,39 

38 U-1.5 

GO TO 40 

39 IM .7 

40 EL=>3.142*U*YY*(DI+HS)/QQ+0.5 
GO TO 42 

41 EL=2.0*CE+3.142*(0.5*HX+DB)+YY*TS*TS/SQRT(TS*TS-BS*BS) 

42 AA«1.571*PE 
AB=3.142*PE 

IF (CM)43,43,44 

43 CM«(AB+SIM(AB))/(SIM(AA)*4.) 

44 IF(CQ)45,45,46 

45 CQ=(0.5*C0S(AA)+AB-SIN(AB))/(4.0*SIN(AA)) 

46 RB=(T1 +23 4.5) *0.003 94*RS 

PR INT3,SS > CC,HC,GE,TS,C1,TT,CW,FS,CP,DF,EL,CF,CM,EC,CQ 
PUNCH1 # VA,EE # EP,PN,F,PX 
PUNCH1,RPM,PI,PF,CK,P0L,DI 
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PUNCHl,DU,CL,SS,HC,SF,QN 
PUNCH1,WL,BK,ZZ,B0,B1, B2 
PUNCH1,B3,BS,H0,HX,HY,HZ 
PUNCHl,HS,HT,HW,QQ,W,RF 
PUNCHl,SC,YY,C,TS,SN,DB 
PUNCH1,CE,SH,SD,TT,SK,RB 
PUNCHl,TP,Dl,FE,RD,RT,COILS 
PUNCHt,T2,RR,SNL,WF,PE,SN1 
PUNCH1,DW1,BP1,BP2,TP1,TP2,ALP 
PUNCHl,DR,WR,DY2,TY2,TY,ALY 
PUNCHl,P1,P2,P3,P4,P5,P7 
PUNCHl,G2,DG2,ALG2,AL1,P8 
PUNCHl,GE,CS,CF,FS,EC,DF 
PUNCHl,RS,GC,PT,Cl,CW,CP 
PUNCHl,EL,CM,CQ,DW,CC,PBA 
PUNCHl,ALZ 
PAUSE 
END 


MAf23 



c 


PASS 2 OUTSIDE COIL LUNDELL 

1 FORMAT(Ell.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 
DIMENSION DA(8),DX(6),DY(8),DZ(8) 

READ 1,VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,CK,POL,DI 
READ 1,DU,CL,SS,HC,SF,QN 
READ 1,WL,BK,ZZ,B0,B1,B2 
READ 1,B3,BS,HO,HX,HY,HZ 
READ 1,HS,HT,HW,QQ,W,RF 
READ 1,SC,YY,C,TS,SN,DB 
READ 1,CE,SH,SD,TT,SK,RB 
READ 1 f TP,D1,FE,RD,RT,COILS 
READ 1,T2,RR,SNL,WF,PE,SN1 
READ 1,DW1,BP1,BP2,TP1,TP2,ALP 
READ 1,DR,WR,DY2,TY2,TY,ALY 
READ 1,P1,P2,P3,P4,P5,P7 
READ 1 t G2,DG2,ALG2,ALl,P8 
READ 1,GE,CS,CF,FS,EC,DF 
READ 1,RS,GC,PT,C1,CW,CP 
READ 1,EL,CM,CQ,DW,CC,PBA 
READ 1,ALZ 
DT-DWt 

IF(ZZ-3.0)49,50,51 

49 SM-TT-BS 
GO TO 53 

50 SM«(3.1416*(DI+2.*HS)/QQ)-B3 
GO TO 53 

51 IF(ZZ-4.0)50,52,49 

52 SM»TT-.94*BS 





53 HM-CL+EL 

IF(DT) 61,61,62 

61 AC-0. 785*DW*DW*S N1 
GO TO 72 

62 ZY-0.0 
DA(1)-0.05 
DA(2)-0.072 
DA(3)-0.125 
DA(4)«0.165 
DA(5)-0.225 
DA(6)-0.438 
DA(7)-0.688 
DA(8)-1.5 

DX(1)"0.000124 
DX(2)-0.00021 
DX(3)-0.00021 
DX(4)«0.00084 
DX(5)-0.00189 
DX(6)-0.00189 
DY(1)-0.000124 
DY(2)«0.000124 
DY(3)-0.00084 
DY(4)-0.00084 
DY(5)-0.00189 
DY(6)-0.00335 
DY(7)-0.00754 
DY(8)-0.03020 
DZ(1)«0.000124 
DZ(2)»0.000124 




DZ(3)«0.000124 
DZ(4)«0.00335 
0Z(5)“0.00335 
DZ(6)»0.00754 
DZ(7)»0.0134 
DZ(8)=0.0302 

63 IF(DT-. 05)201,201,200 

200 JA»0 
JB=0 
JC=0 
JD=0 

64 JA=JA+1 
JB=JB+1 
JC^JC+I 
JD«JD+1 

IF(DT-DA(JA))65,65,64 

201 D=0 

IF(ZY)71»71,54 

65 IF(DW-O.188)66,66,67 

66 CY-DX(JB-I) 

CZ=DX(JB) 

4 

GO TO 70 

67 IF(DW-O,75)68,68,69 

68 CY-DY(JC-I) 

CZ-DY(JC) 

GO TO 70 

69 CY-DZ(JD-I) 

CZ-DZ(JD) 

70 D»CY+(CZ—CY)*(DT~DA(JA-1))/(DA(JA)-DA(JA—1)) 





IF(ZY)71,71,54 
71 AC»(DT*DW-D)*SN1 


72 IF(RT)73,73,74 

73 AS»0.785*RD*RD 
GO TO 55 

74 ZY-1.0 
DT-RT 
DW=RD 

GO TO 63 

54 AS»RT*RD-D 

55 S»PI/(C*AC) 

CY-PT *FE*0.OOOOOI/AS 

FK»RR*CY 

FR»(T2+234.5)*FK*0.00394 
RO0.321*PT *FE*AS 
IF(SH)202,203,202 

203 ET=1 
EB«=1 

GO TO 204 

202 AA«0„584+(SN*SN-1,0)*0.0625*(SD*CL/(SH*HM))**2.0 
AB=(SH*SC*F*AC/(BS*RB))**2.0 
ET«AA*AB*0.00335+1.0 
EB~ET-0.00168*AB 

204 RY=SC*QQ*0.OOOOOI*HM/(PN*AC*C*C) 

RG=RS*RY 

RP=RB*RY 

A=PI*SC*CF/(C*TS) 

PUNCH1,VA,EE # EP,PN,F,PX 
PUNCH1,RPM,PI,PF,CK,POL,DI 


MA-27 



PUNCH1 f DU,CL,SS,HC,SF,QN 
PUNCH1,WL,BK,ZZ,B0,B1,B2 
PUNCH1,B3,BS,HO,HX,HY,HZ 
PUNCHl,HS,HT,HW,QQ,W,RF 
PUNCH1,SC,YY,C,TS f SN,DB 
PUNCHl,CE,SH,SD,TT,SK,RB 
PUNCH1,TP,D1,FE,RD,RT,COILS 

PUNCHl ,T2,RR,SNL,WF,PE,SN1 

PUNCHl,DW1,BP1,BP2,TP1,TP2,ALP 
PUNCHl,DR,WR,DY2,TY2,TY,ALY 

PUNCHI ,P1 ,P2,P3,P i +,P5,P7 

PUNCHl,G2,DG2,ALG2,AL1 
PUNCH1,GE,CS,CF,FS,EC,DF 
PUNCH1,RS,GC,PT,C1,CW,CP 
PUNCHl,EL,CM,CQ,DW,CC,PBA 

PUNCH1,HM,SM,AS,AC,ET,EB 

PUNCH1,S,FK,FR,RC,RG,RP 

PUNCHl,A,P8,ALZ 

PAUSE 

END 
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PASS 3 OUTSIDE COIL LUNDELL 
1 FORMAT(Ell.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 
READ 1,VA,EE,EP,PN,F,PX 
READI ,RPM,PI,PF,CK,POL,DI 
READ 1,DU,CL,SS,HC,SF,QN 
READ 1,WL,BK,ZZ,B0,B1,B2 
READ 1,B3,BS,HO,HX,HY,HZ 
READ 1,HS,HT,HW,QQ,W,RF 
READ 1,SC,YY,C,TS,SN,DB 
READ 1,CE,SH,SD,TT,SK,RB 
READ 1,TP,D1,FE,RD,RT,COILS 
READ 1,T2,RR,SNL,WF,PE,SN1 
READ 1,DW1,BP1,BP2,TP1,TP2,ALP 
READ 1,DR,WR,DY2,TY2,TY,ALY 
READ 1,P1,P2,P3,P4,P5,P7 
READ 1,G2,DG2,ALG2,AL1 
READ 1,GE,CS,CF,FS,EC,DF 
READ 1,RS,GC,PT,C1,CW,CP 
READ 1,EL,CM,CQ,DW,CC,PBA 
READ 1,HM,SM,AS,AC,ET,EB 
READ 1,S,FK,FR,RC,RG,RP 
READ 1,A,P8,ALZ 
IF(PBA-60.0)105,105,108 

105 IF(CS-O. 667 )1 06 ,1 06 , 1 07 

106 FF=0.25*(6.0*CS-1.0) 

107 FF=0.25*(3.*CS+1.0) 

GO TO 75 

108 IF(CF-O. 667 )109,109,11 0 

109 FF»0.05*(24.0*CS-1.0) 
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GO TO 75 
110 FF=0.75 

75 CX=FF/(CF*CF*DF*DF) 

Z-CX*20.0/(PN*QN) 

BT=3.142*0 I /QQ-BO 

Z A=BT*BT/(16.0*TS*GC) 

ZB=0.35*BT/TS 

ZC=HO/BO 

ZD«HX*0.333/BS 

ZE-HY/BS 

IF(ZZ—2.0) 76,77,78 

76 PC=Z*(ZE+ZD+ZA+ZB) 

GO TO 82 

77 PC=Z' t (ZC+(2.0*HT/(BO+BS))+(HW/BS)+ZD+ZA+ZB) 

GO TO 82 

78 IF(ZZ-4.0) 79,80,81 

79 PC=Z*(ZC+(2.0*HT/(B0+B1))+(2.0*HW/(Bl+B2))+(HX*0.333/B2)+ZA+ZB) 
GO TO 82 

80 PC=Z*(ZC+0.62) 

GO TO 82 

81 PC=Z*(ZE+ZD+(0.5*GC/TS)+(0.25*TS/GC)+0.6) 

82 EK=EL/(10.0**(0.103*YY*TS+0.402)) 

IF(DI—8.0) 83,83,84 

83 EK=SQRT(EK) 

84 ZF=.612*L0G(10.0*CS) 

EW=6.28*EK*ZF*(TP**(0.62-(0.228*L0G(ZF))))/(CL*DF*DF) 

87 ZA»3.1416*(DI+HS)/0Q 
IF(ZZ-3.0) 88,89,88 

88 TM=ZA-BS 
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GO TO 90 

89 TM=(3.!41b*(D!+2.*HS)/QQ)-B3 

90 WI=(TM*QQ*SS*HS+(DU-HC)*3.142*HC*SS)*0.283 
IF(WF)445,446,445 

446 WF-2.52E-6*(DR**2.5)*ALP*RPM**1.5 

445 WC».321*HM*QQ*AC*SC 

PUNCH1,VA,EE,EP,PN,F,PX 
PUNCH!,RPM,PI,PF,CK,POL,DI 
PUNCH!,DU,CL,SS,HC,SF,QN 
PUNCH1,WL,BK,ZZ,B0,B1,B2 
PUNCH1,B3,BS,H0,HX,HY,HZ 
PUNCH1,HS,HT,HW,QQ,W,RF 
PUNCH1,SC,YY,C,TS,SN,DB 
PUNCH1,CE,SH,SD,TT,SK,RB 
PUNCH1,TP,D!,FE,RD,RT,COILS 
PUNCH1,T2,RR,SNL,WF,PE,SN1 
PUNCH1,DW!,BP1,BP2,TP1,TP2,ALP 
PUNCH 1 ,DR,WR,DY2,TY2,TY,ALY 
PUNCH!,P!,P2,P3,P4,P5,P7 
PUNCH!,HM,SM,AS,AC,ET,EB 
PUNCH!,S,FK,FR,RC,RG,RP 
PUNCH!,G2,DG2,ALG2,AL! 

PUNCH1,GE,CS,CF,FS,EC,DF 
PUNCH1,RS,GC,PT,C!,CW,CP 
PUNCH!,EL,CM,CQ,DW,CC,PBA 
PUNCH!,FF,CX,PC,EK,EW,TM 
PUNCH! ,A,WI,WC,P8,ALZ 
PAUSE 
END 





c 


PASS 4 OUTSIDE COIL LUNDELL 
1 FORMAT(El 1.5,El 1.5,El 1.5,El 1.5,El 1.5,Ell.5) 
READ1 ,VA,EE,EP,PI'I,F,PX 
READl ,RPM,PI,PF,CK,POL,DI 
READ 1,DU,CL,SS,HC,SF,QN 
READ 1,WL,BK,ZZ,B0,B1,B2 
READ 1,B3,BS,H0,HX,HY,HZ 
READ 1,HS,HT,HW,QQ,W,RF 
READ 1,SC,YY,C,TS,SN,DB 
READ 1,CE,SH,SD,TT,SK,RB 
READ 1,TP,D1,FE,RD,RT,COILS 
READ 1,T2,RR,SNL,WF,PE,SN1 
READ 1,DW1,BP1,BP2,TP1,TP2,ALP 
READ 1,DR,WR,DY2,TY2,TY,ALY 
READ 1,P1,P2,P3,P4,P5,P7 
READ 1,HM,SM,AS,AC,ET,EB 
READ 1,S,FK,FR,RC,RG,RP 
READ 1,G2,DG2,ALG2,AL1 
READ 1,GE,CS,CF,FS,EC,DF 
READ 1,RS,GC,PT,C1,CW,CP 
READ 1,EL,CM,CQ,DW,CC,PBA 
READ 1,FF,CX,PC,EK,EW,TM 
READ 1,A,WI,WC,P8,ALZ 
IF(Pi )400,401,400 

401 P1=3.19*BP1*TP1/AL1 
400 IF(P2)402,403,402 
403 AL2*=TP—(BP1+BP2) /2. 

P2=3.19*(ALP*(TP2+TPl)/2.)/AL2 

402 IF(P3)404,405,404 
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405 R4=AL1+CL/2. 

P3=( (3.*BP1+BP2)/4.) *L0G (R4/AL1) *6.28/3.1416 
404 IF(P4)409,407,409 

407 P4=(3.19*ALP/3.1416)*L0G(1.+(BP1+BP2)/(AL2*2.)) 
IF(PX-4.)408,408,409 

408 P4=P4*1.5 

409 TG=6.E6*EE/(CW*EC*RPM) 

BT1=TG/(QQ*SS*SM) 

FQ=TG*CP/PX 

BC1=FQ/(2.*HC*SS) 

BG1«TG/(3.1416*0l*CL) 

FG»BG1*GE/.00319 
AY2=3.1416*0Y2*TY2 
A2=3.1416*DG2*ALG2 
AP=*BP2*TP2 

AY=3.1416*TY*(DU+TY) 

WQ*= (D U-HC) *1. 42*HC*S S* (BC1 /B K) **2.0*WL 
WT= SM *QQ*SS*HS*0.453*(BT1/BK)**2.0*WL 

132 D2=BG1**2.5*0.000061 
03»(0.0167*QQ*RPM)**1.65*0.000015147 
IF(TS-0.9)133,133,134 

133 D4=TS**1.285*0.81 
GO TO 137 

134 IF(TS-2.0)135,135,136 

135 D4=TS**1.145*0.79 
GO TO 137 

136 D4=TS**0.79*0.92 

137 D7“B0/GC 

IF(D7-1.7)138,138,139 
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138 D5=D7**2.31*0.3 
GO TO 144 

139 IF(D7-3.0)140,140,141 

140 D5=D7**2.0*0.35 
GO TO 144 

141 IF(D7-5.0)142,142,143 

142 D5»D7**1.4*0.625 
GO TO 144 

143 D5=D7**0.965*1.38 

144 06=10.0**(0.932*01-1.606) 

BA=3.142*0 I*CL 

WN=D1*D2*D3*D4*D5*D6*BA 
PUNCH1,VA,EE,EP,PN,F,PX 
PUNCH1,RPI1, PI,PF,CK,POL,DI 
PUNCH1,DU,CL,SS,HC,SF,QN 
PUNCH1,WL,BK,ZZ,B0,B1,B2 
PUNCH1,B3,BS,H0,HX,HY,HZ 
PUNCH 1 ,HS,HT,HW,QQ,V/,RF 
PUNCH1,SC,YY,C,TS,SN,DB 
PUNCH1,CE,SH,SD,TT,SK,RB 
PUNCH1,TP,D1,FE,RD,RT,COILS 
PUNCHI,T2,RR,SNL,WF,PE,SN1 
PUNCH1,DW1,BP1,BP2,TP1,TP2,ALP 
PUNCH1,DR,WR,DY2,TY2,TY,ALY 
PUNCH1,HM,SM,AS,AC,ET,EB 
PUNCH1,S,FK,FR,RC,RG,RP 
PUNCH1,G2,DG2,ALG2,AL1 
PUNCH1,GE,CS,CF,FS,EC,DF 
PUNCH1,RS,GC,PT,C1,CW,CP 


i 



PUNCH1,EL,CM,CQ,DW,CC,PBA 
PUNCH!,FF.CX,PC,EK,EW,TM 
PUNCH!,A,Wl ,WC,WN,WT,WQ 
PUNCH1,TG,FQ,BC1,BT1,BG1,FG 
PUNCH1,P1,P2,P3,P4,P5,P7 
PUNCH1,PX,ALP,A2,G2,AY2,AY 
PUNCH1,DU,DY2,ALY,PT,FK,AS 
PUNCHI,C0ILS,HC,AP,HS,P8,ALZ 
PAUSE 
END 
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C PASS 5 OUTSIDE COIL LUNDELL 
DIMENSION Al(90) 

1 FORMAT(El 1.5,Ell.5,El 1.5,El 1.5,Ell.5,El 1.5) 

888 FORMAT(FIO.O,FIO.O,FIO.O,FIO.O,FIO.O,FIO.O) 

K=>1 

823 READ888,Al(K),AI(K+1),AI(K+2),AI(K+3),AI(K+4),AI(K+5) 
K=K+6 

IF(K-89)823,824,824 

824 DO 825 J-1,20 

READ 1,R1,R2,R3,R4,R5,R6 

825 PUNCH1,R1,R2,R3,R4,R5,R6 
READ1 ,TG,FQ,BC1,BT1,BG1,FG 
READ1 ,P1,P2,P3,P4,P5,P7 
READl ,PX,ALP,A2,G2,AY2,AY 
READ1 ,DU,DY2,ALY,PT,FK,AS 
READ 1,COILS,HC,AP,US,P8,ALZ 
L0AD**1 

C0REL=3.1416*(DU—HC)/(4.*PX) 

X*=BT1 

NA=1 

K-1 

GO TO 802 

806 FT-HS*AT 
X=BC1 
K=2 
NA=1 

GO TO 802 

807 FC=COREL*AT 
FS=FT+FC 
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PL=(2.*(FG+FT)+FC)*(P1+P2+P3+P4)*.001*PX 

PLT=FQ+(PL*2./PX) 

BP*=PLT/AP 
X**BP 
NA=*31 
K=3 

GO TO 802 
808 FP=ALP*AT 

PL7=.001 *? 1 *(FC+FT+FG+F P) 

PL8=P8* (2.*FG+2.*FT+FC)*.001 

PG2=PLT*PX/2.+PL7+PL8 

BG2=PG2/A2 

FG2=BG2*G2/.00319 

Z=FG2+FP+FG+FT+FC 

IF(CO ILS-1.)817,818,817 

817 BY=PG2/AY 
821 NA=61 

X=BY 

K«5 

GO TO 802 
810 FY=AT*ALY 

IF(CO ILS-1.)820,822,820 
820 PL5«.001*P5*(Z+FY) 

818 IF(COILS—1.0)814,815,814 

814 BY2=(PG2+PL5)/AY2 
GO TO 816 

815 BY2-PG2/AY2 

816 X-BY2 
K-4 
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NA=61 
GO TO 002 

809 FY2=((DU-DY2)/6.)*AT 

IF(CO ILS-1.)822,826,022 

826 PL5=.002*P5*(Z+FY2) 

BY=(PG2+PL5)/AY 
GO TO 821 

822 FNL= (FY+FY2+FG2+FP+FG+FT+FC/2.)*2. 

AINL=FNL/(PT*COILS) 

CD=AINL/AS 

EPNL=AINL*FK*COILS 

PUNCH1,TG,FQ,BC1,BT1,BG1,FG 

PUNCH1,P1,P2,P3,P4,P5,P7 

PUNCH1,A2,AY2,AY,BP f BY,BY2 

PUNCH1,BG2,FG2,FNL,AINL,CD,EPNL 

PUNCHI,AP,PL,PLT,FC,FT,P8 

PUNCH1,PL8,ALZ 

PAUSE 

802 IF(AI (NA)-X)830,831,831 

831 NA=NA+3 

835 IF(A!(NA)-X)833,834,834 

833 NA=NA+2 
GO TO 835 

834 AX=A I (NA) 

BB1=AI(NA-2) 

DC=AI(NA+1) 

D=AI(NA-1) 

XX= (AX-BB1)/(.4343*(LOG(DC)-LOG(D+.0001))) 
Y=AX-XX*.4343*L0G(DC) 


MA.-38 



AT=EXP(2.306*(X-Y)/XX) 

GO TO (838,839),LOAD 
838 GO TO (806,807,808,809,810),K 
830 GO TO (836,837),LOAD 
836 PRINT 850, 

850 FORMAT (17HMACHINE SATURATED) 
PAUSE 
END 
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C PASS 6 OUTSIDE COIL LUNDELL 

3 FORMAT(9X F12.5,2X F12.5) 

870 FORMAT(9X F12.5,2X F12.5/) 

1 FORMAT(El 1.5,Ell.5,El 1.5,El 1.5,El 1.5,El 1.5) 
READ1 ,VA,EE,EP,PN,F,PX 
READ1 ,RPM,PI,PF,CK,POL,DI 
READ1 ,DU,CL,SS,HC,SF,QN 
READ1 ,WL,BK,ZZ,B0,B1,B2 
READ1 ,B3,BS,H0,HX,HY,HZ 
READ1 ,HS,HT,HW,QQ,W,RF 
READ1 ,SC,YY,C,TS,SN,DB 
READ1 ,CE,SH,SD,TT,SK,RB 
READ1 ,TP,D1,FE,RD,RT,COILS 
READ1 ,T2,RR,SNL,WF,PE,SN1 
READ1 ,DW1,BP1,BP2,TP1,TP2,ALP 
READ1 ,DR,WR,DY2,TY2,TY,ALY 
READ1 ,HM,SM,AS,AC,ET,EB 
READ1 ,S,FK,FR,RC,RG,RP 
READ1 ,G2,DG2,ALG 2,AL1 
READ1 ,GE,CS,CF,FS,EC,DF 
READ1 ,RS,GC,PT,C1,CW,CP 
READ1 ,EL,CM,CQ,DV/,CC,PBA 
READ1 ,FF,CX,PC,EK,EW,TM 
READ1 ,A,WI,WC,WN,WT,WQ 
READ1 ,TG,FQ,BC1,BT1,BG1,FG 
READ1 ,P1,P2,P3,P4,P5,P7 
READ1 ,A2,AY2,AY,BP,BY,BY2 
READ1 ,BG2,FG2,FNL,AINL,CD,EPNL 
READ1,AP,PL,PLT,FC,FT,P8 
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READ 1,PL8,ALZ 
XR=.07 Q7*A*DF/(BG1*C1) 

XL«XR*(PC+EW+ALZ) 

XD=90.*EC*PI*CM*DF/(2,*PX*(FG+FG2)) 
XQ=CQ*XD/(CM*C1) 

XA=XL+XD 

XB=XL+XQ 

VR=3.1416*DR*RPM/12. 

AGE=((FG+FG2)/(FG))*GE 

PEE=(P1+P2+P3+P4)*PX+P5+P8 

ALF=PEE/CL 

ALA=6.38*DI/(PX*AGE) 

XF=XD*(1.-(Cl/CM)/(2.*CP+ALF*1.273/ALA)) 

S|=PT*PT*PEE*1.E-8*C0ILS 

XU*=XL+XF 

XS».88*XU 

XX=XS 

XY*»XB 

XN“.5*(XX+XY) 

IF(W)414,415,414 
415 X0»0. 

GO TO 422 

414 IF(CS-1.)417,418,417 
418 AKX-1.. 

AKX1*»1. 

GO TO 419 

417 AA=(3.*YY/(PN*QN)) 

AKX«*AA-2. 

IF(AA/3.-.667)420,420,421 





420 AKX1=.75*AA-.25 
GO TO 419 

421 AKX1=.7 5*AA+.25 

419 ABL=(AKX/(CF**2))*.07*ALA 
X0=AKX*(ABL+PC)/AKX1 

XO=XR*(X0+(1,667*(HX+3.*HZ))/(PN*QN*CF**2*DF**2*BS)+.2*EW) 

422 TOS1/(FK*COILS) 

RA=PN*PI*PI*RG/(VA*1000.) 

T A=XN/(628.32*F*RA) 

T5=XS*TC/XA 

T4=2./F 

F SC=XA*.02*(FG+FG2) 

SCR=FNL/FSC 

PRINT3,AC,A,S,XR,HM,XL,RG,XD,RP,XQ,ET,XA,EB,XB,PC,XF,EW,SI,WC,XU 

PRINT3.WI,XS,TP,XX,WR,XY,VR,XN,AS,X0,FK,TC,FR,TA,RC,T5,P1,T4 

PRINT3,P2,TG,P3,FQ,P4,BG1,P5,BT1,P7,BC1,P8,FT,FSC,FC 

PRINT070,SCR,FG 

PUNCH1 f BO,GC,PI,PN,EP,ET 

PUNCH1 # EB,SC,C,XB,XD,PF 

PUNCH1,EE,TG,BT1,FQ f BCl,FG 

PUNCH1,HS,DU,HC,PX,P1,P2 

PUNCH1,P3,P4,P5,P7,BP,ALP 

PUNCH1,BG2,A2,G2,AY2,AY,BY2 

PUNCH1 > DY2,BY,ALY,PT,AS,FK 

PUNCH1 t XA f RG,WF,WQ t V/T f WN 

PUNCH1,SNL,POL,RP,FR,PL,PLT 

PUNCH1,FC,AINL,CD,EPNL,FNL,FT 

PUNCH1,C0ILS,AP,P8,PL8,CK 

PAUSE 
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END 
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PASS 7 OUTSIDE COIL LUNDELL 
DIMENSION GB(4),AE(4),DX(4) 

1 FORMAT(Ell.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 
READ1 ,BO,GC,PI,PN,EP,ET 
READ1 ,EB,SC,C,XB,XD,PF 
READ1 ,EE,TG,BT1,FQ,BC1,FG 
READ1 ,HS,DU,HC,PX,P1,P2 
READ1 ,P3,P4,P5,P7,BP,ALP 
READ! ,BG2,A2,G2,AY2,AY,BY2 
READ1 ,DY2,BY,ALY,PT,AS,FK 
READ1 ,XA,RG,WF,WQ,WT,WN 
READ1,SNL,POL,RP,FR,PL,PLT 
READ1 .FC.AINL.CD.EPNL.FNL.FT 
READ 1,COILS,AP.P8.PL8.CK 
AXX=BO/GC 

IF(AXX-1.)964,965,964 
965 AKSC-2.6 
GO TO 957 

964 IF(AXX-3.75)955,955,956 

955 AKSC=10.**.178/((AXX-1.)**.334) 

GO T0957 

956 AKSC=10.**.11/((AXX-1.)**.174) 

957 XX1=PI*PI*PN 
XX3=3.*EP*PI*PF 
XX2=(ET+EB)/2.-1. 

XX4-AKSC*PI*SC/(C*FG) 

GB(1)“1. 

GB(2)**1.5 
GB(3)-2. 


MAf-44 



GB (4)*=P0L 

AN=ATAN(SQRT(1,-PF*PF)/PF) 
ANI-SIN(AN) 

DO 777 K=1,4 
YB=GB(K) 

AA =ATAN((AN1+XB*YB/100.)/PF) 

AE(K)»COS(AA-AN)+XA*SIN(AA)*YB/100. 
777 DX(K)=.93*XD*YB*SIN(AA)/100. 

PUNCH1,AE(1),AE(2),AE(3),AE(4) 
PUNCHl,DX(1) f DX(2),DX(3),DX(4) 
PUNCH1,B0,GC,PI,PN,EP,ET 
PUNCH1,EB,SC,C,XB,XD,PF 
PUNCH1 f EE,TG,BT1,F0,BC1,FG 
PUNCH1,HS t DU,HC,PX,P1,P2 
PUNCH1 f P3,P4 # P5,P7,BP,ALP 
PUNCH1,BG2,A2,G2,AY2,AY,BY2 
PUNCH1,DY2,BY,ALY,PT,AS f FK 
PUNCH1,XA,RG,WF,WQ,WT,WN 
PUNCH1,XX1,XX2,XX3,XX4,S NL,POL 
PUNCH1 # RP,FR,PL,PLT,FC,FT 
PUNCH1,AINL,CD,EPNL,FNL,C0ILS,AP 
PUNCH1,P8,PL8,CK 
PAUSE 
END 
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PASS 8 OUTSIDE COIL LUNDELL 
DIMENSION Al(90) 

DIMENSION AF.(4),DX(4),BPL(4),BYL(4),BY2L(4),BG2L(4),PLL(4),PTLL(4) 
DIMENSION FFL(4),AlFL(4),CDD(4),EPFL(4) 

1 FORMAT(El 1.5,El 1.5,El 1.5,El 1.5,El 1.5,El 1.5) 

888 FORMAT(FIO.O,FIO.O,FIO.O,FIO.O,FIO.O,FIO.O) 

K=1 

823 READ888,Al (K),AI (K+1),AI (K+2),AI (K+3),AI (K+4),AI (K+5) 

K=K+6 

IF(K-89)823,824,824 

824 READ1 ,AE(1),AE(2),AE(3),AE(4) 

READ! ,DX(1),DX(2),DX(3),DX(4) 

READ1 ,BO,GC,PI,PN,EP,ET 
READ1 ,EB,SC,C,XB,XD,PF 
READ1 ,EE,TG,BT1,F0,BC1,FG 
READl ,HS,DU,HC,PX,P1,P2 
READ1 ,P3,P4,P5,P7,BP,ALP 
READl ,BG2,A2,G2,AY2,AY,BY2 
READl ,DY2,BY,ALY,PT,AS,FK 
READl ,XA,RG,WF,WQ,WT,WN 
READl,XXI,XX2,XX3,XX4,SNL,POL 
READl,RP,FR,PL,PLT,FC,FT 

READl ,AINL,CD,EPNL,FNL,COILS,AP 

READ 1,P8,PL8,CK 

L0AD=2 

DO 900 J-1,4 
AED=AE(J) 

AA=AED*FG+(1.+PF)*FT+FC 
PLL(J)=PL*AA /(FG+FT+FC) 
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P8LL=PLL(J)*PL8/PL 
IF(PF-.195)880,880,881 

881 PR-FQ*CK 
GO TO 882 

880 PR=FQ*(AED-DX(J)) 

882 PTLL(J)=PR+(PLL(J)*2./PX) 
X=PTLL(J)/AP 

BPL(J)=X 
NA=31 
K=1 

GO TO 802 
841 FPL= AT* ALP 

PL7L=P7*.001*AA 

PG2L=(PTLL(J)*PX/2.)+PL7L+P8LL 
BG2L(J)=PG2L/A2 
FG2L=BG2L(J)*G2/.00319 
IF(COILS-1.)855,851,855 

855 BYL(J)=PG2L/AY 
859 NA-61 

X=BYL(J) 

K=3 

GO TO 802 
843 FYL=ALY*AT 

IF(CO ILS-1.)856,857,856 

856 PL5L=P5*.001*(AA+FPL+FYL+FG2L) 

851 IF(C0ILS-1.)852,853,852 

852 X=(PG2L+PL5L)/AY2 
GO TO 854 

853 X=PG2L/AY2 
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854 3Y2L(J)=X 
NA*=61 
K=2 

GO TO 802 

842 FY2L=AT*(DU-[)Y2)/2. 

IF(C0ILS-1.)857,858,857 

858 PL5L=P5*.002*(AA+FPL+FY2L+FG2L) 

BYL(J)=(PG2L+PL5L)/AY 
GO TO 859 

857 FFL(J) = (FG2L+FYL+FPL+fY2L+AA)*2. 
AIFL(J)=FFL(J)/(PT*COILS) 

CDD(J)=AIFL(J)/AS 
900 EPFL(J)=AIFL(J)*FR*COILS 
837 JA=JA/3 

PUNCH 860,JA 

860 FORMAT (13) 

IF(JA)861,862,861 

861 DO 863 J=1,JA 

PUNCH1,PLL(J),PTLL(J),BY2L(J),BYL(J),BG2L(J),BPL(J) 
863 PUNCH1,FFL(J),AIFL(J),CDD(J),EPFL(J) 

862 PUNCH1,EE,TG,BT1,FQ,BC1,FG 
PUNCH1,HS,DU,HC,PX,P1,P2 
PUNCH1,P3,P4,P5,P7,BP,ALP 
PUNCH1,BG2,A2,G2,AY2,AY,BY2 
PUNCH1,DY2,BY,ALY,PT,AS,FK 
PUNCH1 , XA,RG,Y/F,WQ,WT ,WN 
PUNCH1,XX1,XX2,XX3,XX4,SNL,P0L 
PUNCH1,RP,FR ,PF,CO ILS,AP,P8 
PUNCH1,AINL,CD,EPNL,FNL,PL,PLT 
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PAUSE 

802 IF(AI (NA)-X)830,831,031 
831 NA-NA+3 

835 IF(AI(UA)-X)833,834,834 

833 NA=NA+2 
GO TO 835 

834 AX=A1(NA) 

BB1=AI(NA-2) 

DC=AI(NA+1) 

D=AI(NA-1) 

XX= (AX-EJB1) / ( .4343*(LOG (DC)-LOG (D+.0001))) 
Y=AX-X X*.43 43 *L0G(DC) 

AT=EXP(2.306*(X-Y)/XX) 

GO TO ( 83 8, 83 9),LOAD 

838 GO TO (806,807,808,809,810),K 

839 JA=JA+1 

GO TO (841,842,843),K 
830 GO TO ( 836 ,837),LOAD 

836 PRINT 850, 

850 FORMAT (17HMACHINE SATURATED) 

PAUSE 

END 
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C PASS 9 OUTSIDE COIL LUNDELL 

DIMENSION BPL(4),BG2L(4),FFL(4),AIFL(4),CDD(4),EPFL(4),FCUL(4) 
DIMENSION WNL(4),STTL(4),SCUL(4),EDDL(4),T0TL(4),PEFF(4),GB(4) 
DIMENSION PLL(4),PTLL(4),BYL(4),BY2L(4) 

961 FORMAT(F11.3,8X F11.3,F11.3,F11.3,F11.3) 

860 FORMAT (13) 

1 FORMAT (Ell.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 

DO 705 N=1,4 

PLL(N)=0 

ptllOO-o 

BPL(N)=0 
BG2L(N)=0 
BYL(M)-0 
BY2L(N)=0 
FFL(N)=0 
AI FL ( N) = 0 
CDD(N)=0 
EPFL(N)«0 
FCUL(N)=0 
WNL( N)=0 
STTL(N)=0 
SCUL(N)=0 
EDDL(lO B 0 
T0TL(N)=0 
705 PEFF(N)=0 
READ860, JA 
IF (JA) 702,703,702 
702 DO 704 J=1,JA 

READ1 ,PLL(J),PTLL(J),BY2L(J),BYL(J),BG2L(J),BPL(J) 
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704 READ1 

,ffl(j),aifl(j),cdd(j),epfl(j) 

703 READ1 

,EE,TG,BT1,FQ,BC1,FG 

READ1 

,HS,DU,HC,PX,P1 ,P2 

READ1 

,P3,P4,P5,P7,BP,ALP 

READ1 

,BG2,A2,G2,AY2,AY,BY2 

READ1 

,DY2,BY,ALY,PT,AS,FK 

READ1 

,XA,RG,WF,WO,WT,W N 

READ1 

,XX1,XX2,XX3,XX4,SNL,P0L 

READ1 

,RP,FR ,PF,COILS,AP,P8 

READ 1 

,AINL,CD,EPNL,FNL,PL,PLT 

IF(SNL)707,706,707 


707 PUNCH 1,TG,BT1,FQ,BC1,FG,HS 
PUNCH1,DU,HC,PX,P1,P2,P3 
PUNCH1,P4,P5,P7,AP,ALP,A2 
PUNCH1,G2,AY,AY2,DY2,ALY,COILS 
PUNCH1,EE,P8 

706 FEL=AINL*AINL*FK 
T L=F E L+V/T +WQ+W N+WF 
ABX=0 

IF(JA)714,712,714 
714 IF(JA-4) 708,709,700 

709 IF (POL) 708,71 0,708 

710 JA=JA-1 

708 GB(1)=1. 

GB(2)=1.5 
GB(3)=2. 

GB(4)=P0L 

DO 711 K=1,JA 
YB=GB(K) 
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FCUL(K)=AIFL(K)**2*FR*C0ILS 

STTL (K) - ((. 0027 *XA*YB) **2*2.+1.) *WT 

WNL (K) *= ((XX4*YB) **2+1.) *WN 

SCUL(K)=XX1*RP*YB 

ED0L(K)=SCUL(K)*XX2 

T0TL(K)-EDDL(K)+SCUL(K)4WNL(K)+STTL(K)+FCUL(K)+WQ+WF 

711 PEFF(K)=XX3*YB*100./(XX3*YB+T0TL(K)) 

712 IF(POL)958,959,958 

958 PRINT96l,PL,PLL(1),PLL(2),PLL(3),PLL(4) 

PRINT961,PLT,PTLL(1),PTLL(2) f PTLL(3),PTLL(4) 

PRINT961,BP,BPL(1),BPL(2),BPL(3),BPL(4) 

PR I NT961,BG2,BG2L(1),BG2L(2),BG2L(3),BG2L(4) 

PR I NT961,BY2,BY2L(1),BY2L(2),BY2L(3),BY2L(4) 

PR I NT961, BY,BYL(1),BYL(2),BYL(3),BYL(4) 
PRINT961,FNL,FFL(1),FFL(2),FFL(3),FFL(4) 

PRINT961,AINL,AIFL(1),AIFL(2),AIFL(3),AIFL(4) 

PR I NT961,CD,CDD(1),CDD(2),CDO(3),CDD(4) 

PR INT961,EPNL,EPFL(1),EPFL(2),EPFL(3),EPFL(4) 

PR I NT961,WQ,WQ,WQ,WQ,WQ 

PR I NT961,WT,STTL(1),STTL(2),STTL(3),STTL(4) 

PRINT961,ABX,SCUL(1),SCUL(2),SCUL(3),SCUL(4) 

PRINT961,ABX,EDDL(1),EDDL(2),EDDL(3),EDDL(4) 

PR I NT961 ,WN,WNL( 1) 9 V7NL (2) ,WNL(3) ,WNL(4) 
PRINT961,FEL,FCUL(1),FCUL(2),FCUL(3),FCUL(4) 

PR I NT961,WF,WF,WF,WF,WF 

PRINT961,TL,TOTL(1),TOTL(2),TOTL(3),TOTL(4) 

PR I NT961,ABX,PEFF(1),PEFF(2),PEFF(3),PEFF(4) 

PAUSE 

959 PRINT961,PL,PLL(1),PLL(2),PLL(3) 
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PR INT961,PLT,PTLL(1),PTLL(2),PTLL(3) 

PRINT961,BP,BPL(1),BPL(2),BPL(3) 

PRINT96l,BG2,BG2L(1),BG2L(2),BG2L(3) 

PRINT961,BY2,BY2L(1),BY2L(2),BY2L(3) 

PRINT961, BY,BYL(1),BYL(2),BYL(3) 

PRINT961,FNL,FFL(1),FFL(2),FFL(3) 

PRINT961,AIUL,AIFL(1),A!FL(2),AIFL(3) 

PRINT961,CD,CDD(1),CDD(2),CDD(3) 

PRINT961,EPNL,EPFL(1),EPFL(2),EPFL(3) 

PRINT961,WQ,WQ,WQ,WQ 

PRINT961,WT,STTL(1),STTL(2),STTL(3) 

PRINT961,ABX,SCUL(1),SCUL(2),SCUL(3) 

PRINT961,ABX,EDDL(1),EDDL(2),EDDL(3) 

PR I NT961,WN,WNL(1),WNL(2),WNL(3) 

PRINT961,FEL,FCUL(1),FCUL(2),FCUL(3) 

PRINT961,V/F,WF,WF,WF 

PRINT96l,TL,T0TL(1),T0TL(2),T0TL(3) 

PRINT96l,ABX,PEFF(1),PEFF(2),PEFF(3) 

PAUSE 

END 





C PASS 10 OUTSIDE COIL LUNDELL 
DI MENS I ON Al(90) 

1 FORMAT(Ell.5,Ell.5,Ell.5,Ell.5,Ell.5,E11„5) 

888 FORMAT(FI 0.0,FI 0.0,FI 0.0,FI 0.0,FI 0.0,FI 0.0) 

878 FORMAT(F12.5,F12.5,F12.5,F12.5,F12.5,F12.5) 

879 F0RMAT(F1?..5,F12.5,F12.5,F12.5,F12.5,F12.5//) 

K=1 

823 READ888,Al(K),AI(K+1),AI(K+2),AI(K+3),AI(K+4),AI(K+5) 
K=K+6 

IF(K-89)823,824,824 

824 READ1,TG,BT1,FQ,BC1,FG f HS 
READ1 ,DU,HC,PX,P1,P2,P3 
READ1 ,P4,P5,P7,AP,ALP,A2 
READ1 ,G2,AY,AY2,DY2,ALY,COILS 
READ1,EE,P8 

C0REL=3.1416*(DU-HC)/(4.*PX) 

L0AD=1 

YB=.8 

DO 800 U=1,9 

R1=BT1*YB 

R2=FQ*YB 

R3=BC1*YB 

R4=FG*YB 

R5=EE*YB 

X=R1 

NA=1 

K=1 

GO TO 802 
806 FT =HS*AT 
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X=R 3 

K=2 

NA=1 

GO TO 002 

807 FC=COREL Vf AT 
FS=FT+FC 

PL=(2.*(R4+FT)+FC)*(P1+P2+P3+P4)*.001*PX 

PLT=R2+(PL*2./PX) 

BP=PLT/AP 
X=BP 
NA=3 1 
K=3 

GO TO 802 

808 FP=ALP*AT 
PL7=.001*P7*(FC+FT+R4+FP) 

PL8=P8*(2.*R4+2.*FT+FC)*.001 

PG2=PLT*PX/2.+PL7+PL8 

BG2=PG2/A2 

FG2=BG2*G2/.00319 

Z=FG2+FP+R4+FT+FC 

IF(CO ILS—1.)817,810,817 
817 BY=PG2/AY 
821 NA=61 
X=BY 
K=5 

GO TO 802 
810 FY=AT*ALY 

IF(COILS-1.)820,822,820 
820 PL5=.001*P5*(Z+FY’) 





818 IF(CO 1LS—1.0)814,815,814 

814 BY2=(PG2+PL5)/AY2 
GO TO 816 

815 BY2-PG2/AY2 

816 X=BY2 
K=4 
NA=6l 

GO TO 802 

809 FY2=((DU-DY2)/6.)*AT 

IF(CO ILS-1.)822,826,822 
826 PL5=.002*P5*(Z+FY2) 

BY=(PG2+PL5)/AY 
GO TO 821 

822 FNL=(FY+FY2+FG2+FP+R4+FT+FC/2.)*2. 
PRINT878,R5,R1,FT,R3,FC,R4 
PRINT879,PL,PLT,BP,BY2,BY,FNL 
800 YB=YB+.1 
PAUSE 

802 IF(AI(NA)-X)830,831,831 
831 NA=NA+3 

835 IF(AI (NA)-X)833,834,834 

833 NA=NA+2 
GO TO 835 

834 AX=AI(NA) 

BB1=AI(NA—2) 

DC=AI(NA+1) 

D=AI(NA—1) 

XX= (AX-BB1)/(.4343*(L0G(DC)-LOG(D+.00O1))) 
Y=AX-XX*.4343*L0G(DC) 
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AT=EXP(2.306*(X-Y)/XX) 

GO TO (838,839) .LOAD 

838 GO TO (806,807,808,809,810),K 

839 JA=JA+1 

GO TO (841,842),K 
830 GO TO ( 836 , 837 ),LOAD 
836 PRINT 850, 

850 FORMAT (17HMACHINE SATURATED) 
PAUSE 
END 
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SECTION PA 
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HOMOPOLAR INDUCTOR A-C GENERATOR 
COMPUTER PROGRAM AND TEST DATA 



SECTION PA 







6APJ STATOR W.NP.NO | STATOR SLOT _| STATOR STACK 1 PARAMETERS 


HOMOPOLAR COMPUTER DESIGN (INPUT) 


MODEL 


EWO 



REV. A 


PA-01 








































































































































































































SUMMARY OF DESIGN CALCULATIONS - HOMOPOLAR INDUCTOR (OUTPUT) 


MODEL 


EWO 


DESIGN NO. 

i.m’b/u- 


(W) (fl.) 


124) (he) 


dm <7\) 


(27) (T. I/3)| 


(42) (K ,k) 


SOLID CORE LENGTH 


DEPTH BELOW SLOT 


SLOT PITCH 


SLOT PITCH 1/3 DIST. UP 


SKEW FACTOR 


TTcMTTT' 


.89150 


.3?7?9 

.34:192 


1 .00000 


78. poors' 

194.37944 


.o'615"4" 


rjmrj 


CARTER COEFFICIENT 


AIR GAP AREA PER STATOR 


AIR GAP PERM 


EFFECTIVE AIR GAP 


FUND/MAX OF FLD. FLUX 


(67) (K. ) 


( 68 ) ( - ) 


<70 c )(Ao) 


(69) (g * ) 


(71) (Cl ) 


(43) (Kd) 


DIST. FACTOR 


(44) (K p ) 


PITCH FACTOR 


.9$bl3 

. 90436 ' 


.42926 


WINDING CONST. 


(72) (C w ) 


.63W 


POLE CONST. 


(73) (C p ) 


(45) (TT«) 


EFF. CONDUCTORS 


(46) (c c ) 


CONO. AREA 


*>5.25iwr 

TTT75GT 


7.740UT 


END. EXT. ONE TURN 


(48) (Le ) 




DEMAGNETIZING FACTOR 


(74) (C M ) 


(47) (S m ) 


CURRENT DENSITY (STA.) 


(49) (I?» ) 


^K5.?jTQ0(T 


1/2 MEAN TURN LENGTH 


11.05000 


(53) (Rph) 


COLD STA. RES. . 20°C 


oToTTT 


. 47.93 ?- 

MV.bbHW 


CROSS MAGNETIZING FACTOR 


(75) (C.) 


AMP COND/IN 


(128) (A) 


.79807 


REACTANCE FACTOR 


(129) (X) 


(54) (R**) 


HOT STA. RES » X° C 


.01227 


16.12237 

145.50284 

66.59649 


LEAKAGE REACTANCE 


030) (Xg ) 


(55) (EFrop) 


EDDY FACTOR TOP 


,00020 


(56) (EFbot) 


EDDY FACTOR BOT 


1.14130 


REACTANCE OF 
ARMATURE REACTION 


031)(Xod) 


032) (Xoq) 


(62) ( A I ) 


STATOR COND. PERM. 


2.71830 


161.625ZT 


SYN REACT DIRECT AXIS 


033) (X d ) 


(64) (At) 


END PERM. 


H 




82.7T8BTT 


SYN REACT QUAD AXIS 


034) (X q ) 


(65) ( ) 


WT. OF STA COPPER 


26.21200 


(66) ( ) 


WT. OF STA IRON 


(60c* (Pm) 


LEAKAGE PERMEANCE 


58 . 27200 " 

k:mrr 


23 . 1063 (> 

-7^2^ 


FIELD LEAKAGE REACT 


060) (Xf ) 


FIELO SELF INDUCTANCE 


061) (Lf ) 


,00000 


(84u) (P 5 ) 


LEAKAGE PERMEANCE 


(85o) (P 6 ) 


LEAKAGE PERMEANCE 


21.05326 


.00000 


DAMPER 

LEAKAGE REACT 


063) (XD a) 


39.22*17; 


065) (X Dq ) 


UNSAT. TRANS. REACT 


066) (X\iu> 


( 86 a) (P 7 ) 


LEAKAGE PERMEANCE 


i 


Qj ?02 


14.52 12.1 


SAT. TRANS. REACT 


067) (X’ d ) 


(1S3) (a CF ) 


FLD. COND. AREA 


054) (R F) 


COLD FLD RES 20 ° C 


055) (R F ) 


HOT FLD RES » X° C 


056) ( 


WT. OF FLD. COPPER 


,00406 
1.47 900 


34 . 521 00 


SUB. TRANS.REACT DIRECT AX. 


068) (X M j ) 


32.71800 


SUB. TRANS.REACTQUAD AX. 


069) (X” q ) 


2,53 1 9.U 


58.61950 


NEG SEQUENCE REACT 


0 70) (X 2 ) 


11 .328 QQ_ 


2.3 9774 


ZERO SEQUENCE REACT 


(172) (X 0 ) 


057) ( > 


WT. OF ROTOR IRON 


mmi 


>712.50Q3Q 


TOTAL FLUX 


(88) (<£ t ) 


045) (V r ) 


PERIPHERAL SPEED 


076) (T do ) 




632.80000 

47.59600 


FLUX PER POLE 


(92) (d-p) 


GAP DENSITY 


(9S) (B fl ) 


OPEN CIR. TIME CONST. 


,42 57 


97. 916 75 


TOOTH DENSITY 


(91) (B t ) 


o 


(177) (To) 


ARM TIME CONST. 


.01018 


49.18050 


CORE DENSITY 


(94) (B c ) 


(178) (T’ d ) 


TRANS TIME CONST. 


.0 9093 


40.69794 


TOOTH AMPERE TURNS 


(97) (F t ) 


(179) <T" d > 


SUB TRANS TIME CONST. 


.00500 


56671 


CORE AMPERE TURNS 


(98) (F c ) 


-L------ ———1 . 

PERCENT LOAD 




15° | 

1 __ 1 

OPTIONAL 

(0m> (91o) LEAKAGE FLUX 

HDKKERBH 

1 /■'ml ( 202 *) 

IVIRH 

■ I I ■! 1 — 


.noa- 

(fV*)(96o) GAP AMPERE TURNS 

KRociim 

F„l (203) 

IIIIIWHIHI 

wwmxsm* 


.oao 

(B p ) (104h)POLE DENSITY 

■ II Hll 1 — 

Bpl (213b) 

■ 1IH 111 

77.910 

wsMEsmm 

.000 

(B* ) (91c) TOOTH DENSITY 

97.916 I 

Btl (205) 

110.91?. 

■iiiiiim 

fBSEaSMM 

.000 

(B iK ) (113) SHAFT DENSITY 


B.hl 

■ IIIHHIIM 

■mu— 

iiTtnii— 

.000 

(B c ) (94) CORE DENSITY 

wiityH 

B c , (200g) 

m&jvj&mm 


91.5?9 

.000 

(Bvc) (125o)COIL YOKE DENSITY 

■ 1II III— 


■H.H.TOM 


.000 

.000 

(F n |) (127) TOTAL Nl 

inrriiiin 

(Ff| ) (236) 

mimrnm 

iiagwwmi—ii 

.00(5 1 

.000 

(1 f n |) (127o)FIELD AMPS 


(Iffl ) (237) 


■QSKBHH 

.000 

. 000 

(SF ) (127c) CUR.DENS.(FLD.) 

1 T7T47- 

(S ft ) (239) 

■SflUHH 

irrrrnw 

. 000 

• ouo 

(Ef ) 0278) FIELD VOLTS 

uaanuM 

(EMI) (238) 

mamatMmu 


TTJTT0 1 

. UUU 

(l 2 R r )(182) FIELD LOSS I 

m&E&XM i 

mrirgna 


SEtSSUH 


* 000 

(F&W) (183) F&W LOSS 


(F&W) (183) 

mm 

lISEIiSIUSHi 


3 / 0 ./UU 

<W, nl ) (184) STA TOOTH LOSS 



HHiSIilOH 

■ 11 mu— 

. uuu 

. UUU 

(W c ) (185) STA CORE LOSS 

■ llll|i ||||—1 

(Wc) 085) 1 

■ II IHMIIIM 

■TIlTfTP 

^ rouo 

. uuu 

(W Dn 0 (186) POLE FACE LOSS 

■ nFiiri 

mnssM 

■i i mm— 

III 1 W lil 1 — 


. 000 

(W dnl ) (193)OAMPER LOSS 

■■nninH 

|(Wdl|)<2<4> ] 



.ooO 

. 000 

(I 2 R. K194) STATOR CU LOSS 

.000 

fimRPR« 

ffiHHIIM 

mm 11 — 

.000 

.000 

(-) (195) EDDY LOSS 

■mm 

mmssm 

msmsiam 


■HPiatliBHi 

.000 

( _ ) (196) TOTAL LOSSES 

wiumj&mm 

( - ) (247) 

IH—i 

SfECVRISni 


_.uoa. 

( - ) ( r ) RATING (KW) 


( - ) (248) 

maimum 

90.31b 


• UUU 

-n rrn— 

( - ) ( - ) PERCENT EFF. 

■■BJiUtflH 

( - ) (251) 

jjj i 

85". 983" 

. 000 

• UUU 


a 

s 


_J 
<D ( 
< 
5 : 
< 

. > 


PA-03 designer. 


DATE 


REV. A 




































































































































































































































HOMOPOLAR 


NO LOAD SATURATION OUTPUT SHEET 


\lTfMS 

13) (E) 

VOLTS 

(WA)(F g ^ 

AIR CAP A.T. 


m) <p e ) 

CORE a.T. 

(*1> IB f ) 
TOOTH DENSITY 

f*7) (F f ) 

TOOTH A.T. 

VOLTS \ 

POLE DENSITY 

HB9J 

<•»*)' <&**) 
SHAFT DENSITY 

<l«) 

SHAFT a.t. 

IttTdHeyc) 

yO*S DENSITY 

U27» IF,,,) 
TOTAL A.T. (H.L. 


166 .27200 

50.20037 

768.52306 

44.36387 

48.55613 

57.96235 

2.51178 

77.68069 

78.3334C 
71.13030 

4.9913" 

1742.0747 

1 

v<r; 

187.05600 

50.52312 

864.58044 

44.38257 

48.86832 

58.71500 

2.53909 

77.75708 

88.12507 

72.36165 

mm 

■ 

IOO-C. 

207.04000 

50.84580 

960.6 5383 
44.40128 

49.18050 

59.46763 

2.56671 

77.83353 

■RHP 

Wjajm 

- 1 

40.6975 

2199.89541* 

110% 

228.62.400 

51.16864 

1056.71920 

44.42000 

49.49269 

60.22027 
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193.6767b! 

2892.649301 
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270. 19200 
51.81415 
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44.45745 
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61.72556 

2.65137 
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REVOLUTIONS PER MINUTE 


TEST VALUES OF LOSSES FOR A 
12000 RPM, 80 KVA, 400 CPS, 
HOMOPOLAR INDUCTOR A~c’ 
GENERATOR 


12000 


10000 


8000 


6000 



LOSSES IN KILOWATTS 
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TEST EFFICIENCY VALUES FOR A 12000 RPM 
80 KVA, 400 CPS, HOMOPOLAR INDUCTOR A-C 
GENERATOR 



PER-UNIT LOAD (1. 0 P. U. IS 80 KVA) 

















COMPUTED ritCCSdURE FOR 


HOMOPOLAft DSSI iH FALOJLaTIGiE 

Clear core (no switch control). 

Insert output Form #1 into typewriter, set margin for correct 
output, and set typewriter for single space* 

Load pass ff 1 followed by input parameters (output punched cards)# 
Reset and Load pass //2 followed by output from pass #1 (output 
printed plus punched cards)# 

Reset and Load pass #3 followed by output from pass ff2 (output 
printed plus punchedcards)# 

Reset and Load pass ffh followed by output from pass ff3 (output 
printed plus punched cards). 

Reset and Load pass ffS followed by output from pass ff\± (output 
printed plus punched cards). 

Reset and Load pass #6 followed by saturation curve values* and 
output from pass ff$ (output printed and punched cards). 

Reset and Load pass #6A followed by output from pass #6 (output 
punched cards). 

Reset and Load pass #7 followed is/ output from pass y6A (output 
punched cards). 

Reset and Load pass #8 followed by saturation curve va lues* and 
.output from pass #7 (output printed and punched cards). 



12 . Reset and Load pass #9 followed by output from pass #8 (output 
printed plus punched cards if no load saturation curve required). 

13 . If there is punched card output from pass #9, a no load saturation 
curve is required. Insert output Form #2 into typewriter and reset 
margin. Load pass #10 followed by saturation curve values* and 
output from pass #9 (output printed). 


* gjturation curve values are loaded in order shown on Input Form #1. 
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ALL INPUT PARAMETERS ARE IN FORMAT F7.0 (FIG. 1) 


1. 

0 0 8 0 0 9 eh 


III4II 


I 1 I t | 1 11 
22222222 
33333313 
44444444 
55555555 
6(6(66(6 
777777777 
«is 8 mu 


9 9 9 9 S 9 S|S 9 

i ? j * i i if 


10 . 

o e o o g ofo 

• • Ifl N 12 U 

11111 


00000 ole 10008 ofo 

Mia a in •»» n m n it 77 


I 11 


2222222 

3133333 

4 4 4 4 4 4 4 

5 5 5 5 5 5 5 
8(66506 

7 7 7 7 7 7 

mm 


t|i 111111«1111 
222227222222 Al 


9 9 9 9 9(9 

• »• tl « ?] 

»Otl 


100 . 


5 5 5 5 5 5 5 


.001 


o o o o o o|o 

31 33 U 34 

I 111 11|| 


733337337333)3 

44444444444444 


5 5 5 5 5 5 5 


) 


6 ( 666(66666666 

17777777777777 

mmmstm 


9 ? 3 9 9 f9 9 0 9 9 9 M9 

K n !• 1 * » *2 n 24 ft 2 € v ma 


.1 


0 U 0 0 0 Gjo 
1 1 I 1 I 111 


2222222 
3333333 
4444444 
5 5 5 5 5 5 5 
66(66(6 
7 7 7 7 7 7 7 
mom 


9 9 9 3 ? 

M 31 33 33 m 


.01 


0 0 0 0 0 ®jj#80060 0|0 

da « n n 4i it n u n H a 


2222222 
3333333 
44 4 4 4 44 

5 5 5 5 5 5 5 

6 6 6 6 6 H6 
7; 7 7 7 7| 


Ill III I 


9|3 9 9 9 9 9 919 9 


10 . 


11 I 1 I ill I 1 I 1 1 
2 2 2 2 2 Vp2 2 217 2h 


5555555 
66(6(66 
7 7 7 7 7 7 7 

imm 


1.0 


oooooioooooooh 

ifiasin.iicMBKaaiinfr 
l|l 1 1 1 1 I f 1 1 I 1 1 1 


3 3 3 3 3 3|3 3 3 3 3 3 3b 

44444 4|4 44444 4|4 


5555555 
6 ( 6 ( 6 6 6 
7 7 7 7 7 7 7 

mini 


f»uhJi40ii «au 


FIG. 1 


9 9 9 9 9|9 9 9 9 9 9 9B 

44 43 4| 47 4« 4^0 il S* » 5$| 


4 4|05 


1000. 


4 4 4 4 4 4 I 4 


7 7 7 7 7 717 

HI 11 Jo 


9 9 9 9 9 » 




10 . 


2 2 2 2 2 V 2 2 22 22\2 

3 3 3 3 3 3|3 3 3 3 3 3 3b 
4 4 4 4 4 4J4 

5 5 5 5 5 5|5 5 5 5 5 5 5|5 

6 6 6 ( 6 (|( ( ( 6 6 6 6k 
7 7 7 7 7 7 7 
III III I 


959993b 

edtiMM wfc 


000000000 
nnitnannaM 
1 1 1 1 t I 1 1 1 


221222222 

333333333 

444444444 

555555555 

6(666(1(6 

777777777 

III 111111 


999999999 

n 23 * nx 7j mm 


ALL SATURATION CURVE VALUES ARE IN FORMAT F10.0 (FIG 2) 
(ALL SATURATION CURVES MUST HAVE 5 CARDS) 
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HOMDPOLAR INDUCTOR 


CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 


A, a 


(U6) 

A c 

AC 

(68) 

A g 

GA 

(79) 

A 

P 

AP 

(128) 

A 

A 

(91b) 

AT 

ATH 

(9Ua) 

A c 

ACR 

(112) 

^sh 

ASH 

(12Ua) 

A 

y 

AY 

(121+b) 

V 

AYC 

(12Uc) 

V 

AYR 

(UU) 


AB 

(153) 


AS 




(15) 

W 

BV 

(20) 

B 

BK 

(22) 

b o 

BO 

(22) 

b l 

Bl 

(22) 

b 2 

B2 

(22) 

b 3 

B3 

(22) 

b s 

BS 

(57) 

^tm 

TM 

(57a) 

b t 1/3 

SM 

(76) 

hp 

BP 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(78) 

^coil 

BCOIL 

(91) 

B t 

BT 

(91c) 

B 't 

BT 

(94) 

B c 

BC 

(95) 

B g 

BG 

(104b) 

B 

P 

BP 

(113) 

B sh 

BSH 

(125a) 

B yc 

BYC 

(126a) 

B y 

BY 

(135) 

bbo 

WO 

(135) 

bbl 

B 

(200b) 

®PL 

BPL 

(200g) 

B CL 

BCL 

(202c) 

B 'shL 

BSHL 

(205) 

b tl 

BTL 

(213b) 

b pl 

BPL 

(215a) 

b shl 

BSHL 

(228a) 

%CL 

BYCL 

(229a) 

®YL 

BYL 


C, c 


(32) 

C 

C 

(60) 

C x 

CX 

(71) 

C 1 

Cl 
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CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 

(72) 


CW 

(73) 

C P 

CP 

(7h) 


CM 

(7$) 

°q 

CQ 


D, d 


(11) 

d 

DI 

(11a) 


DR 

(12) 

D 

DU 

(35) 


DB 

(78) 

‘Vc 

DYC 

(78) 

d coil 

PCOIL 

(76) 

^coil 

DCOIL 

(78a) 

d, sh 

DISH 

(78a) 

^sh 

DSH 


E, e 


(3) 

E 

EE 

(U) 

®ph 

EP 

(55) 

^top 

ET 

(56) 

E F B ot 

EB 

(127) 

E f 

EF 

(198) 

e d 

ED 

(238) 

*Wl 

EF 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 


F, f 


(5a) 

f 

F 

(96) 


FH 

(96a) 

V 

FGM 

(97) 

f t 

FT 

(98) 

F 

c 

FC 

(106a) 

f p 

FP 

(nil) 

f sh 

FSH 

(125b) 

F 

yc 

FYC 

(125d) 

F 

r yr 

FYR 

(126b) 

F 

y 

FY 

(127) 

f nl 

FNL 

(180) 

F 

sc 

FSC 

(183) 

F & W 

WF 

(198b) 

F dm 

FOIL 

(199) 

F *L 

FGL 


g 

(200) 

f 'tl 

FTL 

(200c) 

f ‘pl 

FPL 

(201) 

f cl 

FCL 

(202d) 

f shl 

FSHL 

(203) 

f gl 

FGL 

(206) 

f tl 

FTL 

(213L) 

f pl 

FPL 

(216a) 

f shl 

FSHL 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(228b) 

F ycL 

FYCL 

(228d) 

F yrL 

FYRL 

( 229 b) 

F yL 

FYL 

(236) 

f fl 

FFL 


JkJSL 


(59) 

g 

GC 

(69) 

g e 

GE 


H, h 


(22) 

h Q 

HO 

(22) 

h i 

HX 

(22) 

h 2 

HY 

(22) 

h 3 

HZ 

(22) 

^8 

HS 

(22) 

\ 

HT 

(22) 

K 

HW 

(21*) 

h c 

HC 

(37) 

^st 

SH 

(38) 

h 'st 

SD 

(76) 

h p 

HP 

(76) 

h‘ 

P 

HP1 

(135) 

‘So 

HD 

(135) 

*b 

DD 

(137) 

‘Si 

H 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 


hi 


(8) 

^■PH 

PI 

(127a) 

X FNL 

FI 

(182) 

l2R F 

PR 

(237) 

t ffl 

FI 

(2U1) 

i 2 r f 

PR 


K, k 


(2) 

K VA 

VA 

(9a) 

K c 

CK 

(16) 

Ki 

RK 

(19) 

k 

WL 

m 

k sk 

FS 

(U3) 

K d 

DF 

(UU) 

s 

CF 

(61) 

K 

X 

FF 

(63) 

Kg 

EK 

(67) 

K s 

CC 


- hz l 


(13) 

1 

L 

(17) 

X s 

SS 

(36) 

1 e2 

CE 

(U8) 

h 

EL 

(U9) 


HM 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(76) 

S 

PL 

(139) 

lb 

SB 

(11*7) 

Hr 

FE 

(161) 

Lp 

SI 


M, m 


(5) 

ID 

PN 


N, n 


04) 


HV 

(30) 

n s 

SC 

(34) 

N st 

SN 

(34a) 

N» 4. 
st 

SNL 

(45) 


EC 

(138) 


BN 

(146) 

% 

PT 


0, o 



P. P 


(6) 

P 

PX 

(9) 

P F 

PF 

(80c) 

P m 

PM 

(84a) 

P 5 

P5 

(85a) 

P 6 

P6 

(86a) 

P 7 

P7 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 


a 


(23) 

Q 

QQ 

(25) 

q 

QN 


R, r 


(7) 

r pm 

RPM 

(53) 

^SPH 

RG 

(5U) 

RgPH (HOT) 

RP 

(15U) 

R f (COLD) 

FK 

(155) 

R f (HOT) 

FR 


S, s 


(127c) 

Sp 

CDD 

(181) 

S CR 

SCR 

(239) 

S FL 

CDD 


... h t 


(78) 

Vr 

TYR 

(78) 

^yc 

TYE 

(78) 

V 

TY 

(176) 

T 'do 

TC 

(177) 

T a 

TA 

(178) 

T, d 

T5 

(179) 

d 

Ti* 
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CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 

* 

&JL— 



I* 1 — 


(11*5) 

V 

r 

VR 




(181*) 

WTNL 

ST 

(185) 

W C 

WQ 

(186) 

W NPL 

WN 

(193) 

W DNL 

WD 

(21*2) 

W TFL 

ST 

(21*3) 

V* 

PP 


PFL 


(21*1*) 

Wdfl 

DL 


X, x 


(50) 

X 

CO 

o 

o 

T1 

(129) 

X 

XR 

(130) 

X L 

XL 

(131) 

X ad 

ID 

(132) 

*aq 

XQ 

(133) 

*d 

XA 

(131*) 

X q 

XB 

(11*2) 

x D °c 

T3 

(150) 

x c ° c 

T2 
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CALCULATION 

NUM3ER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 

(160) 

*T 

IF 

(163) 

X Dd 

XI 

(165) 

X Dq 

X2 

(166) 

X ' D u 

XU 

(167) 

X ’d 

xs 

(168) 

X "d 

XX 

(169) 

x »», 

q 

XY 

(170) 

X 2 

XN 

(172) 

x o 

XO 


-LlJL- 


(31) 

y 

YY 



T 


(26) 

~r 8 

TS 

(27) 

T S l/3 

TT 

(UO) 

^sk 

SK 

(Ul) 

Tp 

TP 

(1U0) 

T b 

TB 

(62) 

A. 

Ai 

PC 

(6k) 

Ae 

EW 

(70c) 

Aa 

AG 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(88) 

-4- 

h 

TG 

(91a) 

K 

PHM 

( 92 ) 

% 

FQ 

( 99 ) 

h 

PH 7 

(112a) 

fei 

PSH 

(118a) 

% 

PH5 

(121a) 

*6 

PH6 

(198c) 

^ mL 

PML 

(200a) 

hi 

PPL 

(200d) 

^5L 

PH 5 L 

(200e) 

hi 

PH6L 

(200f) 

0 CL 

PCL 

(202) 


PH 7 L 

(202b) 

^SHL 

PSHL 

(202e) 

AjL 

PML 

(207a) 

^ 7 L 

PH 7 L 

(213) 

hi 

PPL 

(2U*a) 

^SHL 

PSHL 

(220a) 

hi 

PH6L 

(226a) 

hi 

PH 5 L 




(51) 

4 

RS 

(1U1) 

/ D 

RE 

(15D 

/ F 

RR 
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CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 




(77) 

<< 

PE 




(198a) 

■Or 

AN 


K 


(187) 

h 

D1 

(188) 

*2 

I 

(189) 

K 3 

D3 

(190) 

\ 

DU 

(191) 

s 

D5 

(192) 

h 

D6 
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C PASS 1 HOMOPOLAR INDUCTOR 

DIMENSION DA(8),DX(6),DY(8),DZ(8) 

1 FORMAT(El 1.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 

2 F0RMAT(F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7.0) 
33 READ2,VA,EE,EP,PN,F,PX,RPM,PI,PF,CK 

READ2,POL,DI,DU,CL,HV,BV,SF,WL,BK,ZZ 

READ2,BO,Bl,B2,B3,BS,HO,HX,HY,HZ,HS 

READ2,HT,HW,QQ,W,RF,SC,YY,C,DW,SN 

READ2,SN1,DW1,DB,CE,SH,SD,PBA,SK,T1,RS 

READ2 # GC,GP,Cl,CW,CP,EL,CM,CO,BP,PL 

READ2,HP,HP1,PE,DR,RK,WR,D1,W0,HD,DD 

READ2,H,B,BN,SB,TB,RE,T3,DSH,DISH,DISH1 

READ2,ALH,TYPY,TY,TYE,TYR,DYC,PCOIL,DCOIL,BCOIL,PT 

READ2,FE,RD,RT,T2,RR,SNL,WF,PM,P5,P6 

READ2,P7 

SS«SF*(CL—HV*BV) 

HC—(DU—Dl—2.0*HS)*0.5 
ZY»0.7*HS 
IF(HC-ZY) 33,33,5 

5 QN**QQ/(PX*PN) 

TS-3.142*DI/QO 
IF(ZZ-4.0)29,30,29 

29 TT-(0.667*HS+DI)*3.142/QQ 
GO TO 31 

30 TT-3.1416*(DI+2.*H0+1.32*BS)/QQ 

31 IF(ZZ-1.0)6,6,7 

6 BO-BS 

CC-(5.0*GC+BS)*TS/((5.0*GC+BS)*TS-BS*BS) 

GO TO 8 

PAr.19 
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7 QO(4.44*GC+0.75*B0)*TS 
COQC/(OC-BO*BO) 

8 CS«YY/(PN*0_N) 

TP«3.142*DI/PX 

IF(SK)32,32,92 

32 FS-1.0 
GO TO 34 

92 FS=SIN(1.571*SK/TP)*TP/(1.571*SK) 

34 IF (PBA-60.)9,9,1 0 

9 D-1.0 

GO TO 95 

10 D-2.0 
95 1 -QN 

U-l 

IF(QN-U)36,36,35 

35 U=PX*PN 
XX-U 
N-U 

DO 11 K=1,N 

Z-U/XX 

l-Z 

Z1-I 

IF(Z-Z1)12,12,11 
12 ZY-QQ/XX 
l-ZY 
Z1-I 

IF(ZY-Z1)37,37,11 

11 XX-XX-1. 

36 ZY-QN 
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37 DF-SIN(1.571*D/PN)/(ZY*D*SIN(1.571/(PN*ZY))) 
CF=SIN(YY*1.571/(PN*QN)) 

EC=QQ*SC*CF*FS/C 

DT-DW1 

IF(DT) 13,13,14 

13 AO0.785*DW*DW*SN1 
GO TO 24 

14 ZY=0.0 
DA( 1 )“0.05 
DA(2)«0.072 
DA(3 )**0.125 
DA(4)-0.165 
DA(5)-0.225 
DA(6)»0.438 
DA(7)“0.688 
DA(8)-1.5 
DX(1)=0.000124 
DX(2)«0.00021 
DX(3)“0.00021 
DX(4)=0.00084 
DX(5)“0.00189 
0X(6)-0.00189 
DY(1)-0.000124 
DY(2)“0.000124 
DY(3)“0.00084 
DY(4)“0.00084 
DY(5)“0.00189 
DY(6)-0.00335 
DY(7)“0.00754 
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DY(8)-0.03020 
DZ(1)»0.000124 
DZ(2)«0.000124 
DZ(3)“0.000124 
DZ(4)»0.00335 
DZ(5)«*0.00335 
DZ(6)=0.00754 
DZ(7)»0.0134 
DZ(8)«*0.03 02 

93 IF(DT-.05)94,94,15 

15 JA-0 
JB-0 
JOO 
JD-0 

16 JA-JA+1 
JB-JB+1 
JC-JC+1 
JD=JD+1 

IF(DT-DA(JA))17,17,16 

94 D-0 

IF(ZY)23,23,27 

17 IF(DW-0.188)18,18,19 

18 CY-DX(JB-I) 

CZ-DX(JB) 

GO TO 22 

19 IF(DW-0.75)20,20,21 

20 CY-DY(JC-I) 

CZ-DY(JC) 

GO TO 22 
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21 CY-DZ(JD-I) 

CZ-DZ(JD) 

22 D**CY+(CZ—CY)*(OT-DA( JA-1) ) /(DA( JA)—DA( JA—1) ) 
IF(ZY)23,23,27 

23 AC-(DT*DW-D)*SN1 

24 IF(RT)25,25,26 

25 AS»0.785*RD*RD 
GO TO 28 

26 ZY-1.0 
DT-RT 
DW-RD 
GO TO 93 

27 AS-RT*RD-D 

28 S-PI/(C*AC) 

PUNCH1,VA,EE,EP,PN,F,PX 

PUNCH1 f RPH,PI,PF,CK,POL,DI 

PUNCH1,DU,CL > SS,HC,SF,ON 

PUNCH1 # WL,BK,ZZ,BO,B1,B2 

PUNCH1.B3.BS,HO,HX,HY,HZ 

PUNCH1,HS,HT,HW,QQ,W,RF 

PUNCH1,SC,YY,C,TS,SN,DB 

PUNCHI.CE.SH.SD.TT.SK.TI 

PUNCH1,RS,GC,GP,C1,CW,CP 

PUNCHI.EL.CM.CO.BP.PL.HP 

PUNCH1,PE,DR,RK,CC,WR,D1 

PUNCH1,WO,HD,DD,H,B,BN 

PUNCH1,SB,TB,RE,T3,DSH,DISH 

PUNCH1,ALH.TYPY.PCOIL.DCOIL.BCOI L.TY 

PUNCH1,TYE,TYR,DYC,PT,FE,RD 
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PUNCH1,RT,T2,RR,SNL,WF,CS 

PUNCH1,AS,FS f TP,DF,CF,EC 

PUNCH1,AC,S,HP1,DI SHI,PM,P5 

PUNCH 1,P6,P7,PBA 

PAUSE 

END 


PA~2.li 



C PASS 2 HOMOPOLAR INDUCTOR 

1 FORMAT(E11.5,E11.5,E11.5,E11.5,E11.5,E11.5) 
3 FORMAT(9X F12.5,2X F12.5) 


READ 

1,VA,EE,EP,PN,F,PX 


READ 

1,RPM,PI,PF,CK,POL,D 

1 

READ 

1,DU,CL,SS,HC,SF,QN 


READ 

1,WL,BK,ZZ,BO,B1,B2 


READ 

1,B3,BS,HO,HX,HY,HZ 


READ 

1,HS,HT,HW,OQ,W,RF 


READ 

1,SC,YY,C,TS,SN,DB 


READ 

1,CE,SH,SD,TT,SK,T1 


READ 

1,RS,GC,GP,C1,CW,CP 


READ 

1,EL,CM,CQ,BP,PL,HP 


READ 

1,PE,DR,RK,CC,WR,D1 


READ 

1,WO,HD,DD,H,B,BN 


READ 

1,SB,TB,RE,T3,DSH,DI 

SH 

READ 

1,ALH,TYPY,PCO1L,DCO 

1L,BCOIL,TY 

READ 

1,TYE,TYR,DYC,PT,FE, 

RD 

READ 

1,RT,T2,RR,SNL,WF,CS 


READ 

1,AS,FS,TP,DF,CF,EC 


READ 

1,AC,S,HP1,DISHI,PM, 

P5 


READ1 # P6,P7,PBA 
GA-3.142*DI*CL 
AG=«6.38*DI/(PX*GC*CC) 

GE*»CC*GC 
IF(C1) 44,43,44 

43 C1-(0.649*L0G(PE)+1.359)*((GC/GP)**0.352) 

44 IF(CW)45,45,46 

45 CW-0.707*EE*C1*DF/(EP*PN) 

PAr-2 5 


I 



46 TG=6000000.0*EE/(CW*EC*RPM) 

BG«=TG/GA 

IF(CP)47,47,48 

47 CP«(GC/GP)**0.4l*PE*(L0G(GC/TP)*.0378+1.191) 

48 FQ=TG*CP/PX 
IF(ZZ-3.0)49,50,51 

49 SM-TT-BS 
GO TO 53 

50 SM»(3.1416*(DI+2.*HS)/QQ)-B3 
GO TO 53 

51 IF(ZZ-4.0)50,52,49 

52 SM*=TT-. 94*BS 

53 ! F {EL) 54,54,62 

54 I F(RF) 55,55j,6 1 

55 i F(PX-2.0 ) 56,56,57 

56 0=1.3 

GO TO 60 

57 IF(PX—4,0) 58,58,59 

58 u«1.5 

GO TO 60 

59 U=1.7 

60 EL=3.l42*U*YY*(Di+HS)/0Q+0.5 
GO TO 62 

61 EL«=2.0*CE+(3.142*(0.5*HX+DB))+(YY*TS*TS/(SQRT(TS*TS-BS*BS))) 

62 HM-CL+EL 
RX«RS*0.000001 

RB=*(T 1+234.5)*0.003 94*RX 
IF(SH)37,38,40 
38 ET=*1 
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EB-1 

GO TO 39 

40 AA«0.584+(SN*SN-1.0)*0.0625*(SD*CL/(SH*HM))*2.0 

AB«(SH*SC*F*AC/(BS*RB*1000000.0))**2.0 

\ 

ET «AA*AB*0.003 3 5+1.0 
EB-ET-0.00168*AB 
39 IF(CM)63,63,64 

63 AA-SIN(3.142*PE) 

AB-SJN(1.571*PE)*4.0 

CH-(3.142*PE+AA)/AB 

64 PRINT3,SS,CC,HC,GA,TS,AG,TT,GE,FS,C1,DF,CW,CF,CP,EC,EL,AC,CM 
PUNCH1,VA,EE,EP,PN,F,PX 

PUNCH1,RPM,PI,PF,CK,POL,DI 

PUNCH1,DU,CL,SS,HC,SF,QN 

PUNCH1,WL,BK,ZZ,B0,B1,B2 

PUNCH1,B3,BS,H0,HX,HY,HZ 

PUNCH1,HS,HT,HW,QQ,W,GE 

PUNCH 1,SC,YY,C,TS,BG,TG 

PUNCH1,FQ,TT,HM,SH,GA,EC 

PUNCH1,GC,C1,CW,CP,RX,RB 

PUNCH1,EL,CM,CQ,BP,HP,HP1 

PUNCHt,PL,PE,DR,RK,CC,WR 

PUNCH1,D1,WO,HD,DD,H,B 

PUNCH1,BN,SB,TB,RE,T3,DSH 

PUNCH1,DISH,ALH,TYPY,PCOIL,DCOIL,BCOIL 

PUNCH!,TY,TYE,TYR,DYC,PT,FE 

PUNCH1,RD,RT,T2,RR,SNL,WF 

PUNCH1,CS,AS,ET,TP,DF,CF 

PUNCH1,EB,AC,S,AG,SM,DISH1 
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PUNCH1,PM,P5,P6,PM,PBA 

PAUSE 

END 
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C PASS 3 HOMOPOLAR INDUCTOR 

1 FORMAT(E11.5,El 1.5,El 1.5,El 1.5,El 1.5,El 1.5) 
3 FORMAT(9X F12.5.2X F12.5) 

READ1, VA,EE,EP,PN,F,PX 

READ1, RPM,PI,PF,CK,POL,DI 

READ1, DU,CL,SS,HC,SF,QM 

READ 1, WL,BK,ZZ,B0,B1,B2 

READ1, B3,BS,HO,HX,HY,HZ 

READ1, HS,HT,HW,QO,W,GE 

READ1, SC,YY,C,TS,BG,TG 

READ1 ,FQ,TT,HM,SM,GA,EC 

READ 1,GC,C1,CW,CP,RX,RB 

READ 1,EL,CM,CQ,BP,HP,HP1 

READ 1,PL,PE,DR,RK,CC,WR 

READ 1,D1,WO,HD,DD,H,B 

READ 1,BH,SB,TB,RE,T3,DSH 

READ 1,DISH,ALH,TYPY,PCOIL,DCOIL,BCOIL 

READ 1,TY,TYE,TYR,DYC,PT,FE 

READ 1,RD,RT,T 2,RR,S NL,WF 

READ 1,CS,AS,ET,TP,DF,CF 

READ 1,EB,AC,S,AG,SM,DI SHI 

READ1,PM,P5,P6,P7,PBA 

FGML-.45*EC*PI*CM*DF/PX 

A-PI*SC*CF/(C*TS) 

RY“SC*00*HM/(P N*AC*C*C) 

RG«RX*RY*2. 

RP=RB*RY*2. 

1F(CQ)69,69,70 
69 AA-1.571*PE 
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AB«3.1416*PE 

CQ=(0.5*C0S(AA)+AB-SIN(AB))/(4.0*SIN(AA)) 
70 XR=.0707*A*DF/(C1*BG) 

ZX=YY/(PN*QN) 

IF(ZZ-5.)350,351,350 

351 FF=1.0 
GO TO 75 

350 IF(PB A-60. ) 3 52,3 53,3 52 

353 IF(ZX-.667)354,355,355 

355 D=.75 
Z=.25 

GO TO 74 

354 D=1.5 
Z—.25 
GO TO 74 

352 I F (ZX-.667)356,357,357 
357 FF=>.75 

GO TO 75 

356 0=1.2 
Z=-.05 

74 FF=D*ZX+Z 

75 CX=FF/(CF*CF*DF*DF) 

Z-CX*20.0/(PN*QN) 

BT=3.142*01/OQ-BO 

ZA«BT*BT/(16.0*TS*GC) 

ZB=0.35*BT/TS 

ZC-HO/BO 

ZD«HX*0.333/BS 

ZE-HY/BS 
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IF(ZZ-2.0) 76,77,78 

76 PC»Z*(ZE+7.D+ZA+ZB) 

GO TO 82 

77 POZ*(ZC+(2.0*HT/(BO+BS))+(HW/BS)+ZD+ZA+ZB) 

GO TO 82 

78 IF(ZZ-4.0) 79,80,81 

79 PC“Z*(ZC+(2.0*HT/(B0+B1))+(2.0*HW/(Bl+B2))+(HX*0.333/B2)+ZA+ZB) 
GO TO 82 

80 PC«Z*(ZC+0.62) 

GO TO 82 

81 POZ*(ZE+ZD+(0.5*GC/TS)+(0.25*TS/GC)+0.6) 

82 WO.642*SC*Q0*AC*HM 
ZA-3«1416*(DI+HS)/OQ 
IF(ZZ-3.0) 88,89,88 

88 TM-ZA-BS 
GO TO 86 

89 TM»(3.1416*(DI+2.*HS)/QQ)-B3 

86 WI-(TM*QQ*SS*HS+(DU-HC)*3.142*HC*SS)*.566 
AN-0.0 

100 AN-AN+0.005 
AL-COS(AM) 

IF(PF-AL) 100,100,101 

101 VR-0.262*DR*RPM 
AP«BP*PL*RK 
FH-BG*GE/0.00319 
ZG-PT*FE*0.000001/AS 
FK«RR*ZG 

FR-(T2+234.5)*FK*0.00394 

RC-.321*PT*FE*AS 
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PRINT3,S,CQ,HM,A,RG,XR 

PUNCH1,VA,EE,EP,PN,F,PX 

PUNCH1,RPM,PI,PF,CK,POL,DI 

PUNCH1,DU,CL,SS,HC,PC,ON 

PUNCH 1,WL,BK,ZZ,BO,PE 

PUNCH1,XR,BS,CQ,HX,HY,HZ 

PUNCH1,HS,WC,AC,QQ,W,GE 

PUNCH1,SC,YY,C,TS,BG,TG 

PUNCH1,F0,TT,EL,AG,GA,FGML 

PUNCH1,RG,GC,RP,C1,TP,CP 

PUNCH1,DF,CM,CF,BP,HP,HP1 

PUNCH 1,PL,EB,DR,RK,CC,WR 

PUNCH1,01,WO,HD,DD,H,B 

PUNCH 1,BN,SB,TB,RE,T3,DSH 

PUNCH1,DISH,ALH,TYPY,PC0IL,DC0IL,BC0IL 

PUNCH 1,TY,TYE,TYR,DYC,PT,FE 

PUNCH1,RD,RT,T2,RR,SNL,WF 

PUNCH1,CS,AS,ET,SM,H0,HT 

PUNCH1,B1,B3,SM,DI SHI,PM,P5 

PUNCH1,TM,WI,AN,VR,AP,FH 

PUNCH1,FK,FR,RC,AL,P6,P7 

PAUSE 

END 
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c PASS 4 HOMOPOLAR INDUCTOR 

1 F0RMAT(E11.5,El 1.5,El 1.5,Ell.5,El 1.5,El 1.5) 

3 FORMAT(9X F12.5,2X F12.5) 

READ1, VA,EE,EP,PN,F,PX 

REAOI, RPM,PI,PF,CK,POL,DI 

READ1, DU,CL,SS,HC,PC,QN 

READ 1,WL,BK,ZZ,BO,PE 

READ 1,XR,BS,CO,HX,HY,HZ 

READ1, HS,WC,AC,00,W,GE 

READ1, SC,YY,C,TS,BG,TG 

READ 1,FO,TT,EL,AG,GA,FGML 

READ 1, RG,GC,RP,C1,TP,CP 

READ 1,DF,CM,CF,BP,HP,HP1 

READ 1,PL,EB,DR,RK,CC,WR 

READ 1,D1,WO,HD,DD,H,B 

READ 1,BN,SB,TB,RE,T3,DSH 

READ 1,DISH,ALH,TYPY,PCOIL,DCOIL,BCOIL 

READ 1,TY,TYE,TYR,DYC,PT,FE 

READ 1,RD,RT,T2,RR,SNL,WF 

READ 1,CS,AS,ET,SM,HO,HT 

READ 1,B1,B3,SM,DISH1,PM,P5 

READ 1,TM,WI,AN,VR,AP,FH 

READ 1,FK,FR,RC,AL,P6,P7 

EK=EL/(10.0**(0.103 * Y Y*T S+0.4 02)) 

IF(DI-8.0) 83,83,84 

83 EK-SQRT(EK) 

84 ZF-.612*L0G(10.0*CS) 

EW«6.28*EK*ZF*(TP**(0.62-(0.228*LOG(ZF))))/(CL*DF*DF) 

XL«(PC+EW)*XR*2. 
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XD«=XR*AG*C1*CM 


XQ=XR*CO*AG 
IF(PM)341,340,341 

340 PM«10.02*DR*CL/(PX*(HP1+GC)) 

341 IF(P5)343,342,343 

342 P5=1„675*(DC0IL-PCOIL)*(DCOIL+PCOIL)/BCOIL 

343 IF(P6)345,344,345 

344 P6=2„5*(PC0IL-DI)*(PCOIL+DI)/BCOIL 

345 IF(P7>347,346,347 

346 P7=2.5*(DI+DISH1)*(DU-DI)/(DIJ-DISH1) 

347 RL=(P5+P6+PM*PX/4.)/CL 

XF=>( 1 .0—Cl/(( 1,273*RL/AG+2.0*CP)*CM) )*XD 
S|=PT*PT *(PX*.3 927*CP*AG*PL+RL)*1•E-8 
IF (BN) 3 07,306,3 07 

306 X1*=0 
P2«=0 
X2»0 

GO TO 308 

307 IF(DD)103,103,102 

102 ZG-0.62 
GO TO 104 

103 ZG=0.333*H/B 

104 ZF-HD/W0+ZG+.5 
BD»ZF*6.38 

BE*»(BH-(BP-1 ,0)*TB)*2.127/GE 

P MBD+BE)*RL*COS((BN-1.0)*T B*1.572/TP)/(BD+BE+RL) 
X1=XR*P1 

P2=(ZF+GC/TB)*20.0*TB/TP 

X2-XR*P2 
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308 XA=XL+XD 
XB=XL+XQ 
XU=XL+XF 
XS=0.88*XIJ 

PR I NT3,RP,XL,ET,XD,EB,X0_ 

PR INT3,PC,XA,EW,XB,WC,XF,WI,SI,PM,XI,P5,X2,P6,XU,P7,XS 

PUNCH1,VA,EE,EP,PN,F,PX 

PUNCH1,RPM,PI, PF,CK,P0L,DI 

PUNCH 1,DU,CL,SS,HC,PC,QN 

PUNCH1,WL,BK,ZZ,B0,XD,XQ 

PUNCH 1,XR,BS,PE,HX,HY,HZ 

PUNCH1,HS,AC,Q0,W,GE 

PUNCH1,SC,YY,C,TS,BG,TG 

PUNCHl,FQ,TT,EW,AG,GA,FGML 

PUNCH1,RG,GC,RP,C1,AP,P2 

PUNCHl,DF,WF,CF,FH,BP,PM 

PUNCH 1,P5,P6,P7,HP,HP1 

PUNCH 1,PL,EB,DR,RK,CC,WR 

PUNCHl,D1,W0,SI,DD,H,B 

PUNCH1,BN,SB,TB,RE,T3,DSH 

PUNCHl,D ISH,ALH,TYPY,PC0IL,DC0IL,BC0IL 

PUNCHl,TY,TYE,TYR,DYC,PT,VR 

PUNCHl,RD,RT,SNL,CS,AS,ET 

PUNCHl,FK,FR,XA,XB,AN,AL 

PUNCHl,RC,HD,XS,XL,X1,X2 

PUNCHl,SM,DISHI,TP 

PAUSE 

END 
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C PASS 5 HOMOPOLAR INDUCTOR 

1 FORMAT(E11.5,Ell.5,El 1.5,El 1.5,El 1.5,El 1.5) 
3 FORMAT(9X F12.5.2X FI2.5) 

READ1, VA,EE,EP,PN,F,PX 

READ 1, RPM,PI,PF,CK,POL,DI 

READ1, DU,CL,SS,HC,PC,QN 

READ 1, WL,BK,ZZ,BO,XD,XQ 

READ 1,XR,BS,PE,HX,HY,HZ 

READ 1,HS,AC,QQ,W,GE 

READ1, SC,YY,C,TS,BG,TG 

READ 1,FQ,TT,EW,AG,GA,FGML 

READ1, RG,GC,RP,C1,AP,P2 

READ 1,DF,WF,CF,FH,BP,PM 

READ 1,P5,P6,P7,HP,HP1 

READ 1,PL,EB,DR,RK,CC,WR 

READ 1,D1,W0,SI,DD,H,B 

READ 1,BN,SB,TB,RE,T3,DSH 

READ 1,DISH,ALH,TYPY,PCOIL,DCOIL,BCOIL 

READ 1,TY,TYE,TYR,DYC,PT,VR 

READ 1,RD,RT,SNL,CS,AS,ET 

READ 1 # FK,FR,XA,XB,AN,AL 

READ 1,RC,HD,XS,XL,X1,X2 

READ 1,SM,DI SHI,TP 

IF (BN) 105,105,106 

105 XX-XS 
XY-XB 

GO TO 107 

106 XX-XL+X1 
XY-XL+X2 
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107 XN=(XX+XY)*0.5 
TC=SI/FK 

RA=PN*PI*PI*RP*0.001/VA 
TA-XN/(628.4*F*RA) 

T5=XS*TC/XA 
IF(F-60.0) 108,108,109 

108 T4=0.035 
GO TO 110 

109 T4-0.005 

110 IF(WF)111,111,112 

111 WF=DR**2.5*(RPM**1.5)*PL*0.00000252 

112 IF(ZZ-1.0)301,301,302 

301 BO=BS 

302 GT-BO/GC 

IF(GT-1.0)304,304,303 

304 AA-2.6 

GO TO 115 

303 IF(GT-3.75)113,114,114 

113 AA-1 0.0**0.178/( (GT-1,0)**0.334) 

GO TO 115 

114 AA-10.0**0.11/((GT-1.0)**0.174) 

115 GF-AA*P!*SC/(C*FH) 

305 IF(SC-1.0)121,121,122 

120 A5-0.0 

GO TO 129 

121 AX-1.0 
AY-1.0 

GO TO 125 

122 AX-3.0*YY/(PN*QN)-2.0 
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IF(CS-0.667)123,124,124 

123 AY=1.5*YY/(PN*QN)-0.25 
GO TO 125 

124 AY=.75*YY/(PN*QN)+0.25 

125 A3=AX*P2/AY 

A4=0.07*AX*AG/(CF*CF) 

IF(AX)120,120,126 

126 IF(BN)127,127,128 

127 A5=A4 

GO TO 129 

128 A5=(A4+A3)/(A3*A4) 

129 IF(W) 130,130,131 

130 XO-O.O 

GO TO 132 

13 1 AA**(3.0*HZ+HX)*1,667/(PN*QN*CF*CF*DF*DF*BS) 
X0=((PC+A5)*AX/AY+AA+0.2*EW)*XR 
132 AA=WO/(GC*CC) 

VT=0 

IF(AA)148,147,148 
148 IF(AA-0.65)145,147,146 

145 VT=LOG(10.0*AA)*(-0.242)+0.59 
GO TO 147 

146 VT«0.327-(AA*0.266) 

147 EZ«(ET+EB)*0.5-1.0 
AA-PN*PI*PI 
PU=AA*RG 
PV-AA*RP 

VV-EP*PI*PF*.003 

FSC-XA*FH*0.01*2 
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ACR«=(DU-2.*HC)*3.1416*PE*SS/PX 

ATH=QQ*SS*SM 

PRINT3,AS,XX,FK,XY,FR,XN,RC,XO,WR,TG,VR,FQ 
PRINT348,BG 
348 FORMAT(23X,FI2.5) 

PUNCH1,VA,EE,EP,PN,F,PX 

PUNCH 1,RPM,PI,PF,CK,POL,TB 

PUNCH 1,BO,GC,HP,HP 1,PL,DR 

PUNCH 1,SB,RE,T3,PT,T5,T4 

PUNCH 1,WO,DD,H,BN,GF,VT 

PUNCH 1,SNL,TS,CC,BG,FK,AP 

PUNCH 1,FO,FR,XD,FH,GA,HC 

PUNCH 1,AN,AL,XA,WF,TC,TA 

PUNCH 1,AS,HS,B,GE,BP,XB 

PUNCH1,EZ,PU,VV,FSC,PV,HD 

PUNCH 1,QN,DSH,DISH,ALH,TYPY,PCOIL 

PUNCH1,DC0IL,BCOIL,TY,TYE,TYR,DYC 

PUNCH1,C1,DI,CL,QQ,GA,SM 

PUNCH 1,DU,FGML,ATH,PM, ACR, P5 

PUNCH1,P6,P7,TT,SS,BK,WL 

PUNCH1,TG,DISH1,TP,D1 

PAUSE 

END 
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C PASS 6 HOMOPOLAR INDUCTOR 

DIMENSION AI(90) 

1 F0RMAT(E11.5,Ell.5,El 1.5,El 1.5,El 1.5,El 1.5) 

250 F0RMAT(9X F12.5,2X F12.5/) 

888 FORMAT(FI 0.0,FI 0.0,FI 0.0,FI 0.0,FI 0.0,FI0.0) 

3 FORMAT(9X F12.5,2X F12.5) 

K-1 

823 READ888,AI(K),AI(K+1),AI(K+2),AI(K+3),AI(K+4),AI(K+5) 
K-K+6 

IF(K-89)823,199,199 


199 READ 

1,VA,EE,EP,PN,F,PX 

READ 

1,RPM,PI,PF,CK,P0L,TB 

READ 

1,B0,GC,HP,HP1,PL,DR 

READ 

1,SB,RE,T3,PT,T5,T4 

READ 

1,WO,DD,H,BN,GF,VT 

READ 

1,SNL,TS,CC,BG,FK,AP 

READ 

1,FQ,FR,XD,FH,GA,HC 

READ 

1,AN,AL,XA,WF,TC,TA 

READ 

1,AS,HS,B,GE,BP,XB 

READ 

1,EZ,PU,VV,FSC,PV,HD 

READ 

1,QN,DSH,DISH,ALH,TYPY,PC0IL 

READ 

1,DCOIL,BCOIL,TY,TYE,TYR,DYC 

READ 

1,C1,D1,CL,QQ,GA,SM 

READ 

1,DU,FGML,ATH,PM,ACR,P5 

READ 

1 # P6,P7 # TT,SS,BK,WL 

READ 

1,TG,D1 SHI,TP,D1 


PHM-PM*FH*.001 

FGM»FH+PX*PHM*GE/(.00319*GA) 
BP-(FQ+PHM)/AP 
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X=BP 
NA«=3 1 
K=1 

GO TO 802 

805 FP=HP*AT 
PH5=P5*(FGM+FH)*.00l 
PH6=P6*(FGM+FH)*.001 
PH7=P7*FGM*.OOl 
PSH=(PX*FO/2.)+PX*PHM+PH7 
BT=(TG+PHM*PX)/ATH 

X=BT 

NA*=1 

K=2 

GO TO 802 

806 FT=HS*AT 
BC«(FQ+PHM)/ACR 
X«=BC 

NA=1 

K=3 

GO TO 802 

807 FC=HC*AT 

ASH*= (DSH**2-D I SH**2)*. 7854 

BSH-PSH/ASH 

NA=31 

K=4 

X-BSH 

GO TO 802 

808 FSH«*ALH *AT 
AY-TY*(DU+TY)*3.1416 
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IF(TYPY-2.)816,815,822 
822 ALY=1.334*CL 
GO TO 821 

816 AYR=0 
AYC-0 

ALY“BCOIL+.667*CL 
ALYR=0 
ALYOO 
GO TO 817 
815 ALY=.667*CL 
821 AYODYC+TYE 

AYR=TYR*(DU+2.*TY)*3.1416 

ALYC=*BCOI L 

ALYR-DYC-DU 

817 Z=PSH+PH5+PH6 
X-Z/AY 

BY*=X 
NA=61 
K»5 

GO TO 802 

809 FY“ALY*AT 
IF(TYPY-1.)818,819,818 

818 X-Z/AYC 
BYC=X 
NA-61 
K-6 

GO TO 802 

810 FYOALYC*AT 
X-Z/AYR 
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NA=*61 
K=7 

GO TO 802 
811 FYR»ALYR*AT 

820 FNL=2.*(FGM+FT+FC+FP)+FSH+FY+FYC+FYR 
PRINT3,TC,BT,TA,BC,T5,FT,T4,FC 
PRINT250,FSC,FH 
PUNCH1,ATH,FT,FGML, PM, FNL,FH 
PUNCH1,GE,PX,GA,FO,XA,PHM 
PUNCH1,AP,HP,D1,P5,P6,P7 
PUNCH 1,ALH,BT,HS,ACR,HC,AY 
PUNCH1,ALY,AYC,ALYC,AYR,ALYR 
PUNCH1,PT,FK,FR,AS,VA,EE 
PUNCH1,PN,F,RPM,PI,PF,P0L 
PUNCH!,BO,GC,DR,SB,RE,T3 
PUNCH1,WO,DD,H,BN,GF,VT 
PUNCH1,SNL,TS,CC,BG,WF,TB 
PUNCH1,B,BP,XB,EZ,PU,VV 
PUNCH!,FSC,PV,HD,ON,XD,EP 
PUNCH 1,SM, DU, SS,O.Q,BK,WL 
PUNCH1,TT,BC,BSH,ASH,TG,BYC 
PUNCHl,TP,AN,AL,TYPY,C1,FGM 
PUNCH1,CK 
PAUSE 

802 IF(AI(NA)-X)830,831,831 
831 NA-NA+3 

835 1F(AI(NA)-X)833,834,834 
833 NA-NA+2 
GO TO 835 

PA-43 



834 AA-AI (NA) 


BB1“AI(NA-2) 

DC-AI(NA+1) 

D-AI(NA-I) 

XX-(AA-BB1)/(.4343*(L0G(DC)-L0G(D+.0001))) 
Y«AA-XX*.4343*LOG(DC) 
AT-EXP(2.306*(X-Y)/XX) 

GO TO (805,806,807,808,809,810,811),K 
830 PR I NTS50 

850 FORMAT(17HMACHINE SATURATED) 

PAUSE 
819 FYC-0 
FYR-0 
BYC-BY 
GO TO 820 
END 
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C PASS 6a HOMOPOLAR INDUCTOR 

1 FORMAT(E11.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 
READ 1,ATH,FT,FGML,PM,FNL,FH 
READ 1,GE f PX,GA,FQ,XA,PHM 
READ 1 t AP,HP,Dl,P5,P6,P7 
READ 1,ALH,BT,HS,ACR,HC,AY 
READ 1,ALY,AYC,ALYC,AYR,ALYR 
READ 1 ^PTjFKjFRjASjVAjEE 
READ 1,PN,F,RPM,PI,PF,POL 
READ 1,BO,GC,DR,SB,RE,T3 
READ 1,WO,DD,H,BN,GF,VT 
READ 1,SNL,TS,CC # BG,WF,TB 
READ 1,B,BP,XB,EZ,PU,VV 
READ 1,FSC,PV,HD,QN,XD,EP 
READ l,SM t DU f SS,QO,BK f WL 
READ 1,TT,BC,BSH f ASH,TG,BYC 
READ 1,TP,AN # AL,TYPY,C1,FGM 
READ 1,CK 
IF (RE >309,309,311 

309 WD-0.0 
WU-0.0 
GO TO 178 

311 FS1-2.0*QN*PN*F 
FS2»2.0*FS1 
M **0 

150 IF(H-1)151,152,178 

151 RM-RE 

GO TO 153 

152 RM»RE*(T3+234.5)/254.5 
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153 AA=(FS1/RM)**0.5*DD*0.32 
AB-(FS2/RM)**0.5*DD*0.32 
I F(AA-2.5) 160,160,161 

160 VI-1.0-0.15*AA+0.3*AA*AA 
GO TO 162 

161 V1-AA 

162 IF(AB—2.5) 163,163,164 

163 V2-1.0-0.15*AB+0.3*AB*AB 
GO TO 165 

164 V2-AB 

165 IF(H-B) 167,166,167 

166 VOO.75/VI 
GO TO 169 

167 IF(DD) 166,168,166 

168 VC«H/(3.0*B*V1) 

169 VS-HD/WO+VT+VC 
VG«TB/(CC*GC) 

Q1“1•0—(1•0/(((B0*0.5/GC)**2.0+1.0)**0.5)) 
QZ-BO/TS 

Q2»1.05*SIN(0Z*2.844) 

1F(QZ-0.37)170,170,171 

170 Q3-0.46 
GO TO 172 

171 Q3-0.23*SIN(10.46+QZ-2.1)+0.23 

172 Q4-S1N(6.283*TB/TS-1.571)+1.0 
Q5-SI N( 12.566*TB/TS-1.570+1.0 
IF(H)173,173,174 

173 AB-0.785*DD*DD 
GO TO 175 
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174 AB«H*DD 

175 W2«PX*BN*SB*RM*1.246/(AB*1000.) 
VO-(Q2/(2.0*VS+(VG/0_4)))**2.0*V1 
W5-(Q3/(2.0*VS+(VG/Q5)))**2.0*V2 
WD-(TS*BG*Q1*CC)**2.0*W2*(W3+W5) 

IF(M-1)176,176,177 

176 WU-WD 

177 GO TO 150 

178 PUNCH1,ATH,FT,FGML,PM,F NL,FH 
PUNCH1,GE,PX,GA,FQ,XA,PHM 
PUNCH1,AP,HP,P5,P6,P7 
PUNCH1,ALH,BT,HS,ACR,HC,AY 
PUNCH1, AL Y,AYC,ALYC,AYR,AL YR 
PUNCH1,PT,FK,FR,AS,VA,EE 
PUNCH1,PN,F,RPM,PI,PF,POL 
PUNCH1,B0,GC,DR,SB,RE,T3 
PUNCH1,WO,DD,H,BN,GF,VT 
PUNCH1,SNL,TS,CC,BG,WF,TB 
PUNCH1,B,BP,XB,EZ,PU,VV 
PUNCH1,FSC,PV,HD,QN,XD,EP 
PUNCHl,SM,DU,SS,QQ,BK,WL 
PUNCH1,TT,BC,BSH,ASH,TG,BYC 
PUNCH1,TP,AN,AL,TYPY,WU,WD 
PUNCH1,C1,FGM,D1,CK 

PAUSE 
END ‘ 
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D ! HE MS I ON G X (h ) , YA (4 } , ED ( 4 ) , EGX (4) 

1 FORMAT(E11.5,E11 . 5 , E11 . 5,E11 . 5,E11 . 5,E11. 


,-i ft 

r<£ ml* 

1,ATH 

i > FT,rGMI.,PM 

,FNL,FH 

ft lAD 

1,6E, 

PX,GA,FQ,XA 

, PHM 

READ 

1, AP , 

HP,P5,Pt>,P7 


R L AD 

1 , ALH 

,BT,HS,ACR, 

HC, AY 

READ 

1,ALY 

, AYC , AL. YC , AYR , AL YR 

READ 

1#PT. 

FK,FR,AS,VA 

• EE 

READ 

1,PN, 

F,RPM,PI,PF 

, POL 

READ 

1,B0, 

GC,DR,SB,RE 

#T3 

READ 

1 ,wo. 

DD # HjBNjGF f 

VT 

READ 

1 , SNL 

,TS,CC,BG,WF,TB 

READ 

i,B,BP,X8,EZ,PU, 

vv 

READ 

1,FSC 

,PV,HD,QN,Xi 

D,EP 

READ 

1 , SM, 

DU,SS,QQ,BK 


READ 

1,TT, 

BC,BSH,ASH, 

TG,BYC 

READ 

1,TP, 

AN,AL,TYPY,WU,WD 

READ 

1.C1, 

FGM,D1,CK 



D 2=BG**2.5*0.000061 

D3»(0.0167*QQ*RPM)**1.65*0.000015147 
IF(TS-O.9) 133,133,134 

133 D4-TS**1.285*0.81 
GO TO 137 

134 IF(TS—2.0) 135,135,136 

135 D4-TS**1.145*0.79 
GO TO 137 

136 D4-TS**0.79*0.92 

137 D7-B0/GC 
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IF(D7-1.7) 138,138,139 

138 D5-D7**2.31*0.3 
GO TO 144 

139 IF(D7—3.0) 140,140,141 

140 D5=D7**2.0*0.35 
GO TO 144 

141 IF(D7-5.0) 142,142,143 

142 D5-07**1.4*0.625 
GO TO 144 

143 D5-D7**0.965*1.38 

144 D6=10.0**(0.932*C1-1.606) 

WN-D1*D2*D3*D4*D5*D6*GA 

WT=(S H)*Q0*SS*HS*0.453*(BT/BK)**2.0*WL 
WQ»(DU-HC)*1.42*HC*SS*(BC/BK)**2.0*WL 
YA(1)-100. 

YA(2)=66.66667 
YA(3)»50. 

!F(POL)310,320,310 
320 ED(4)*=0 
GX(4)-0 
FGX(4)=0 
JA-3 

GO TO 330 
310 YA(4)=100./POL 
JA-4 

330 DO 99 K-1,JA 

AA-ATAN((XB/YA(K)+SIN(AN))/AL) 

BB-AA-AN 

ED(K)-XA*SIN(AA)/YA(K)+COS(BB) 





r OX ( K.}*-•:M|. v \ f;-i, /YA(K) 

iF(PF~„95)21.3,213,212 

212 GX(K)«rO*CK 
GO TO 99 

2 13 GX(K)»(rD(K)-{. 33 *XD*SIN(AA}/YA(K)))*F0 
99 CONTINUE 


|FfPOL)820,021,820 

820 AJ=4 

GO 10822 

821 AJ=3 

822 PUNCH1,ED(1),ED(2),ED(3),ED(4),AJ,BC 
P U N C H1 , G X (1) , G X (2) , G X ( 3 ) , G X (4 ) , F G M 
PUNCH 1,FGX(1),FGX(2),FGX(3),FGX ( 4 ) 
PUNCH1,FT,PF,PM,FH,PX 
PUNCH1,GE,AP,HP,P5,P6 

PUNCH 1,P7,PHM, ASH ,AL.H,BT,ATH 
PUNCH 1,HS, HC ,ACR,TYPY,AY,AYC 
PUHC H1 , AYR , At. Y , At.YC , AL YR , B P , B S H 
PUNCH 1,BYC,F NL,TG,G A, FQ , TP 
PUNCH 1, POL,PU,F K, PV f FR ,WU 
PUNCH1,WD,PT,AS,WF,WQ,GF 
PUNCH1 ,WT,XA,VV,WN,EZ,SNL 
PUNCH1,FH,PX,GE,AP,HP 
PUNCH1,P5 # P6,P7,ASH,ALH,ATH 
PUNCH1,HS,HC,ACR,TYPY,AY,AYC 
PUNCH 1, AYR , AL Y, AL YC , AL YR, BC, BT 
PUNCH1,TG,GA,FQ,TP,EP,PM 
PAUSE 
END 
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c PASS 8 HOMOPOLAR INDUCTOR 

4 FORMAT (F11.3,8X Fl1.3,F11.3,F11.3,F11.3) 

888 FORMAT (FI O.O, FIO.O,FIO.O,FIO.O,FIO.O,FtO.O) 

DIMENSION FFL(5),BSHLL(4),BCL(4),BTL(4),BPLL(4),BYCL(4),FGLL(4) 
DIMENSION GX(4),ED(4),PMLL(4),AI(90),FGX(4) 

1 FORMAT (Ell. 5, Ell. 5, Ell. 5, Ell. 5, Ell. 5, Ell. 5) 

K=1 

823 READ888,Al(K),AI(K+1),AI(K+2),AI(K+3),Al(K+4),AI(K+5) 

K=K+6 

IF(K-89)823,199,199 

199 READ 1,ED(1),ED(2),ED(3),ED(4),AJ,BC 
READ 1,GX(1),GX(2),GX(3),GX(4),FGM 
READ1 ,FGX(1),FGX(2),FGX(3),FGX(4) 

READ 1,FT,PF,PM,FH,PX 
READ1 ,GE,AP,HP,P5,P6 
READ 1,P7,PHM,ASH,ALH,BT,ATH 
READ 1,HS,HC,ACR,TYPY,AY,AYC 
READI,AYR, ALY,ALYC,ALYR,BP,BSH 
READ1,BYC,FNL,TG,GA,FO,TP 
DO 886 J=1,4 
PMLL(J)»0 
FFL( J+1 )*»0 
BSHLL( J)**0 
BCL(J)“0 
BTL(J)“0 
BPLL(J)“0 
BYCL(J)»0 
886 FGLL(J)*»0 
JA»AJ 
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DO 885 
PPL=GX(J) 

EDD»ED(J) 

FGML*=FGX( J) 

FTL*»FT*( 1 .+PF) 

PML»PM*(FGML+FH*EDD)*.OOl 

ZZ-PX*GE/(.00319*GA) 

FGL-FH*EDD+PML*ZZ 

BPL“(PPL+PML)/AP 

NA-31 

K=1 

X-BPL 

GO TO 802 

870 FPL»AT*HP 
Z-FGL+FTL+FPL 
PH5L-P5*.002*Z 
PH6L“P6*.002*Z 
PH7L»P7*.001*Z 

PSHL=(PPL*PX/2.)+PX*PML+PH7L 

BSHL-PSHL/ASH 

X-BSHL 

NA-31 

K-2 

GO TO 802 

871 FSHL-ALH*AT 

pml-pm*(fgml+fgl)*.ooi 

PHLL(J)-PML 

PPLL-PPL+PML 

X-PPLL/AP 
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BPLL(J)=X 
NA**3 1 
K=3 

GO TO 802 

872 FPL«AT*HP 
FGL“FH*EDD+PML*ZZ 
FGLL(J)=FGL 
X=(TG+PX*PML)/ATH 
BTL( J)*»X 

NA=1 

K»4 

GO TO 802 

873 FTL-AT*HS 
Z-FTL+FGL+FPL 
PH7L»P7*Z*.OOl 

PSHL«(PPLL*PX/2.)+PX*PML/2.+PH7L 

X=PSHL/ASH 

NA=3 1 

K=5 

BSHLL(J)=X 
GO TO 802 

874 FSHL=ALH*AT 
Z“2.*Z+FSHL 
PH5L=P5*Z*.OOl 
PH6L»P6*Z*.001 
PCL=PPLL+(PH5L+PH6L)/PX 
X=PCL/ACR 

BCL( J)**X 
NA-1 
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K=6 

GO TO 002 


875 FCL=AT*HC 

IF(TYPY-1.)880,881,880 

881 PY-PSHL+PH6L+PH5L 
GO TO 882 

880 PY«*PSHL+PH6L 

882 X=PY/AY 
NA=61 
K«=7 

GO TO 802 

876 FYL=AT*ALY 

IF(TYPY-1.)883,884,883 

883 PY=PSHL+PH6L+PH5L 
X=PY/AYC 
BYCL(J)=X 

NA=61 
K=8 

GO TO 802 

877 FYCL=AT*ALYC 
X=PY/AYR 
NA=61 

K-9 

GO TO 802 

878 FYRL«AT*ALYR 

889 FFL(J+1)-2,*(FGL+FTL+FCL+FPL)+FSHL+FYL+FYCL+FYRL 
885 CONTINUE 

PR INT4,PHM,PMLL(1),PMLL(2),PMLL(3),PMLL(4) 

PR INT4,FGM,FGLL(1),FGLL(2),FGLL(3),FGLL(4) 
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PR INT4,BP,BPLL(1),BPLL(2),BPLL(3),BPLL(4) 

PRINT4,BT,BTL(1),BTL(2),BTL(3),BTL(4) 

PUNCHi t BT,BTL(1),BTL(2),BTL(3),BTL(4) 

PRINT4,BSH,BSHLL(1),BSHLL(2),BSHLL(3),BSHLL(4) 

PRINT4,BC,BCL(1),BCL(2),BCL(3),BCL(4) 

PUNCH1,BC,BCL(1),BCL(2),BCL(3),BCL(4) 
PRINT4,BYC,BYCL(1),BYCL(2),BYCL(3),BYCL(4) 
PUNCH1,FNL,FFL(2),FFL(3),FFL(4),FFL(5) 

DO 887 K»1,8 
READ 1,R1,R2,R3,R4,R5,R6 
887 PUNCH1,R1,R2,R3,R4,R5,R6 
PAUSE 

884 FYCL-0 
FYRL-0 
BYCL(J)-X 
GO TO 889 

802 IF(AI(NA)-X)830,831,831 
831 NA-NA+3 

835 IF(AI(NA)-X)833,834,834 

833 NA-NA+2 
« 

GO TO 835 

834 AA-AI(NA) 

BB1-AI(NA—2) 

DC-A1(NA+1) 

D-AI(NA-1) 

XX«(AA-BB1)/(.4343*(LOG(DC)-LOG(D+.0001))) 
Y-AA-XX*.4343*L0G(DC) 

AT-EXP(2.306*(X-Y)/XX) 

GO TO (870,871,872,873,874,875,876,877,878),K 
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830 GO TO 885 
PAUSE 
END 
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C PASS 9 HOMOPOLAR INDUCTOR 

DIMENSION PR(5),FI(5),PS(5),G(5),DL(5),PP(5),EX(5),ST(5),VA(5) 
DIMENSION P(5),E(5),PZ(5),SP(5),FFL(5),CDD(5),EF(5) 

DIMENSION BTL(5),BCL(5),WQL(5) 

1 FORMAT(E11.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 

4 FORMAT (F11.3.8X F11.3,F11.3,F11.3,F11.3) 
READ1,BT,BTL(2),BTL(3),BTL(4),BTL(5) 
READ1,BC,BCL(2),BCL(3),BCL(4),BCL(5) 
READ1,FNL,FFL(2),FFL(3),FFL(4),FFL(5) 

READ 1,POL,PU,FK,PV,FR,WU 

READ 1,WD,PT,AS,WF,WQ,GF 

READ 1,WT,XA,VV,WN,EZ,SNL 

READ 1,FH,PX,GE,AP,HP 

READ 1,P5,P6,P7,ASH,ALH,ATH 

READ 1,HS,HC,ACR,TYPY,AY,AYC 

READ 1,AYR,ALY,ALYC,ALYR,BC,BT 

READ 1,TG,GA,FQ,TP,EP,PM 

PRINT4,FNL,FFL(2),FFL(3),FFL(4),FFL(5) 

FFL(1)=FNL 
BTL(1)=BT 
BCL(1)«BC 
G( 1 )*»0 
G(2)=1 
G(3)**1.5 
G(4)»2. 

G(5)«POL 

PW-PU 

FW*»FK 

WW»WU 
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DO 183 M=1,5 
FI (M)*=FFL(M)/PT 
EF(H)=FI(M)*FW 
CDD(M)=FI(H)/AS 
UA«G(M) 

PR ( M) **F I ( M) *F I (M) *FW 
IF (FI (M)) 198,197,198 
198 PS(M)=PW*UA*UA 

WQL (M) =*WQ* (BCL(M) /BC) **2 
X=*WF+WOL(M) 

GM = (GF*IJA)**2.0+1.0 

ST(h)=(2.0*(0.0027*XA*UA)**1.8+1.0)*WT*(BTL(M)/BT)**2 
VA(M)*»VV*UA 
181 DL(M)*»GM *WW 
PP(M)»GM *WN 
EX(M)=EZ*PS(M) 

SP(M)**PP(M)+DL(M)+PR(M)+PS(M)+EX(M)+ST (H)+X 
P(MMSP(M)/1000.)+VA(M) 

IF(GM)185,184,185 

184 PZ(M)*=0 
E(M)~0 

GO TO 186 

185 PZ(M)=(SP(M)/P(M))*.1 
E (M) *» 100.0-PZ (M) 

186 FW=*FR 
WW-WD 

183 PW-PV 

PRINT4,FI(1),FI(2),FI(3),FI(4),FI(5) 

PRINT4,EF(1) # EF(2),EF(3),EF(4),EF(5) 
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PRINTS,CDD(1),CDD(2),CDD(3),CDD(4),CDD(5) 
PRINT4, PR(1),PR(2),PR(3),PR(4),PR(5) 
PRINT4, WF,WF,V/F,WF,WF 
PRINT4, ST(1),ST(2),ST(3),ST(4),ST(5) 
PRINT4,WQL(1),WQL(2),WQL(3),W0_L(4),WQL(5) 
PRINT4, PP(1),PP(2),PP(3),PP(4),PP(5) 
PRINT4, DL(1),DL(2),DL(3),DL(4),DL(5) 
PRINT4, PS(1),PS(2),PS(3),PS(4),PS(5) 
PRINT4, EX(1),EX(2),EX(3),EX(4),EX(5) 
PRINT4, SP(1),SP(2),SP(3),SP(4),SP(5) 
PRINT4, VA(1),VA(2),VA(3),VA(4),VA(5) 
PRINT4, E(1),E(2),E(3),E(4),E(5) 

1F(SNL)191,1 91,190 

190 PUNCH 1,FH,PX,GE,AP,HP 
PUNCH1,P5,P6,P7,ASH,ALH,ATH 
PUNCH1,HS,HC,ACR,TYPY,AY,AYC 
PUNCH1,AYR,ALY,ALYC,ALYR,BC,BT 
PUNCH1,TG,GA,F0,TP,EP,PM 

191 PAUSE 

197 PS(M)-0 

G^O 
ST (M)**0 
WQL(M)«0 
X-0 

VA(H)-0 
GO TO 181 
END 
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c 


PASS 10 HOMOPOLAR INDUCTOR 
DIMENSION AI(90) 


245 F0RMAT(F12.5,F12.5,F12.5,F12.5,F12.5,F12.5) 

246 F0RMAT(F12.5,F12.5,F12.5,F12.5,F12.5,F12.5//) 

1 FORMAT(El1.5,E11.5,E11.5,E11.5,E11.5,E11.5) 

888 FORMAT(F10.0,FI 0.0,FI 0.0,FI 0.0,FI 0.0,FI0.0) 

K-1 

823 READ888,Al (K),AI (K+1),AI(K+2),AI(K+3),AI(K+4),AI(K+5) 
K-K+6 

IF(K-89)823 , 199,199 
199 READ 1,FH,PX,GE,AP,HP 

READ 1,P5,P6,P7,ASH,ALH,ATH 

READ 1,HS,HC,ACR,TYPY,AY,AYC 

READ 1,AYR,ALY,ALYC,ALYR,BC,BT 

READ 1,TG,GA,FQ,TP,EP,PM 

R-.7 

FF-TG 

FG-FH 

DO 247 L-1,9 

R-R+.1 

TG-FF*R 

FH-FG*R 

V-1.732*EP*R 

PHM-PM*FH*.001 

FGM-FH+PX*PHM*GE/(.00319*GA) 

BP-(FQ+PHM)/AP 

X-BP 

NA-31 

K-1 


PA-60 



GO TO 802 

805 FP»HP*AT 
PH5=P5*(FGM+FH)*.OOl 
PH6-P6*(FGM+FH)*.001 
PH7=P7*FGM*.001 

PSH« (P X*F0_/2 .) +PX*PHM+PH7 

BT=(TG+PHM*PX)/ATH 

X-BT 

NA-1 

K-2 

GO TO 802 

806 FT=HS*AT 

BC*= (FQ+PHM) / ACR 

X»BC 

NA=1 

K=3 

GO TO 802 

807 FOHC*AT 
BSH»PSH/ASH 
NA«3 1 

K=*4 
X=BSH 
GO TO 802 

808 FSH«ALH *AT 
817 Z-PSH+PH5+PH6 

X-Z/AY 

BY-X 

NA-61 

K-5 
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GO TO 802 

809 FY-ALY*AT 
IF(TYPY-1.>818,819,818 

818 X-Z/AYC 
BYOX 
NA-61 
K=6 

GO TO 802 

810 FYOALYC*AT 
X-Z/AYR 
NA=61 

K»7 

GO TO 802 

811 FYR»ALYR*AT 

820 FNL-2.*(FGM+FT+FC+FP)+FSH+FY+FYC+FYR 
PRI NT245,V,FGM,BC,FC,BT,FT 
247 PRI NT246,BP,FP,BSH,FSH,BYC,FNL 
PAUSE 

802 IF(AI(NA)-X)830,831,831 
831 NA-NA+3 

835 IF(AI(NA)-X)833,834,834 

833 NA*»NA+2 
GO TO 835 

834 AA-AI(NA) 

BB1-AI(NA-2) 

DC-AI(NA+1) 

D-AI(NA—1) 

XX-(AA~BB1)/(.4343*(L0G(DC)-L0G(D+.000l))) 
Y-AA-XX*.43 43*L0G(DC) 
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AT«EXP(2.306*(X-Y)/XX) 

GO TO (805,806,807,808,809,810,811),K 
830 PRINT850 

850 FORMAT(17HMACHINE SATURATED) 

PAUSE 
819 FYOO 
FYR-0 
BYC-BY 
GO TO 820 
END 
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SECTION RA 


PERMANENT-MAGNET A-C GENERATOR 
COMPUTER PROGRAM AND TEST DATA 



SECTION RA 




GAP STATOR WINDING STATOR SLOT STATOR STACK 



MTR'L DAMPER BAR ROTOR STACK 
























































































































































STATOR 


1 ■. > K 


(OUTPUT) 


MODEL 


... i *0 




[CARTER COEFFICIENT 


'(42) (K a k) 


(43) <K d ) 


(44) (K. ) 


SKEW FACTOR 


DIST. FACTOR 


PITCH FACTOR 






■ # 5 91 




-_JL Sl^blkL 




.41900 


±isim 


AIK GAP PERM 


EFFECTIVE AIR GAP 


FUHO^MAX OF FLO. FLUX 


WIMPING CONST. 


POLE CONST. 


(j6?) <£,) j 

wr ( t ] 

i ? 0 a) ( 0 )~U 


(a * ) 


(71) (C ) ) 


(72) (C w I 


(73) (C ? ) 


1(43) ( J 


EFF.CONDUCTORS 


625.94000 


3.66032 


END. EXT. ONE TURN 


(46) (L F ) 


K4«) (o c ) 


COND. AREA 


[(47) (S . ) 


CURRENT DENSITY (STA.) 


ilQ 2 Q 2 L 


2MJL. 1QQQQ 


.84000 


DEMAGNETIZING FACTOR 


(74) (C M ) 


id.QflQ, 


CROSS MAGNETIZING FACTOR 


2*LJ£aJJ 


K<») ( ,) 


1/3 MEAN TURN LENGTH 


5.28530 


265.62000 


AMP COND/IN 


028) (A) 


(33) <R p h ) 


COLD STA. RES. • 20° C 


l3iiq. 


.^S2k5. 


REACTANCE FACTOR 


(129) (X) 


H54) (R p h ) 


HOT STA. RES. • X C 


Jk lLSL 


12.89 Qn 


LEAKAGE REACTANCE 


(130) (X q ) 


K33) (EFtop) 


EDDY FACTOR TOP 


1 T nnnnn 




SYN REACT DIRECT AXIS 


(133) (Xd) 


1(56) (EFbot) 


EDDY FACTOR BOT 


(62) ( |) 


STATOR COND. PE RM. 


l.ononn 


. DQ .QQD , 


9.. S 3 7 2Q 


J1QQQQ 


DAMPER 

LEAKAGE REACT 


(163) (Xpd) 


(165) (XDq) 


1(64) ( c ) 


END PERM. 


,Z..7 a il.Q. 


...12^.8.9Q.QQ. 


UNSAT. TRANS. REACT 


(166) (X *du) 


K65) ( ) 


WT. OF STA COPPER 


mn —r 


2.227 IQ 


11-34300 


SUB. TRANS. REACTDIRECT AX. 


(166) (X M a ) 


WT. OF STA IRON 


3 519 010 


11.31*300 


SUB. TRANS. REACT QUAD AX. 


069) (X” q ) 


K41) 


( p ) 


POLE PITCH 


■■U9633Q.. 


UOlLll 


NEG. SEQUENCE REACT 


(170) (X 2 ) 


WOf) P I 


PERMEANCE IN STATOR 


46,79500. 




ZERO SEQUENCE REACT 


(172) (Xq ) 
(«•) ( 0 | > 


K510) P, 


PERMEANCE OUT STATOR 


23.60000 


328.8&QQQ 


TOTAL FLUX 


KS07) P r 


PERMEANCE MAGNET 


90 05 5 


.01000 


FLUX PER POLE 


(W) ( 0 „) 


I»1D fc 


PERMEANCE AIR GAP 


162.55QQQ 


17.175QQ 


GAP DENSITY 


(93) (Bp) 


»1S7) ( _ )_ 


WT. OF ROTOR IRON 
WT. OF MAGNETS 


■Q0Q00 


4Q^LQ0. 


TOOTH DENSITY 


(91) <B t ) 


-2.-24Q5Q 


31.2.82 QQ 


CORE DENSITY 




V*5) ,V r ) 


PERIPHERAL SPEED 


7788.70000 


18.47.701 


SHORT CIRCUIT AMPS 


l 5 ? 2 ! 1 L d 


vo 

OJ 

%o 

S 


U 

P 


MO 



< 

X 

u 


*/V> (102o) POLE FLUX 

36,01000 


(Bp) (103.) POLE DENSITY 

>6 69200 

(F&W) (163) FtW LOSS 

77.99700 

77-99700 

(F 8 .W) (183) 

(Wmi) (184) STA TOOTH LOSS 

... 7.73540 

8.02617 

(W,„| ) (242) 

(W e ) (185) STA CORE LOSS 

10-32300 

10.32300 

(W e ) 085) 

(Wpnl) 086) POLE FACE LOSS 

59.56300 

98.92551 

(Wpf| ) (243) 

(Wdnl) 093) DAMPER LOSS 

.00000 

.00000 

(Wdll ) (244) 

(l*R,)(194) STATOR CU LOSS 

-00000 

35.76800 

(12 R.) (245) 

( -) (195) EDDY LOSS 

-00000 

-00000 

( - ) (246) 

( - ) (196) TOTAL LOSSES 

15 *;, 61840 

231.03968 

( - ) (247) 

( - ) ( - ) RATING (KW) 

,00000 

- 99590 

( - ) C48) 

( - ) ( - ) RATING Si LOSSES 

t 15561 

1-22693 

( - ) (249) 

( - ) ( - ) PERCENT LOSSES 

_ 100-00000 

18-83056 

( - ) (250) 

(-> (-) PERCENT EFF. 

_nnona.. 

_ 81.16944 

<-) (251) 


VOLTS • 0 LOAD AMPS 

99,78900 

VOLTS # 1/4 LOAD AMPS 

98,92 900 

VOLTS • 1/2 LOAD AMPS 

07 £0700 

VOLTS • 3/4 LOAD AMPS 

: 81200 j 

VOLTS • 4/4 LOAD AMPS 

93-52800 

VOLTS . 5/4 LOAD AMPS 

90.72900 

VOLTS • 3/2 LOAD AMPS 

_ 87-37800 


* 

MO 

ui 

a 

o 


REV. 3 


T) A ~Z 


OP » f'T Akirc 5 outo ft ajt c 



























































































INPUT PARAMETERS FOR PERMANENT MAGNET DESIGN 


1.05 

115. 

66.5 

3. 

400. 

5.54 

W 

1.625 

0 . 

- 0 . 

.92 

.19 

.266 

.226 

.031 

0 . 

36. 

1 . 

0 . 

18. 

5. 

.25 

.25 

0 . 

0 . 

60. 

1.075 

.419 

J>£7- 


.84 

1.625 

1.625 

.71 

3.716 

1 . 

0 . 

0 . 

0 . 

0 . 

0 . 



PO/PM 

CURVE 


340. 

0 . 

1500. 

20. 

1100. 

35. 

1000. 

85. 

500. 

150. 

385. 

187 

280. 

340. 

200. 



6. 

8000. 

5. 

1.0 

3.75 

15. 

77.4 

3. 

.090 

.186 

0 . 

0 . 

.51 

.05 

0 . 

1 . 

.0508 

1 . 

1 . 

0 . 

0 . 

V/l 

• 

.694 

.017 

.02 

M 

1.38- 

1 . 

.T8Z 

.8T2 

0 . 

7. 

0 . 

0 . 

0 . 

-o. 

0 . 

0 . 

.9 

10. 


1250. 

32. 

640. 

130. 

315. 

257. 


RA-4 



JALCTJl 

5 ERMA* 


SHORT- 


fHE TI 


_i ; : 

*TED AND TEjST PEI 'FORMA VCE FC R A | 

SWT MAGNET A~C GENERATOR BEFORE | 

3IRCUATING THE jTEST POINT jAT SH^RT 
C-MAS-I-OBTAX (JED A^ 6 QQQi ^REM AND- DliiCE 

IERAT0R DEMjAGNET^ZED SLIGHTLY UPON 
’IRCU^T, THE CALlJULATEt) POINT AND' 

5T POINT cannot He expected !to MA^CH. 


4-4*0 






T“ 80 " 


—f-*o- 


i 

i 


















Qj_ 'TO.A J l a ') b! 0 - A/ 3d a 


TESTS SHOWING THE EFFECT 
OF SHORT-CIRCUITS ON THE 
OUTPUT VOLTAGE OF THE 
EXAMPLE PERMANENT MAGNET 
A-C GENERATOR 



O / 3 ^ 6 7 3 


6PEEP IN THOU$^nD 


RA-6 





















TERMINAL VOLTS LINE-to-LINE 


TEST 1 ERFORMANCE OF A P.M.G. AT 8000 RPM 
GENERATOR WAG SUPPLYING A RESISTIVE LOAD 
AND HAD BEEN SHORT-CIRCUIT STABILIZED AT 
6000 RPM 



OUTPUT AMPERES PER PHASE 


RA-7 









200 





COMPUTER PROCEDURE FOR 


PERMANENT MAGNET DESIGN CALCULATIONS 


1. 

2 . 

3 . 


c 


a 


9* 


10 . 


Clear core (no switch control). 

Insert output Form #1 into typewriter, set margin for correct 
output, and set typewriter for single spacing* 

Load pass #1 followed by input parameters (output punched cards). 
Load pass #2 followed by output from pass ftl (output printed plus 
punched cards)» 

Load pass if 3 followed by output from pass 42 (output pnnteo. plus 
punched cards). 


Load pass followed. 

by 

output 

from pass #3 

Load pass #5 followed 


output 

from pass 

Load pass #6 followed 

by 



and output from pass i 

f 5 

(output 

. punched car 

Lo a d pas:; f7 fol.1 o ’ ’ c i 

[ •,/ 

Mat pat 

i rom 

r-unched cares). 





output punched :crds). 
output p onehed cards)» 
FO/rIN values^ 


■ i r,. 


h f'ju 


Load pa s? 


: oxi 


we a 


by 


outT- at 


"om 


at crintec 


o 


HA-10 



ALL INPUT PARAMETERS ARE IN FORMAT F7.0 (FIG. I) 


• I • III sit 


aaaaada 


999999 $|9 99999 19 99999 19 99999 jig 

' * * 4 s i 7 |i *** 


a a a a a « a a a a i ala 


a a a a a a a 


100 • #001 # 1 *01 10 # 2 # 0 1000 I 10 

:;^^!HTH7;;T;;TT^HSTT;5757T5!;;HT;5H^TH^5^^f^:?^ !!!; ’ , * , . 

.... . ........ 

5555S5555555555555S555S55555555555 5|5 55555 5p 
566SS(SE6S(SSGSS6SS6S(6SE6S66 
M 7 7 7 7 7|7 7 7 7 7 7 , ,., 

7 7 7 7 7 7|7 7 7 7 7 7 7|7 7 7 7 7 7 7|7 7 7 7 7 7 7|7 7 7 7 7 7 7b 


I I I I I III 


5 5 5 5 5 5|5 5 5 5 5 5 5p 5 5 5 5 5 5j5 5 5 5 5 5 5|5 5 5 5 5 5 5 5 5 5 

6 6 G S 6 6|S 6 6 C C G 'fijs S 6 6 6 S S|S C 6 6 6 6 616 6 6 6 6 6 G|6 6 G 6 6 G s|s 6S566GS66 

77777 7|7 77777 ^7 77777 7|7 77777 7|7; 


I I I 0 I lit 




777777777 

tint ib 11 111 i|i mu ijiiiiiiiiii 

8 9 9 9 9 S|9 ! 


•5 « 17 U (I 70*71 


7? 73 N ft* 77 76 79 W 


FIG. 1 


ALL SATURATION CURVE VALUES ARE IN FORMAT F10.0 (FIG 2 ) 
(ALL SATURATION CURVES MUST HAVE 5 CARDS) 


100 . 

908010000 o|a 

hhiimi 

<1 I 1 11 I 1 1 ill 


10 . 

ooooeooo o|o 

a o h ft n 17 

1 11111111 I 1 


22222222222 

33333333333 

44444444444 

S55555555S5 

6666GGG6666 

77777777777 

lllllllltll 


9 9 9 9 9 9 9 9 9 9|9 

miiim 

KSOI 


I MR 




2222222222 

3333333333 

4444444444 

555555S555 

6666666666 

7777777777 

iiiimii 1 


999999999B 

« is 14 is n n « a job 


1. 

00000000 do 


mtiiiiii 

9999999999 

RRNIiRIIRRN n 


100 . 

ooooooueolo 


31 C 33 >4 75 3G JJ J* JI 


11111111111111111111 

22222222222222222222 

33333333333333333333 

44444444444444444444 

SS5S55S555555S555S55 

66666666666666666666 

77777777777777777777 


1 1 n 1111 ik 

99999999 M9 

BHM»anMn« 4 n 


10 . 

00100800 0|0 
:uuu««<i < 1.1 Jsi 

I I 1 1 1 1 1 1 III 


2222222222 

3333333333 

4444444444 

S5S55SSS5S 

011(116(11 

7777777777 

Itlllllllt 


9 9 9 99 9 9 9919 

C 41 44 « « *J 41 «t it 


.01 

00000000009000000000000000000 

” u* a a S7 a a m e a a u m m u a a m n n n m 75 n n n n m 
1 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 1 1 1 j 1 | | 1 | | j j j 


22222222222222222222222222222 

33333333333333333333333333333 

44444444444444444444444444444 

SSS5S55555SSS5SS5S55555555555 

(( 66666 ( 6 ( 66 ( 66666 ( 6666666(66 

77777777777777777777777777777 

lllllllllllllllllllllllllllll 




FIG. 2 
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PERMANENT MAGNET GENERATOR 


CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYNBOL 


A, a 


(U6) 

A 

c 

AC 

(66) 

A 

g 

GA 

(79) 

a p 

AP 

(128) 

A 

A 

( 501 ) 

a l 

Al 

(502) 

a 2 

A2 

(516) 

Aip 

Y 

(519) 

\ 

AX 

(520) 

Al 

Al 


B, b 


(15) 

by 

BV 

(20) 

B 

BK 

(22) 

b o 

BO 

(22) 

b l 

Bl 

(22) 

b 2 

B2 

(22) 

b 3 

B3 

CNJ 

C\J 

b s 

BS 

(57) 

b tm 

TM 


RA-12 




CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(57a) 

b tl/3 

SM 

(76) 

b P 

BP 

(91) 

B t 

TE 

(9U ) 

B 

c 

BX 

(95) 

B g 

FH 


C, c 


(32) 

C 

C 

(60) 

C x 

cx 

(71) 

C 1 

Cl 

(72) 

Cw 

cw 

(73) 

c 

CP 


p 


(7U) 

c 

CM 


m 


(75) 

C q 

CQ 

(508) 

c 

ACM 


D, d 


(11) 

d 

DI 

(11a) 


DR 

(12) 

D 

DU 

(35) 

db 

DB 


RA-13 




CALCULATION 

NUMBER 


ELECTRICAL 

symbol 


FORTRAN 

SYMBOL 




(3) 

E 

EE 

(b) 

E PH 

EP 

(55) 

Ep (top) 

ET 

(56) 

Ep (bot) 

EB 

(516) 

e nl 

ENL 

(525) 

®FL 

EL 


F, f 


(5a) 

f 

F 

(96) 

F 

g 

FH 

(183) 

F Sc W 

WF 


-Si-jBL 


(69) 

g e 

GE 

(59) 

g min 

GC 

(59g) 

Smax 

GP 


H, h 


(22) 

h 

o 

HO 

(22) 

h l 

HZ 


RA-lii 



CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 

(22) 

h 2 

HY 

(22) 

h 3 

HZ 

(22) 

h s 

HS 

(22) 

h t 

HT 

(22) 

h 

w 

HW 

(2U) 

h 

c 

HC 

(37) 

^st 

SH 

(38) 

h 'st 

SD 

(137) 

h bl 

H 

(76) 

\ 

HH 

(76) 

h f 

HF 

(519a) 

h 

AH 


I, i 


(8) 

I PH 

PI 

(2U5) 


PS 


K, k 


(2) 

^A 

f 

VA 

(16) 

Ki 

RK 

(19) 

k 

WL 

0*2) 

K sk 

FS 


RA-15 


CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 

(U3) 

K d 

DF 

(hh) 

fa 

CF 

(61) 

fa 

FF 

(63) 

k e 

EK 

(67) 

K s 

CC 


■fa. A. 


(13) 

1 

L 

(17) 

fa 

SS 

(36) 

fa2 

CE 

(U8) 

fa 

EL 

(U9) 

fa 

HM 

(76) 

fa 

PL 

(139) 

fa 

SB 


M, m 


(5) 

m 

PN 


N, n 


(li*) 

°v 

HV 

(30) 

n s 

SC 

Oh) 

N . 
st 

SN 


RA-16 




CALCULATION 

NUMBER 

iLlOTRlCAL 

SYMBOL 

FOk'l'h.-.: 

SYMBOL 

(3Ua) 

N» . 

*> w 

SftL 

(45) 

n e 

EC 

(136) 

"b 

BN 


P> P 


(6) 

P 

PX 

(9) 

P F 

PF 

(5oo) 

P 1 

PI 

(503) 

P 2 

P2 

(SOU) 

P 3 

P3 

(5o5) 

? 3l 

PS.1 

(506) 

P s2 

PS2 

(507) 

P 

m 

PM 

(509) 

p i 

pi: 

(510) 

p o 

PO 

(5U) 

P 

g 

PG 


Q, q 


(23) 

Q 

QQ 

(25) 

q 

QN 


FLA-17 


CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 


-SiJL- 


(7) 

r pm 

RPM 

(53) 

r sph 

RG 

(51*) 

r SPH (Hot) 

RP 


S, s 


(1*7) 

S s 

S 


JjL v 


(11*5) 

v r 

VR 


W, w 


(181*) 

W TNL 

WT 

(165) 

W L 

WQ 

(186) 

W NPL 

WN 

(153) 

Wdnl 

WD 

(21*2) 

W TFL 

ST . 

(21*3) 

W PFL 

PP 

(21*1*) 

W DFL 

DL 


RA-18 


CALCULATION 

NUMBER 

ELECTRICAL 

SYMBOL 

FORTRAN 

SYMBOL 


X, x 


(50) 

X 

00 

o 

o 

T1 

(129) 

X 

XR 

(13C) 

X e 

XL 

042) 

Vc 

T3 

(163) 

X Dd 

XI 

(165) 

X Dq 

X2 

(166) 

x ’eu 

XJ 

(167) 

X ’d 

xs 

(168) 

X, ’d 

XX 

(170) 

x 2 

XN 

(523) 

X d olunu 

XD 


Y, y 


(31) 

y 

YY 


cr_ 


(26) 

T s 

1 

TS 

(27) 

^sl/3 

TT 

(40) 


SK 

(41) 

T P 

TP 

vl40) 

t; 

TB 


RA-19 



CALCULATION 

ELECTRICAL 

FORTRAN 

NUMBER 

SYMBOL 

SYMBOL 


A 


(62) 

Ai 

PC 

(64) 

Ae 

EW 

(70c) 

Aa 

AG 

(80b) 

Asl 

SL 

(81b) 

Ati 

TL 

(82b) 

A el 

ES 


-JL- 


(88) 

< h 

TG 

(92) 

h 

FQ 




(5D 

A 

RS 

(14D 

/d 

RE 




(77) 


PE 


■O' 


(198a) 

re- 

AN 


RA-20 


CALCULATION 

NUMBER 


ELECTRICAL 

SYMBOL 


FORTRAN 

SYMBOL 


K 


(167) 

K 1 

D1 

(168) 

k 2 

D2 

(189) 

K 3 

D3 

(190) 

K U 

DU 

(191) 

% 

D5 

(192) 

k 6 

D6 


RA-21 


C PASS 1 PERMANENT MAGNET 

DIMENSION DA(8),DX(6),DY(8),DZ(8) 

1 FORMAT(Ell.5,Ell.5,El 1.5,Ell.5,El 1.5,El 1.5) 

2 FORMAT(F7.0,F7.0,F7.0,F7.0,F7.0,F7.0,F7o0,F7.0,F7.0,F7.0) 
33 READ2,VA,EE,EP,PN,F,PX,RPM,PI,PF,DI 

READ2,DU,CL,HV,BV,SF,WL,BK,ZZ,BO,Bl 
READ2,B2,B3,BS,HO,HX,HY,HZ,HS,HT,HW 
READ2,QQ,W,RF,SC,YY,C,DW,SN,SN1,DW1 
READ2,DB,CE,SH,SD,PBA,SK,T1,RS p GC,GP 
READ2,C1,CW,CP,EL,CM,CQ,BH p BP,HH,HF 
READ2,PL,ALN,PE,DR,RK,WR,D1 p WO,HD,DD 
READ2,H,B,BN,SB,TB,RE,T3,WF,ACM,AH 
SS=SF*(CL—HV'^BV) 

HC=(DU-DI-2.0*HS)*0.5 

ZY=0.7*HS 

IF(HC-7.Y) 33,33,5 

5 QN=QQ/(PX*PN) 

TS=3.142*01/OQ 

IF(ZZ-4.0)29,30,29 

29 TT=(0.667*HS+DI )*3.142/QQ 
GO TO 31 

30 TT=3.1416*(DI+2.*H0+1.32*BS)/QQ 

9 

31 IF(ZZ-1.0)6,6,7 

6 BO=BS 

CC= (5.0*GC+BS)*TS/((5.0*GC+BS)*TS-BS*BS) 

GO TO 8 

7 QC=(4.44*GC+0.75*B0)*TS 
COQC/(QC-BO*BO) 

8 CS=YY/(PN*QN) 


RA-22 



TP-3.142*01/PX 
I F (SK)32,32,92 
32 FS=1.0 
GO TO 34 

92 FS=SIN(1.571*SK/TP)*TP/(1 0 571*SK) 

34 I F (PBA-60.) 9, 9,1 0 
9 D=1.0 

GO TO 95 

10 D=2.0 
95 l-QN 

U= I 

IF(QN-li)36,3b,35 

35 U=PX*PI'I 
XX=U 
N=U 

DO 11 K= 1, N 
Z=U/XX 
I =Z 

z 1 = 1 

IF(Z-Z1)12,12,11 
12 ZY=QQ/XX 

i=zy 

Z1 = I 

IF(ZY-Z1)37,37,11 

11 XX=XX-1. 

36 ZY=QN 

37 DF=SIN(1.571*D/PN)/(ZY*D*S!N(1.571/(PN*ZY))) 
CF=SIN(YY*1.571/(PN*QN)) 

EC=QQ*SC*CF*FS/C 


RA-23 



DT=DW1 

IF(DT) 13,13,14 

13 AC= 0 o 785*DW*DW*SN1 
GO TO 28 

14 ZY-0.0 

DA(1)=0.05 
DA(2)=0.072 
DA(3) = 0.1 25 
DA(4) = 0.165 
DA(5)=0.225 
DA(6)=0.438 
DA(7)=0.688 
DA(8)=1.5 
DX(1)=0.000124 
DX(2)=»0.00021 
DX(3)“0.00021 
DX(4)«0.00084 
DX(5)=0.00189 
0X(6)=0.00189 
DY(1)=0.000124 
DY(2)*=0.000124 
DY(3)=0.00084 
DY(4)=0.00084 
DY(5)=0.00189 
DY(6)=0.00335 
DY(7)=0.00754 
DY(8)=0.03020 
DZ(1)=0.000124 
DZ(2)=0.000124 

RA-24 


I 



DZ(3)=0.000124 
0Z(4)=0.00335 
DZ (5) = 0.00335 
DZ(6)=0.00754. 

DZ(7)=0.0134 
DZ(8)=0.0302 

93 IF(DT-.05)94,94,15 

15 JA=0 
JB=0 
JC=0 
JD=0 

16 JA=JA+1 
JB=JB+1 
JOJC+1 
JD=JD+1 

IF(DT-DA(JA))17,17,16 

94 D=0 

GO TO 23 

17 IF(DW-0.188)18,18,19 

18 CY=DX(JB—1) 

CZ=DX(JB) 

GO TO 22 

19 IF(0W-0.75)20,20,21 

20 CY*=DY( JC—1) 

CZ-DY(JC) 

GO TO 22 

21 CY-DZ(JD-I) 

CZ**DZ (JD) 

22 D»CY+(CZ-CY)*(DT-DA(JA-1))/(DA(JA)-DA(JA-1)) 
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23 AO(DT*DW-D)*SU1 
28 S=PI/(C*AC3 

PUNCH1 f VA,EE,EP p PN,F,PX 
PUNCH 1,RPM,PI,PF,D!,AC,S 
PUNCH 1 ,DU,CL,SS,HC,SF ,0N 
PUNCH1,WL,BK,Z7.,B0,B1,B2 
PUNCH1,B3,BS,H0,HX,HY,HZ 
PUNCH1,HS,HT,HW,00,W,RF 
PUNCH1,SC,YY,C,TS,SN,DB 
PUNCH1,CE,SH,S0,TT,SK,T1 
PUNCH1,RS,GC,GP,C1,CW,CP 
PUNCH 1,EL,CM,CO,BH,BP,HH 
PUNCH 1,HF,PL,ALN,PE,DR,RK 
PUNCH!,CC,WR,D1,WO,HD,DD 
PUNCH 1,H,B,BN,SB,TB,RE 
PUNCH1,T3,WF,CS,FS,TP,DF 
PUNCH1,CF,EC,ACM,AH,PBA 
PAUSE 
END 
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C PASS 2 PERMANENT MAGNET 

1 F0RMAT(E11.5,El 1.5,El 1.5,El 1.5,El 1.5,El 1.5) 
3 FORMAT(,9X F12.5,2X F12.5) 

READ1, VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,DI,AC,S 
READ 1, DU,CL,SS,HC,SF,ON 
READ 1, WL,BK,ZZ,B0,B1,B2 
READ1, B3,BS,H0,HX,HY,HZ 
READ1, MS,HT,HW,QO,W,RF 
READ 1, SC,YY,C,TS,SN,DB 
READ1, CE,SH,SD,TT,SK,T1 
READ1, RS,GC,GP,C1,CW,CP 
READ 1, EL,CM,CQ,BH,BP,HH 
READ 1,HF,PL,ALN,PE,DR,RK 
READ1, CC,WR,D1,WO,HD,DD 
READ 1, H,B,BN,SB,TB,RE 
READ 1,T3,WF,CS,FS,TP,DF 
READ 1,CF,EC,ACM,AH,PBA 
GA*=3.142*01 *CL 
AG=6.38*DI/(PX*GC*CC) 

GE=CC*GC 

IF(C1) 44,43,44 

43 C1=(0.649*L0G(PE)+1.359) *((GC/GP)**0.352) 

44 IF(CW)45,45,46 

45 CW=0.707*EE*C1*DF/(EP*PN) 

46 TG=6000000.0*EE/(CW*EC*RPM) 

BG=TG/GA 

IF(CP)47,47,48 

47 CP=(GC/GP)**0.41*PE*(L0G(GC/TP)*.0378+1.191) 
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48 FQ=TG*CP/PX 
IF(ZZ-3.0)49,50,51 

49 SM=TT-BS 
GO TO 53 

50 SM=(3.1416*(DI+2.*HS)/QQ)-B3 
GO TO 53 

51 IF(ZZ—4.0)50,52,49 

52 SM=TT-.94*BS 

53 TE=TG/(00*SS*SM) 

BX=0.5*F0_/(HC*SS) 

IF(EL) 54,54,62 

54 IF(RF) 55,55,61 

55 IF(PX—2.0) 56,56,57 

56 U=1.3 

GO TO 60 

57 IF(PX—4.0) 58,58,59 

58 U=1.5 

GO TO 60 

59 U=1.7 

60 EL=3.l42*U*YY*(DI+HS)/QQ+0.5 
GO TO 62 

61 EL=2.0*CE+(3.142*(0.5*HX+DB))+(YY*TS*TS/(S0RT(TS*TS-BS*BS))) 

62 HM=CL+EL 

RY=SC*QQ*HM/(PN*AC*C*C) 

RX=RS*0.000001 

RB=(T1+234.5)*0.00394*RX 

RG=>RX*RY 

RP=»RB*RY 

IF(SH)37,38,40 
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38 ET=1 
EB=1 

GO TO 39 

40 AA=0.584+(SN*SN-1.0)*0.0625*(SD*CL/(SH*HM))*2.0 
AB=(SH*SC*F*AC/(BS*RB*1000000.0))**2.0 
ET=AA*AB*0.00335+1.0 
EB=ET-0.0016 8*AB 

39 IF(CM)63,63,64 

63 AA“S1N(3.142*PE) 

AB=SIN(1.571*PE)*4.0 
CM=(3.142*PE+AA)/AB 

64 A=PI*SC*CF/(C*TS) 

PRINT3,SS,CC,HC,GA,TS,AG,TT,GE,FS,C1,DF,CW,CF,CP,EC,EL,AC,CM 

PUNCH 1,VA,EE,EP,PN,F,PX 

PUNCH1,RPM,PI,PF,DI, ACM 

PUNCH 1,DU,CL,SS,HC,SF,QN 

PUNCH1,WL,BK,ZZ,BO,B1,B2 

PUNCH1,B3,BS,H0,HX,HY,HZ 

PUNCH 1,HS,HT,HW,QO,W,GE 

PUNCH1,SC,YY,C,TS,BG,TG 

PUNCH1,FQ,TE,BX,TT,HM,SM 

PUNCH1,RG,GC,RP,C1,CW,CP 

PUNCH1,EL,CH,CQ,BH,BP,HH 

PUNCH1,HF,ALN,PL,PE,DR,RK 

PUNCH1,CC,WR,D1,WO,HD,DD 

PUNCH1,H,B,BN,SB,TB,RE 

PUNCH1,T3,WF,CS,ET,AH,EC 

PUNCH!,TP,DF,CF,EB,AC,S 

PUNCH!,AG,A,SM,PBA 
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PAUSE 


END 
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PASS 3 PERMANENT MAGNET 

1 F0RMAT(E11.5,El 1.5,El 1.5,Ell.5, El 1.5, El 1.5) 

3 FORMAT(9X F12.5,2X F12.5) 

READ 1, VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,DI,ACM 
READ 1, DU,CL,SS,HC,SF,QN 
READ1, WL,BK,ZZ,BO,B1,B2 
READI, B3,BS,HO,HX,HY,HZ 
READ 1, HS,HT,HW,QO,W,GE 
READ1, SC,YY,C,TS,BG,TG 
READ1, FO,TE,BX,TT,HM,SM 
READ1, RG,GC,RP,C1,CW,CP 
READ 1, EL,CM,CQ,BH,BP,HH 
READ 1,HF,ALN,PL,PE,DR,RK 
READ1, CC,WR,D1,WO,HD,DD 
READ 1, H,B,BN,SB,TB,RE 
READ 1,T3,WF,CS,ET,AH,EC 
READ1, TP,DF,CF,EB,AC,S 
READ 1,AG,A,SM,PBA 
IF (CO)69,69,70 

69 AA=1.571*PE 
AB=3.14l6*PE 

* 

CQ=(0.5*C0S(AA)+AB-SIN(AB))/(4.0*SIN(AA)) 

70 XR=.0707*A*DF/(C1*BG) 

IF(ZZ-5.)350,351,350 

351 FF=1.0 
GO TO 75 

350 IF(PBA-60.)352,353,352 

353 IF (CS-. 667 ) 354 ,355,355 



355 D=„75 
Z=. 25 

GO TO 74 
354 D=1.5 
Z=-.25 
GO TO 74 

352 IF(CS-.667)356,357,357 
357 FF=.75 
GO TO 75 

356 D=1.2 
Z=-.05 

74 FF»D*CS+Z 

75 CX-FF/(CF*CF*DF*DF) 

Z=CX*20.0/(PN*QN) 

BT=3.142*01/QQ-BO 
ZA=BT*BT/(16.0*TS*GC) 

ZB=0.35*BT/TS 

ZC=HO/BO 

ZD=HX*0.333/BS 

ZE-HY/BS 

IF(ZZ—2.0) 76,77,78 

76 PC=Z*(7.E+ZD+ZA+ZB) 

GO TO 82 

77 PC=Z*(ZC+(2.0*HT/(BO+BS))+(HW/BS)+ZD+ZA+ZB) 

GO TO 82 

78 IF(ZZ-4.0) 79,80,81 

79 PC-Z*(ZC+(2.0*HT/(B0+B1))+(2.0*HW/(B1+B2))+(HX*0.333/B2)+ZA+ZB) 
GO TO 82 

80 PC=Z*(ZC+0.62) 
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GO TO 82 

81 PC=Z*(ZE+ZD+(0.5*GC/TS)+(0.25*TS/GC)+0.6) 

82 EK=EL/(10.0**(0.103*YY*TS+0.402)) 

IF(DI-8.0) 83,83,84 

83 EK=SQRT(EK) 

84 ZF=.612*L0G(10.0*CS) 

EW=6.28*EK*ZF*(TP**(0.62-(0.228*L0G(ZF))))/(CL*DF*DF) 
IF(PN-3.0)85,86,86 

85 ZC=0.1*DI*SIH(3.0*YY/(PN*QN))*1.57/(PX*GE*CF) 

GO TO 87 

86 ZC=0.0 

87 XL=(PC+EW+ZC)*XR 
XD=XR*AG*C1*CM 
XQ=XR*CO*AG 

WC=0.3 21*SC*QQ*AC*HM 
PRINT3,S,C0,HM,A,RG,XR,RP,XL 
PUNCH 1,VA,EE,EP,PN,F,PX 
PUNCH1,RPM,PI,PF,DI,EC,PE 
PUNCH1,0U,CL,SS,HC,PC,QN 
PUNCH1,WL,BK,ZZ,BO,XD,AH 
PUNCH1,XR,BS,XL,HX,HY,HZ 
PUNCH1,HS,WC,AC,OQ,W,GE 
PUNCH1,SC,YY,C,TS,BG,TG 
PUNCH1,FO,TE,BX,TT,EW,AG 
PUNCH 1,RG,GC,RP,C1,TP,CP 
PUNCH 1,DF,CM,CF,BH,BP,HH 
PUNCH1,HF,PL,ALN,EB,DR,RK 
PUNCH1,CC,WR,D1,WO,HD,DD 
PUNCH 1,H,B,BN,SB,TB,RE 
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PUNCH1,T3,ACM,WF,CS 
PUNCH 1,ET,SM 
PAUSE 
END 
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PASS 4 PERMANENT MAGNET 
3 FORMAT(9X F12.5,2 X F12.5) 

1 FORMAT(E11.5,E1t.5 f El1.5,E11.5 p E11„5,El1 0 5) 
READ 1, VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,DI,EC,PE 
READ 1, DU,CL,SS,HC,PC,QN 
READ 1,WL,BK,Z7.,BO,XD,AH 
READ 1, XR,BS,XL,HX,HY,HZ 
READ 1, HS,WC,AC,QO,W,GE 
READ 1, SC,YY,C,TS,BG,TG 
READ 1, FO,TE,BX,TT,EW,AG 
READ 1, RG,GC,RP,C1,TP,CP 
READ 1, DF,CM,CF,BH,BP,HH 
READ 1,HF,PL,ALN,EB,DR,RK 
READ 1, CC,WR,D1,WO,HD,DD 
READ 1, H,B,BN,SB,TB,RE 
READ 1,T3,ACM,WF,CS 
READ 1,ET,SM 
ZA=3.1416*(DI+HS)/()Q 
IF(ZZ—3.0) 08,89,88 

88 TM=ZA-BS 
GO TO 90 

t 

89 TM=(3.1416*(DI+2.*HS)/QQ)-B3 

90 WI=(TM*Q0*SS*HS+(DU-HC)*3.142*HC*SS)*0.283 
AN=0.0 

100 AN=AN+0.005 
AL=COS(AN) 

IF(PF-AL) 100,100,101 

101 VR=0.262*DR*RPM 
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AP=BP*PL*ACM 
FH=BG*GE/0.00319 
IF (BN) 3 07,306,3 07 

306 X1=0 
P2=0 
X2=0 

GO TO 308 

307 IF(DD) 1 03,1 03,1 02 

102 ZG=0.62 
GO TO 104 

103 ZG=0.333*H/B 

104 BD=(HD/W0+ZG+0.5)*6.38 
BE=(BH-(BN-1.0)*TB)*2.127/GE 

P1=(BD+BE)*RL*C0S((BN-1,0)*TB*1.572/TP)/(BD+BE+RL) 
X1=XR*P1 

P2*=(HD /WO+Z G+0.5+G C /T B ) *2 0 . 0*T B /T P 
X2=XR*P2 

308 XU=XL+XF 
XS«=0.88*XIJ 
IF(BN)105,105,106 

105 XX=XS 

GO TO 107 

106 XX=XL+X1 

107 XN=(XX+XY)*0.5 

110 1F(WF)111,111,112 

111 WF=DR**2.5*(RPM**1,5)*PL*0.00000252 

112 WQ=(DU-HC)*1.42*HC*SS*(BX/BK)**2.0*WL 
WT=(SM)*QQ*SS*HS*0.453*(TE/BK)**2.0*WL 
PUNCH1,VA,EE,EP,PN,F,PX 
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PUNCH1,RPM,PI ,PF,DI ,PC,V/C 
PUNCH1,DU,CL,SS,HC,PC,0N 
PUNCH1 ,WL,BK,7.Z,B0,XD,WI 
PUNCH1,XR,BS,HX,HY,HZ,EC 
PUNCH1,HS,AC,OQ,W,GE,HD 
PUNCH!,SC,YY,C,TS,BG,TG 
PUNCH1,FQ,TE,BX,TT,EW,AG 
PUNCH1,RG,GC,RP,C1,AP,P2 
PUNCH1,OF,CF,FH,BP,HH,AL 
PUNCH 1,HF,ALN,PL,EB,DR,RK 
PUNCH1,CC,WR,D1,WO,HD,DD 
PUNCH1,H,B,BN,SB,TB,RE 
PUNCH1,T3,VR,WT,ACM,AH 
PUNCH 1,W0,WF,CS,ET,TP,X1 
PUNCH1,XX,XN,XL,CM,CP,PE 
PUNCH1,X2,XU,XS,AN 
PAUSE 
END 
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C PASS 5 PERMANENT MAGNET 

1 FORMAT(El1.5,E11.5,E11.5,E11.5,E11«5,E11o5) 
3 FORMAT(9X F12.5.2X FI2.5) 

READ1, VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,DI,PC,WC 
READl, DU,CL,SS,HC,PC,QN 
READ 1,V/L,BK,ZZ,BO,XD,WI 
READ 1,XR,BS,HX,HY,HZ,EC 
READ 1,HS,AC p QO,W,GE,HD 
READ 1, SC,YY,C,TS,BG,TG 
READl, FO_,TE,BX,TT,EW,AG 
READl, RG,GC,RP,C1,AP,P2 
READ 1,DF,CF,FH,BP,HH,AL 
READ 1,HF,ALN,PL,EB,DR,RK 
READ 1, CC,WR,D1,W0,HD,DD 
READ1, H,B,BN,SB,TB,RE 
READ 1 ,T3,VR,WT,ACM,AH 
READ 1,WQ,VVF,CS,ET,TP,X1 
READ 1,XX,XN,XL,CM,CP,PE 
READ 1,X2,XU,XS,AN 
GT=BO/GC 

IF(GT-1.0)304,304,303 

0 

304 AA=2.6 
GO TO 115 

303 IF(GT-3.75)113,114,114 

113 AA=10.0**0.178/((GT-1.0)**0.334) 

GO TO 115 

114 AA=10.0**0.11/((GT-1,0)**0.174) 

115 GF=AA*PI*SC/(C*FH) 
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305 IF(SC-1.0)121,121,122 

120 A5=0.0 

GO TO 129 

121 AX=1 .0 
AY=1.0 

GO TO 125 

122 AX=3.0*YY/(PN*0_N)-2.0 
IF(CS- 0 .667)123,124,124 

123 AY=1.5*YY/(PN*QN)-0.25 
GO TO 125 

124 AY=.75*YY/(PN*QN)+0.25 

125 A3=AX*P2/AY 

A4»0.07*AX*AG/(CF*CF) 

I F (AX ) 12 0, 12 0 , 126 

126 1F(BN)127,127,128 

127 A5=A4 

GO TO 129 

128 A5=(A4+A3)/(A3*A4) 

129 I F(W) 130,130,131 

130 X0=0.0 

GO TO 132 

131 AA=(3. 0 *HZ+HX)* 1 . 667 /(PN*QN*CF*CF*DF*DF*BS) 
X0=((PC+A5)*AX/AY+AA+0.2*EW)*XR 

132 D2=BG**2. 5*0.000061 
D3=(0.0167*QQ*RPM)**1.65*0.000015147 
IF(TS— 0 .9) 133,133,134 

133 D4=TS**1.285*0 .81 
GO TO 137 

134 IF(TS—2.0) 135,135,136 
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135 D4=TS**1.145*0.79 
GO TO 137 

136 d4=TS**0.79*0.92 

137 D7“B0/GC 

IF(D7—1.7) 138,138,139 

138 05=07**2.31*0.3 
GO TO 144 

139 IF(07-3.0) 140,140,141 

140 D5=D7**2.0*0.35 
GO TO 144 

141 IF(D7-5.0) 142,142,143 

142 D5“D7**1.4*0.625 
GO TO 144 

143 D5=D7**0.965*1.38 

144 D6=10.0**(0.932*01-1.606) 

BA**3.142*01 *CL 

WN*»D 1 *D2*D3*D4*05*D6*BA 
UY«(SL+ES+TL)*ALN*0.00638 
VT=0. 

AA=WO/(GC*CC) 

IF(AA)148,147,148 
148 IF(AA-0.65)145,145,146 

145 VT=L0G(10.0*AA)*(-0.242)+0.59 
GO TO 147 

146 VT=0.327-(AA*0.266) 

147 UZ=(DU-HC)*0.7850/PX 
EZ=(ET+EB)*0.5-1.0 
AA=PN*PI*PI 
PU=AA*RG 
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WRM=BP*PL*HF*.283*PX 
PV=AA*RP 

VV=EP*PI*PF*.003 
PUNCH1,VA,EE,EP,PN,F,PX 
PUNCH1,RPM,PI,PF,TB,BO,GC 
PUNCH1,HH,HF,ALN,DR,SB,RE 
PUNCH1,T3,W0,DD,H,BN,GF 
PUNCH1,VT,TS,CC,BG,AP,FQ 
PUNCH1,TE,BX,FH,W0.,WT,AN 
PUNCH1,RP,WF,HS,B,GE,BP 
PUNCH1,WN,UZ,EZ,PU,VV,XD 
PUNCH1,PV,HD,QN,XO,TG,VR 
PUNCH1,PC,WC,WR,WI,TP,DI 
PUNCH1,XI,X2,XU,XS,ET,EB 
PUNCH1,XX,PL,ACM,AH,EW,CM 
PUNCH1,PE,AG,XN,CP,CL,EC 
PUNCH1,XL,WRM,DF 
PAUSE 
END 
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PASS 6 PERMANENT MAGNET 
DIMENSION Al (30) 

1 FORMAT (Ell .5, El 1.5, El 1.5,EH.5,Ell.5,Ell. 5) 

888 FORMAT (F10.0,F10.0,F10.0,F10,0,F10.0,F10.0) 

K=1 

823 READ888,AI(K),AI(K+1),AI(K+2),AI(K+3),AI(K+4),AI(K+5) 
K=K+6 

IF(K-29)823,199,199 
199 READ 1,VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,TB,BO,GC 
READ 1,HH,HF,ALN,DR,SB,RE 
READ 1,T3,W0,DD,H,BN,GF 
READ 1,VT,TS,CC,BG,AP,FQ 
READ 1,TE,BX,FH,WQ,WT,AN 
READ 1,RP,WF,HS,B,GE,BP 
READ 1,WN,UZ,EZ,PU,VV,XD 
READ 1,PV,HD,QN,X0,TG,VR 
READ 1,PC,WC,WR,WI,TP,DI 
READ 1,X1,X2,XU,XS,ET,EB 
READ 1,XX,PL,ACM,AH,EW,CM 
READ 1,PE,AG,XN,CP,CL,EC 
READ 1,XL,WRM,DF 

t 

PD=FQ/AP 

Z=C0S(3.1416/PX) 

TAN=SIN(3.1416/PX)/Z 

A2=((DI/2.-G C-HH-HF)*TAN-BP/2.)*Z*2. 

Al=((DI/2.-GC-HH)*TAN-BP/2.)*Z*2. 

IF(A2-.2)401,401,400 
401 P1=3.19*PX*PL/3.1416 
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GO TO 402 

400 P1=PL*PX*(1.-(A2/(A1-A2))*L0G(1.+(A1-A2)/A2))*3.19 

402 P2=1,66*(!.+1,23*L0G(1./(1.-PE)))*PL 

P3=HF*1,66*(1.+1.23*L0G(1,+2./((A1+A2)/BP))) 
PS1=6.38*HH*PL/(TP-BP) 

PS2=2.*HH*P2/Pl_ 

PM=PL*BP*ACM/(2.*HF) 

PI I“PS1+PS2+P1+P3 
P0=PII+P2 
PG=.785*AG*CP*CL 
PW*=P I I+PG 
X=P0/PM 
NA= 1 

GO TO 802 

805 ENL=.001*(EP*Y/PD)*(P0/PM+AH)*(PW/PM-PII/PM)/(PW/PM+AH) 
PUNCH!,VA,EE,EP,PN,F,PX 
PUNCH1,RPM,PI,PF,TB,B0,GC 
PUNCH1,HH,HF,ALN,DR,SB,RE 
PUNCH1,T3,W0,DD,H,BN,GF 
PUNCH 1,VT,TS,CC,BG,F0.,Y 
PUNCH1,TE,BX,FH,WQ,WT,AN 
PUNCH1,WF,HS,B,GE,BP,PW 

» 

PUNCH !,V/N,UZ,EZ,PU,VV, PD 
PUNCH1,PV,HD,QN,XO,TG,VR 
PUNCH!,PC,WC,WR,WI,TP,DI 
PUNCH1,X!,X2,XU,XS,ENL,PG 
PUNCH1,XX,PL,AH,EW,RP,PM 
PUNCH1,PE,AG,XN,CP,CL,EC 
PUNCH!,XL,WRM,DF,CM,ET,EB 
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PUNCH1,F0,BG,TG,BX,TE 
PUNCH1,P I I»P0 
PAUSE 


802 IF(AI(NA)-X)830,831,831 
831 NA**NA+3 

835 IF(AI (NA)-X)833,834,834 

833 NA=NA+2 
GO TO 835 

834 AA=AI(NA) 

BB1=AI(NA-2) 

DC*=AI (NA+1) 

0=AI(NA-1) 

XXX«(AA-BB1) /(.4343*(LOG(DC)-LOG(D+.OOOl))) 
Y-AA-XXX*.43 43*L0G(DC) 
Y=EXP(2.306*(X-Y)/XXX) 

GO TO 805 
830 PR1NT850 

850 FORMAT(17HMACHINE SATURATED) 

PAUSE 

END 
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C PASS 7 PERMANENT MAGNET 
DIMENSION EL( 6 ) 

3 FORMAT(9X F12.5.2X F12.5) 

1 FORMAT(E11.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 

READ 1,VA,EE,EP,PN,F,PX 
READ 1,RPM,PI,PF,TB,BO,GC 
READ 1,HH,HF,ALN,DR,SB,RE 
READ 1,T3,WO,DD f H,BN,GF 
READ 1,VT,TS,CC,BG,FO_,Y 
READ 1,TE,BX,FH,WO,WT,AN 
READ 1,WF,HS,B,GE,BP,PW 
READ 1,WN,UZ,EZ,PIJ,VV,PD 
READ 1,PV,HD,QN,XO,TG,VR 
READ 1,PC,WC,WR,WI,TP,DI 
READ 1,X1,X2,XU,XS,ENL,PG 
READ 1,XX,PL,AH,EW,RP,PM 
READ 1,PE,AG,XN,CP,CL,EC 
READ1 ,XL,WRM,DF,CM,ET,EB 
READ 1,FO_,BG,TG,BX,TE 
READ 1 # PII,PO 
PW=PI I+PG 

AX=Y*(PO/PM+AH)/(PIl/PM+AH) 

Al=Y*(PO/PM+AH)/(PW/PM+AH) 

AZ=P I 

326 HD 1 =.45*A7.*EC*CM*DF/(PX*HF) 

AISC=(ENL*(AX-A1-HD1)/(AX-A1))/SQRT(RP**2+(XL*EP*.01/PI)**2) 
D=AISC-AZ 
IF(D)322,323,323 
322 D=-D 
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323 IF(D-.05)324,325,325 
325 AZ=(AlSC—AZ)/2.+AZ 

GO TO 326 

324 XD*=ENL/AI SC 
XDPU=XD*PI*100./EP 
TPF=SIN( AN)/COS(AN) 

SPF=SIN(1,5708-AN) 

CPF=COS(1,5708-AN) 

Z=. 25 

DO 327 K=1,6 
AZ=Z*PI 

X=(AZ*XD-AZ*RP*TPF) 

B=X*SPF/ENL 

B=ATAN(B/SQRT(1,-B**2)) 
EL(K)=F.NL*COS(B)-X*CPF-AZ*RP/PF 
327 Z=Z+.25 

PRINT3,ET,XDPU,F.B,X1,PC,X2,EW,XU,WC,XS 
PRINT3,WI,XX,TP,XN,PI I,XO,PO,TG,PM,FQ 
PRINT3,PG,BG,WR,TE,WRM,BX f VR,AISC 
PRINT321,F0 
PRINT3,PD 

321 FORMAT (/9X FI2.5) 

PUNCH1,F,PX,TB,B0,GC,SB 

PUNCH 1,RE,T3,W0,DD,H,BN 

PUNCH 1,TS,CC,BG,B,HD,ON 

PUNCH1,EZ,XDPU,WT,WF,WQ,VV 

PUNCH1,ENL,EL(1),EL(2),EL(3),EL(4),EL(5) 

punchi,el(6),pu,pv,gf,wn 

PAUSE 
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END 
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PASS 8 PERMANENT MAGNET GENERATOR 

DIMENSION PR(2),PS(2),G(2),DL(2),PP(2),EX(2),ST(2),VA(2) 
DIMENSION P(2),E(2),PM(2),SP(2),EL(6) 

1 FORMAT(El 1.5,Ell.5,Ell.5,Ell.5,Ell.5,Ell.5) 

3 FORMAT(9X,F12.5,2X,F12.5) 

READ 1,F,PX,TB,BO,GC,SB 

READ 1,RE,T3*V0,DD,H,BN 

READ 1,TS,CC,BG,B,HD,QN 

READ1 ,EZ,XDPU,WT,WF,WQ,VV 

READ 1,EWL,EL(1),EL(2),EL(3),EL(4),EL(5) 

READ 1,EL(6),PU,PV,GF,WN 
IF(RE)309,309,311 
309 WD=»0. 0 
WU«0.0 
GO TO 178 

311 FS1=2.0*0_N*PN*F 
FS2=2.0*FS1 
M=0 

150 IF(M-1)151,152,178 

151 RM=RE 

GO TO 153 

152 RM=RE*(T3+234.5)/254.5 

153 AA=(FS1/RM)**0.5*DD*0.32 
AB=(FS2/RM)**0.5*DD*0.32 
IF(AA-2.5) 160,160,161 

160 VI=1.0—0.15*AA+0.3*AA*AA 
GO TO 162 

161 V1=AA 

162 IF (AB—2.5) 163,163,164 
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163 V2=1.0-0.15*AB+0.3*AB*AB 
GO TO 165 

164 V2=AB 

165 IF(H-B) 167,166,167 

166 VC=0.75/VI 
GO TO 169 

167 IF(DD) 166,168,166 

168 VC=H/(3.0*B*V1) 

169 VS=HD/V'/O+VT +VC 
VG=TB/(CC*GC) 

Q1 «=1 .0—(1.0/((( B0*0. 5/GC )**2.0+1 ,0)**0.5) ) 
QZ*=BO/TS 

Q2=1.05*3 IN(0Z*2.844) 

IF(QZ-0.37)170,170,171 

170 Q3=0.46 
GO TO 172 

171 Q3“0.23*SIN(10.46*0Z-2.1)+0.23 

172 Q4»SIN(6.283*TB/TS-1.571)+1.0 
Q5»SIN(12.566*TB/TS-1.571)+1.0 
t F(H)173,173,174 

173 AB=0.785*DD*DD 
GO TO 175 

174 AB=H*Dt) 

175 W2-PX*BM*SB*RM*1.246/(AB*1000.) 
W3=(Q2/(2.0*VS+(VG/Q4)))**2.0*V1 
W5«(Q3/(2.0*VS+(VG/Q5)))**2.0*V2 
WD=(T S*BG*Q1 *CC)**2.0*W2*(W3+W5) 

M=H+1 

IF(^1)176,176,177 
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176 WU=WD 

177 GO TO 150 

178 G(1)=0 
G(2)=1 
PW=PIJ 
WW=WU 

DO 183 M=1,2 
UA=G(M) 

PS ( M) ■=PW*UA*UA 
X=V/F +V/0 

GM=(GF*UA)**2.+1. 

ST(M)=(2.0*(0.0027*XDPU*UA)**1.8+1.0)*WT 

VA(M)=VV*IJA 

DL(M)«GM *WW 

PP(M)-GM*WN 

EX(M)-EZ*PS(M) 

SP(M)=PP(M)+DL(M)+PS(M)+EX(M)+ST(M)+X 
P(M)=*(SP(M) / I 000. )+VA(M) 

IF(GM)185,184,185 

184 PM(M)=0 
E(H)=0 

GO TO 186 

185 PM(M)-(SP(M)/P(M))*.1 
E(M)=100,—PM(M) 

186 WW=WD 
183 PW=PV 

PR INT3,WF,WF 
PRINT3,ST(1),ST(2) 

PR I NT3 f V/0,WQ 
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PRINT3,PP(1),PP(2) 

PRINT3,RL(1),DL(2) 

PRINT3,PS(1),PS(2) 

PRINT3,EX(1),EX(2) 

PRINT3,SP(1),SP(2) 

PRINT3,VA(1),VA(2) 

PRINT3,P(1) f P(2) 

PRINT3,PM(1),PM(2) 

PR INT3,E(1),E(2) 

PRINT321.ENL 

PR INT4,EL(1),EL(2),EL(3),EL(4),EL(5),EL(6) 
4 FORMAT(9X,F12.5) 

321 FORMAT (/9X F12.5) 

PAUSE 

END 
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POLE-FACE LOSSES IN SOLID POLE GENERATORS 


IN BRIEF 

Pole-face losses in solid rotor generators can limit the rotational speed, the 
output, or both. The following article discusses the calculation of, and the 
design limits imposed by the pole face losses. 

GENERAL STATEMENT 

Paradoxically, the solid rotor machine is proposed for high speeds and high 
air-gap flux densities when it cannot operate at speeds or loadings comparable 
to those of the laminated pole generators, without the penalty of high losses in 
the solid pole faces. 

A solid-pole Lundell generator that operates over a speed range of two to one 
may require de-rating at the maximum speed because of pole-face losses. 

This de-rating is required because of the pole face losses under load - not the 
no load losses. 

When the wide speed range generator operates at its maximum speed, the gap 
density is reduced and the reduction of losses due to the reduced flux density 
is greater than the loss increase due to the increase in tooth ripple frequency. 
The result at no load is a reduction in pole face losses. At full load though, 
the same stator current flows and the armature mmf is the same as it was at 
the lower speed. 

Because of the doubled slot frequency, the pole face load losses at the top speed 
will be of the order of 3 times as great as at the 1/2 speed condition. The actual 
increase in pole face losses under load depends upon the slot opening, air gap 
length, ampere loading, and gap density. 

SOURCE OF THE FORMULAE 

T. Spooner and I. F. Kinnard "Surface Iron Losses with Reference to Laminated 
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Materials" Trans. AIEE, Vol. 43, 1924, pp 262-281. 

The above reference points out that the following factors influence pole face 
losses: 

1. Air gap induction Bg 

2. Field form C. 

b 

3. Ratio of slot width to single air gap — 

4. Tooth frequency f^. 

*5. Tooth pitchT' or slot pitch 

s 

6. Resistivity of material o°) 

7. Thickness of individual laminations t 

8. Hysteresis Coefficient f 

9. Insulation between laminations 

10. Effect of punching 

* We use different symbol than Spooner and Kinnard 
In the equation - 

Pole face losses = W g = Kj Kg Kg K 4 Kg K 0 (Bore area), the hysteresis co¬ 
efficient is not used, the interlamination resistance is assumed high, no burrs 
are assumed to be shorting the laminations and the effect of work-hardening 
is assumed to be removed by annealing. 

In the paper, Spooner and Kinnard summarize that as a rough approximation, 
surface losses vary as: 

h 1 to 2 

\2 /t \ 1 • 5 s , . 1 

(B g ) , (£ t ) , — (r s ) 

In the text of the paper, higher exponents are used and the constant K^ is ad¬ 
justed to give the required accuracy. The exponents used in the curves given 
in the design manuals of this report are about 30% higher than the rough ap¬ 
proximation given by Spooner and Kinnard, and account for interlaminated 
currents caused by imperfect insulation. 
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SOLID-POLE P.F. LOSSES 

For this study, tests were made on solid-pole generators using 4130 pole steel. 
Tests correlated within 5% of calculated losses when the constant Kj = 7.0 
was used. Load losses in the pole face are calculated just as for laminated 
poles. 

Note: If the rotor is not relieved at the leading and trailing edges of the pole, 
the losses will be higher than those calculated with Kj = 7.0. Such a condition 
might exist when 316 steel or Inconel X is welded between poles and the weld¬ 
ment is cylindrical. In such cases, the armature reaction flux causes high 
losses in the interpolar area. This discussion and calculation assumes that 
best conditions exist. 

NO-LOAD POLE FACE LOSSES 

The pole face losses are a function of the stator bore area, and when the stator 
tooth pitch is constant, the losses are a function of the rotor diameter. They 
are also functions of the gap density and under load they are functions of the 
stator loading. Since the RPM determines the tooth ripple frequency, the P.F. 
losses are a function of the RPM. 

The flux density and speed of a specific, solid pole generator, are limited by 
the amount of heat that can be dissipated from the surface of the poles so with¬ 
out considering the output of the generator, the effect of slot pitch and air gap 
density on pole face losses can be observed, while the generator operates at 
no-load. 

With a specified air gap density and specified slot pitch, the pole face losses 
can be made high enough to cause the machine to fail and this can happen with¬ 
out load being applied. It can happen to a laminated pole or to a solid pole 
generator. 

The no load pole face losses in a solid pole generator are of the order of six 
(6) times as great, for the same stator design, as those in the pole faces of a 
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laminated rotor. For this reason only the design limits for the solid pole gen¬ 
erators are investigated here. 

b 

No load pole face losses are = fn (Bg, —, f^, C j) where 

2 

Bg = air gap density, Kilolines/in 

b = slot opening, inches 
s 

g = air gap length, inches 
f^ = tooth ripple, cycles/sec 

'Y = tooth pitch on stator, inches 

C ^ = fundamental of field form (a ratio) 

For most of the following discussion, the following values are assumed fixed: 
b 

— = 2.0 

g 

'Y = .3” 
s 

Cj = 1.0 

Since = . 3, now the number of teeth is a function of d. 
s 

No. Teeth = ^ 

, _ 7Td RPM 

\ " .3 ‘ "W 

PFL at NL = K 1 K 2 K g K 4 K g K g (Bore Area) 

K x = 7.0 

2 

K« = fn (B ) = .8 for 45 Kilolines/in gap density 
& g 
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/. 63 ~ 

Kg = fn (f T ) 



K 4 = fn (r s ) = • 175 


b 

K k = fn (—) = 1.45 
5 g' 

Kg = fn (C^ = .22 

PFL at NL = (7. 0) . 8 (Kg) (. 175) 1. 45 (. 22) (Bore Area) 
= .313 Kg (Area) for Bg * 45 Kilolines/in 2 
= .235 Kg (Area) for Bg = 40 Kilolines/in 2 
= . 11 Kg (Area) for Bg = 30 Kilolines/in 2 


= . 039 Kg (Area) for Bg = 20 Kilolines/in 2 


Rotor 

Dia. 

In. 


c 

No-Load Pole Face Loss in Watt/in* 

II 

CO 

RPM 

_i 

Bg = 45 

Bg = 40 

Bg = 30 

Bg = 20 

4 

2,000 

.72 

.54 

.25 

.09 

If 

4, 000 

2.28 

1.7 

.8 

.284 

ft 

8, 000 

7.2 

5.4 

2.5 

.9 

tl 

12, 000 

14. 1 

10.6 

4.9 

1.76 

If 

24, 000 

44 

33 

15.4 

5.5 

If 

48, 000 

128 

96 

45 

16 

6 


1.4 

1.05 

.49 

. 17 

ft 

4, 000 

4.4 

3.3 

1.54 

.55 

ft 

8, 000 

14 

10.5 

4.9 

1.75 

ft 

12, 000 

27 

20.2 

9.5 

3.37 

ft 

24, 000 

86 

64.5 

30 

10.7 

ft 

48, 000 

272 

204 

95 

34 
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POLEFACE LOSSES IN A SOLID POLEFACE AT 

NO LOAD AS A FUNCTION OF SPEED WHEN 

THE GAP DENSITY B g = 45 KL/f N ^ 

-2,0 SLOT PITCH *X =* 30 
a S 



ROTOR SREED_, THOUSAND R PM 




CURVE SA-] 




POLE FACE LOSSES IN A SOLID POLE FACE 
AT NO LOAD AS A FUNCTION OF SPEED AND 


AIR GAP DENSITY FOR A 4.0 DIA. ROTOR. 



ROTOR SPEED, THOUSAND RPM 

'•-8 CURVE S 4-2 

I 







POLE FACE LOSSES IN A SOLID POLE FACE 
AS A FUNCTION OF SPEED 
FOR A 6.0" DIA. ROTOR. 


AT NO LOAD 

AND air-gap DENSITY 
3 = 2,0 SLOT PITCR, 2 / =.30 



ROTOR SPEED, THOUSAND RPM 

CURVE SA-3 
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DISSIPATING THE POLE FACE LOSSES 


If the pole face losses are to be dissipated to gas similar to air or if the ma¬ 
chine is force ventilated with air, the best coefficient that could be expected is 
0.10 watts/in 2 of surface/°C rise above the air. See curve5A4from Luke and 
curveSA5from A. D. Moore. 

2 

If 25 watts/in is to be dissipated from the rotor surface, the temperature rise 

25 o 

of the rotor surface in this case will be -jq = 250 C above the cooling air or 
gas. Similarly, if 50 watts/in 2 must be dissipated, then the temperature rise 
of the rotor will be 500°C above the cooling air. 

2 

If room temperature cooling air can be used in sufficient quantity, 40 watts/in 
might possibly be dissipated from the generator rotor surface. If a heat trans¬ 
fer coefficient of . 10 w/in 2 /°C can be obtained, the rotor temperature will be 
approximately 800°F or 425°C. 

A cooling coefficient more likely to be obtained is . 08 watts/in 2 /°C and this 
value has been used on curveSA^to show probable rotor temperatures associ¬ 
ated with various no-load pole face loss levels. 

EFFECT OF REDUCING GAP DENSITY 

If a high-speed generator is designed and the PF losses are too great at say 
45 kilolines/in air gap density, reducing the flux density will reduce the pole- 
face, no-load losses by the ratio of gap densities raised to the power 2.5. The 
factor for gap density is Kg- 


Bg 

K 2 

20 

. 10 

25 

. 18 

30 

.28 

40 

.6 

45 

.8 
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WATTS PER SQ. INCH PER DEGREE CENT 


SURFACE HEAT DISSIPATION FROM 
A GENERATOR ROTOR. 



AVERAGE AIR VELOCITY THROUGH 
MACHINE IN THOUSANDS OF FEET 
PER MINUTE. 

FROM LUKE: HEATING OF RAILWAY 

MOTORS. AIEE TRANS 1922 VOL 41 
PP 165-173. 
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CU RYE 5A-4- 



WATTS PER SQ. INCH PER DEGREE 


u 


WATTS DISSIPATED FROM THE ROTOR OF AN 
ELECTRICAL MACHINE WHEN FORCED-AIR 
COOLING IS USED. 

FROM A. D. MOORE: FUNDAMENTALS OF 
ELECTRICAL DESIGN. McGRAW HILL 1927 
PAGE 134. 



AVERAGE COOLING AIR VELOCITY THROUGH 
THE MACHINE AIR-GAP - THOUSANDS OF 
FEET PER MINUTE. 
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NO LOAD POLE FACE LOSSES FOR A SOL ID-ROTOR 
GENERATOR HAVING 6.0 DIA. ROTOR AND A TOOTH 
PITCH OF ,3'.' bS/g — 2.0 



3 4 6 8 10 20 30 40 60 80 


ROTOR SPEED-THOUSAND RPM 


CURVE 




EFFECT OF VARYING SLOT PITCH 


When the slot pitch is varied, the pole face loss varies as the slot ripple fre¬ 
quency variation raised to the 1. 65 power and varies as the tooth pitch variation 
raised to the 1. 26 power. The no-load pole-face loss will decrease as the tooth 
or slot pitch is increased. 

The net change in the loss as the slot pitch increases is approximately (ratio 
of slots)*'*’'* (ratio of slot pitches)*' . 

In the tabulation for diameter = 4.0", Bg = 40 Kl/in^, ^ = .3", RPM = 

2 s 

12000, the NLPFL = 10.6 watts/in of stator bore. 

If the slot pitch is increased to . 6", the NLPFL will be reduced to 8. 1 watts/ 
in^ of stator bore area. 


EFFECT OF SLOT OPENING 


Completely closing the slots can have only limited effect in reducing pole face 
losses since the bridge closing the slot must saturate before the machine can 
operate. The saturated bridge represents an effective opening in the slot and 

bo 

the ratio — is some definite value. In considering the effect of closing the slot 

S |j| 

with a bridge, a smaller effective opening can be assumed such as — = 1. 0. 

& 

The pole face loss at no load varies approximately as the ratio 



So, by reducing the slot opening to 1/2 its original width, the no-load losses 
can be reduced to about 1/2. 4 of the initial no-load value. This applies to the 
no load losses only and the pole-face losses under load may be much higher 
because of the reduced slot opening . Load losses in the pole-face are a function 
of the stator loading and of the slot opening to air gap ratio (bs/g). This can be 
expressed as a function of bs/g and a function of (Xad) . 
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LOAD LOSSES IN THE POLE FACE 


Under load, the pole face losses are increased. This increase is due to the 
armature mmf wave and resulting flux wave which has a saw tooth shape. 

The higher the current loading in the stator, the more pronounced the saw 
tooth shape of the flux wave becomes when the machine is fully loaded. This 
effect is a function of the ratio - 


series conductors/slot x I , K 
__ ph sc 

_ no load air gap ampere turns 


(NLPFL) 

b 


where K gc is a modifying factor describing the effect of — . The smaller the 
opening, the sharper the step or saw-tooth in the armature reaction flux wave 


E.I. Pollard treated the pole face load losses in his paper "Load Losses in 
Salient Pole Synchronous Machines", AIEE Trans., Vol. 54, 1935, pp 1332- 
1340. His work shows that if the slot opening to air gap ratio is made too small, 
the result may be an increase in total pole face losses. See curve SA-7 


In the rotor sizes 4. 0" dia. to 8. 0" dia., an air gap length of . 030" is reason¬ 
able to assume. If the following values are assumed: 

ampere loading = 900 ampere wires per inch of stator bore periphery, 

bs/g = 2.0, Bg = 40 Kl/in , then the air gap ampere turns are about 400 and 
from curve 


PF load losses 


K sc x ph N s 


L AGAT J 


NLPFL 


1.4 (270) 
3(50 


NLPFL 


PF Load Losses = . 91 x no load PF losses 

where N = conductors/slot 
s 


J ph N s = 900 x • 3 = ampere wire loading x slot pitch 
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LOAD LOSSES IN THE POLE FACE 
FROM E.l. POLLARD 1 . LOAD LOSSES IN 
SALIENT POLE SYNCHRONOUS MACHINES 
AIEE TRANS. VOL. 54 1935 PPI332 1340. 



16 


CURVE SA-7 




SOLID POLE FACE 

POLE FACE LOSSES AT NO LOAD AND AT FULL 
LOAD FOR 4.0" DIA. ROTOR AT VARIOUS GAP 
DENSITIES. A=900 t%^. = 2.o T =.Z 
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CURVE SA-8 



POLE FACE 

LOSSES AT 

NC 

LOAD AND AT FULL 

LO 

FOR 6.0" DIA. ROTOR 

AT 

VARIOUS GAP DENSIT 

IES 

A= 900 

= 2.0 


r-3 



LOSS LIMIT 




2.5 5 10 20 30 40 50 

ROTOR SPEED, THOUSAND RPM 

• a .i 8 CURVE SA —c 9 








POLEFACE LOSSES AT NO LOAD AND AT FULL LOAD 
FOR 8.0" DIA. ROTOR AT VARIOUS GAP DENSITIES 



1.5 a5 5 10 15 20 30 40 

ROTOR SPEED, THOUSAND RPM 

C U RVE SA-lO 
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POLE FACE LOSSES AT NO LOAD AND AT FULL 
LOAD FOR 12.0" DiA, ROTOR AT VARIOUS GAP 
DENSITIES, A=900 = 2.0 r -30 



ROTOR SPEED, THOUSAND RPM 


3A-20 


curve: s a-ii 






POLEFACE LOSSES AT NO LOAD AND AT FULL LOAD 
FOR ROTORS OF VARIOUS DIAMETERS = 30 Kl/IN 2 

A = 900 AV y/|N SLOT PITCH-.4 =20 



ROTOR SPEED., THOUSAND RPM 


CURVE S A.-/*. 








ROTOR DIA. VS ROTOR SPEED FOff SOLID POLE FACE ALTERNATORS 
LIMITED TO 20 WATTS/ IN 2 POLE FACE LOSSES BASED ON STATOR 



ROTOR DIAMETER, INCHE5 CURVE SA" 
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ROTOR DIAMETER^ INCHES CURVE 5A-14 


AGAT = 400 = ampere turns drop across main air gap 


K = loss factor from curve 
sc 


The pole-face load losses are added to the no-load pole-face losses. 


When the air gap density is decreased, the ampere turns drop across the air 

gap decreased in direct proportion and the pole-face load losses increase as 

the square of the change in gap density. For example, the factor just calcu- 

lated above for a 40 Kl/in gap density was . 91 (NLPFL). When the gap density 

2 

in the same machine is reduced to 30 Kl/in but the ampere loading is main- 

2 

tained at 900 AW/in , the loss eqn. becomes - 


PFLL = 


1.42 (270) 

m — 


NLPFL = 1.63 (NLPFL) 


4 

or PFLL = ( 3 ) (.91) NLPFL =1.63 (NLPFL) 


We have already discussed how the reduction of the air gap flux density re¬ 
duces the pole face no load loss per square inch of stator bore area by 


2 5 2 2 

and hi the case of a reduction from 40 Kl/in to 30 Kl/in the re- 


B 


2.5 


. gl/ 3 

duction in NLPFL is (^) = . 488. 


The full load PFL increased by (^) or 1. 78 so the total load loss was reduced 
from 1. 0 to 1. 78 (. 488) = . 87 a reduction of 13%. 


SLOTTING THE ROTOR TO REDUCE THE POLE-FACE LOSSES 

Figure 1 shows graphically what happens when a solid pole-face is slotted to 
reduce pole face losses, when the slot approximates 50% of the pole surface. 

Assuming that the air gap flux density is unity, the initial pole face density at 
its surface was also unity. 
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When the rotor was slotted to the extent that 50% of the surface was re¬ 
moved, the density in 1/2 the remaining surface or 1/4 of the initial sur¬ 
face became 1.74 times the initial value or 1.74 per unit. The surface 
density in the remaining 1/4 of the initial surface increased to 2 times the 
surface density or 2.0 per unit. The average increase in total pole face 
losses based upon the increase in gap density at the surface of the pole 
is 2.4 times. However, the increased air gap density in the region of 
the pole surface reduces the effect of the slot ripple and tends to reduce 
the load losses. These effects can nearly enough offset the effect of 
grooving to make the pole face losses about the same with or without 
grooving. 

Experiments made in support of this study showed insignificant results when 
the rotor was grooved approximately with .020-. 03 wide grooves, and only 
50% of the initial surface was left. 

Figure SA1 shows the results obtainable when the rotor grooving is about 
10 to 12% of the surface area. This amount of grooving does not change 
the rotor surface loss a significant amount. To realize a satisfactory re¬ 
duction in pole face losses, the lands formed by the grooving should be 
thin; .020" width is a good maximum for 400 cps generators. 

To satisfy the requirement for lamination-thickness lands, and only 12% 
loss of surface, the grooves should be only about .0025" wide and the re¬ 
quirement for such an extremely narrow slot imposes a difficult manufactur¬ 
ing problem. Some methods that can be used to cut the narrow grooves are: 
diamond sawing with narrow saw or wire, electric-discharge machining with 

small diameter wire used as electrodes, and electron-beam machining. 
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DISTRIBUTION OF SURFACE FLUX DENSITIES 
WHEN A ROTOR POLE-HEAD HAS 50% OF ITS 
SURFACE CUT AWAY BY GROOVING--ASSUME ? 

AN AIR-GAP DENSITY OF 50 KILOLINES PER INCH' 4 
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STATOR SURFACE 



WITH . 025" LANDS AND . 005" WIDE SLOTS, 
THE POLE SURFACE FLUX DENSITY IS 109% 
THAT OF A SMOOTH POLE FACE 


FIGURE. SA ~Z 



FLUX-PLOTTING 
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GRAPHICAL FLUX ANALYSIS 


Graphical flux analysis is the quickest and most direct solution to many field 
problems. Irregular or complex fields that yield slowly to mathematical anal¬ 
ysis, if at all, can be solved graphically. Manual plots are used to solve heat- 
flow, air-flow, dielectric field, and flux field problems. 

In this study, only flux field distributions are mapped, and only the simplest 
case is considered — that is, the case where the iron surfaces are equipoten- 
tial surfaces and the space between the iron surfaces has the permeability of 
air. The same total field potential exists across the air space regardless of 
any change in dimension or configuration. 

The permeability of air is constant, so if the field gradient per unit of linear 
measurement (ampere turns per inch) across the air space changes, the flux 
density must change by the same ratio. 

In making flux plots, the flux is considered to consist of tubes of flux. In a 
three dimensional field, a single tube is conceived of as looking like this: 



In a two dimensional field, the depth of the tube is constant and the tube looks 
like this : 



If the number of lines of flux in a tube is held constant, and the depth of the 
tube is constant, the sides of the tube of flux must converge or diverge in di¬ 
rect ratio to the change in field gradient. 

By choosing the scale of field gradient and flux density, the area of flux tube 
enclosed by the sides of the tube and the equipotential field gradient lines be¬ 
come a square. 
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The lines of flux must cross equipotential lines at right angles and enter the 
iron at right angles. The corners of the areas so defined are right angles and 
if the squares are successfully divided into fourths, the smaller squares be¬ 
come more closely true squares. In any of the squares, the two dividing lines 
used to divide the square into four squares will be closely equal in length. 



To start a map or flux plot, draw the area to be mapped as large in scale as 
practical. Ink the boundaries of the iron so erasures will not remove them. 
Also, ink in the lines of symmetry. 


SA-31 



Draw in the middle an equipotential line and then start the plot in the most 
irregular region or at a corner and a line of symmetry. 



Start the maps with large sweeping curves and close the map before trying to 
perfect the detailed areas. 

Use 2H pencil for easy erasures with minimum smudging. 

Make large squares and divide only where necessary to check map accuracy 
or to obtain accurate density. 

Intersections between flux lines and equipotential lines must be right angles 
at all times or the map can npt be made accurately. 

When the tubes of equal flux are divided, they are most easily divided into 
multiples of two. 

A map can not be forced. However, useable accuracy can be obtained without 
a precise map, and the operator should exercise his own judgment as to the 
accuracy required. 
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To move lines a small amount, try widening the line with pencil and removing 
the unwanted portion with a sharp-pointed eraser. 



To determine, by means of a freehand flux plot, the no load field form of an 
electromagnetic machine: 

1. Make the flux plot in the manner described in the literature, with 
all intersections of flux lines and equipotential lines at right ang¬ 
les and with all curvilinear squares capable of being divided into 
smaller, curvilinear squares. 
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2. The distance from centerline between poles to the centerline of 
the pole is 90 electrical degrees. Divide this arc on the stator 
surface into 10° increments. 



3. At any point on the stator surface, the distance from the surface 
to the first equipotential line is proportional to the flux density at 
that point. A plot of flux density, therefore, is a plot of the ratio 
of distance from the stator surface to the first equipotential line. 
Where the equipotential lines have been further divided, the dis¬ 
tance ratios increase proportionately. 


sa-34 


4. The maximum density (usually at the pole head centerline) can be 
used as a one per unit; 1. 0 or 100%. 



Cos 0° = 1.000 
Cos 10° = . 985 

Cos 20° = . 940 

Cos 30° = . 866 

Cos 40° = . 766 

Cos 50° = . 643 

Cos 60° = . 500 

Cos 70° = . 342 

Cos 80° = . 173 

Cos 90° = . 000 


Ref: Mathematics of Modern Engineering, Vol 1, pp 73-92, Doherty and 
Keller. 
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References to use for a study of Flux-Mapping techniques are: 

"Graphical Flux Analysis in Transformer Design", M. G. Leonard, Electro- 
Technology, Oct., 1961, pp 122-128. 

"Fundamentals of Electrical Design", A.D. Moore, McGraw-Hill Book. 

"Electromagnetic Devices", H.C. Roters, John-Wiley & Sons Book. 

Notes on Air-Gap and Interpolar Induction, F.W. Carter, IEE Proc. (British), 
Vol 29, 1900, p 925. 

I 

Air Gap Induction, F.W. Carter, El. World, Vol. 38, p 884. 

Sketches of Magnetic Fields in Iron, Th. Lehmann, Rev. Gen. El., Vol. 17, 
1926 

Mapping Magnetic and Electrostatic Fields, A.C. Moore, El. J., Vol. 23, 
1926, p 355. 

Fundamental Theory of Flux Plotting, A.R. Stevenson, G. E. Rev., Vol. 29, 
1926, p 797. 

Graphical Determination of Magnetic Fields, R.W. Wieseman, AIEE Trans., 
Vol. 46, 1927, p 430. 

Graphical Determination of Magnetic Fields, E.E. Johnson and C.H. Green, 
AIEE Trans., Vol. 46, 1927, p 136. 

Graphical Determination of Magnetic Fields, A.R. Stevenson and R.H. Park, 
AIEE Trans., Vol. 46, 1927, p 112. 

The Interpolar Fields of Saturated Circuits, T. Lehmann, AIEE Trans., Vol. 
46, 1927, p 1411. 

A Practical Application of Graphical Flux Mapping, J.F. Calvert, El. J., 

Vol. 24, 1927, p 543. 
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Graphical Flux Mapping, J.F. Calvert and A. M. Harrison, El. J., Vol. 25, 
1928. Theory and General Discussion, p 147; Fields of Non-Salient Pole 
Synchronous Machines, p 179; D-C Motors and Generators, p 399; D-C 
Motor, Salient Pole Synchronous Machine, Universal Motor, etc., p 510. 

Analytical Determination of Magnetic Fields, B.L. Robertson and I. A. Terry, 
AIEE Trans., Vol. 48, 1929, p 1242. 

Magnetic Fields in Machinery Windings, J.F.H. Douglas, AIEE Trans., Vol. 
54, 1935, p 959. 


SA-37 



THE MAXIMUM A RATIO FOR ROTATING COIL LUNDELL GENERATORS 
d 

In Brief 

The maximum practical A ratio for a rotating coil Lundell generator is .3. 

d 

Discussion 

The flux circuit of a rotating-coil Lundell generator can be represented by 
the mmf drops below: 


F gap F teeth F core F teeth F gap 



The rotor shaft completes the flux circuit, and the area of the shaft 
limits the useful stator length by limiting the amount of flux that can be 
carried inside the excitation coil. 

When a maximum length machine is made, about . 50% of the stator bore 
cross-section area is used up or made unavailable because of the thickness 
of the cantilevered poles at their bases. 

Another 20% of the stator bore cross-section area must be used for an 
excitation coil and only about 30% of the area can be utilized for shaft. 
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Assumptions are: 


Area Shaft = .3 HlL$L 

4 

Gap Density = 35 Kl/in 2 

Shaft Density at FL = 100 Kl/in 2 

Pole Embrace = .6 

Total Flux at FL = 2.5 x useful flux. (This assumes that the 
leakage flux is 1-1/2 times the useful flux). 

Then: 

Total Pole Flux = ^d^(35) C p Kilolines for all poles (this flux 
links the stator winding). 

Half of this flux passes through the shaft. 

Total Shaft Flux at FL = 

0 sh = 2. 5 = g..5(. 6) /(35) 

Ct 

= 26.2 ^d t 

0 sh also = 100 t 3^L (. 3) = 7.5 <d 2 = B gh (Area of Shaft) 

26.2 'ir'&jL = 7.5 ^d 2 
/= 2XT d = - 286d 

This can be rounded off to .3. 
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THE MAXIMUM A RATIO FOR TWO, INSIDE, STATIONARY COIL LUNDELL 
d ’ 

GENERATORS (OR BECKY-ROBINSON GENERATORS) 

Si Brief 

The maximum practical ratio for a Becky-Rob in son type of stationary- 
coil Lundell Generator is approximately: 



Discussion 

The flux circuit of a two-coil Lundell generator of the Becky-Robinson 
type is represented by the mmf drops shown below. 


F gap F teeth F core F teeth F gap 



F g Pole 

F skirt 
F gap 2 
F yoke 1 
F gap 3 


Note that the flux in one parallel branch circuit must pass through the 
rotor skirt, or flux collector ring, into the yoke and back through the 
shaft. 
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The combined shaft and flux collector ring areas must carry the equi¬ 
valent of the total flux that passes into and back out of the stator. The 
areas of the combined skirt or collector ring and shaft cannot greatly 
exceed 60% of the total cross-section area of the stator bore because 
of the room needed for the excitation coil. 

The following assumptions are now made: 

Shaft and skirt area r .6 1 

4 

Shaft and skirt density - 100 KZ/in^ 

Air-Gap density = 35 K//in^ 

Total flux in shaft - 2. 5 times useful flux 

The combined areas of the rotor skirt plus shaft should be equal to 
the combined shaft areas at the two ends or skirt area = shaft area. 

2. 5 x shaft flux - . 6 <f lf^ (100) 

4 

= 15 ¥ d 2 

Useful flux = _( 15 ) ¥d 2 = 6 ¥d 2 
2. 5 

Also, Useful flux = Tjl /(35) Cp K^in shaft 

2 
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6 tff' d 2 = 35^d^ (Cp = .65) 
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THE DEVELOPMENT OF EQUATIONS DESCRIBING THE WEIGHTS 
OF ELECTROMAGNETIC PARTS FOR THREE GENERATOR TYPES 


A weight breakdown is made for the electromagnetic parts of 1) two, 
inside, stationary-coil Lundell generators (Becky-Robinson Gene¬ 
rators), 2) two, outside-coil Lundell generators (Rice Generators) 
and 3) Homopolar Inductors. 

hi order to have simple but realistic relationships between the three 
generator types mentioned above, some simplifying assumptions have 
been made. 

1. In all cases, the tooth width is equal to the slot width. 

2. A pole embrace of . 65 is used and the pole constant is 
therefore approximately . 65. 

3. The gap density Bg is 50 Kl/w? for one set of calculations, 
and 35 K^/in^ for another set of calculations. 

4. To keep the machine reactances of the 4 pole, 6 pole, 8 pole 
generators approximately constant for comparison purposes, 
we have assumed that the gap density is constant and that 
ampere wires divided by poles times air gap length is also 
constant. 
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—- = constant. 

P g e 


For 4-pole a-c generators A = 600 
For 6-pole a-c generators A = 700 
For 8-pole a-c generators A = 800 

As a matter of academic interest, the gaps are: 

8-pole generator gap = g 
6-pole generator gap = 1.16 g 
4-pole generator gap = 1. 5 g 

In the programs note that all the reactance formulae are 
multiplied by a common factor designated as X. 

A Ka 

X - -—- = AK. When Bg is constant, then 

2 C x B g (10) 

\ d A 

x ad = X Cj C m /\ a = AK’ p — = constant and = constant. 

The constant-reactance approximation as stated is inaccurate 
by the variation in d when the number of poles is changed. 

This variation in diameter (4^ 10%) can be compensated for by 
further changing the length of the air-gap to keep the reactance 
constant. The assumption is realistic enough to allow useful 
weight comparisons between the three types of generators. 
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SLOT DEPTH (h s ) 

If a stator loading of 800 ampere wires/inch of stator is assumed and the 
slot width is equal to the tooth width there are 800 amperes in 1/2 inch of 
stator bore periphery. For 8,000 ampere/in 2 current density in the stator 
winding and a slot fill factor of 2/3, ^ = . 150 in 2 of slot is needed. 

The slot must be about -■ = . 3" deep. For a machine with a stator bore 

diameter of 6”, a . 3” deep slot is . 05 d and this slot depth estimate has 
been used in all calculations. 


CORE DEPTH (h c ) 

The useful flux per pole in the stator is 0 P = — '^g C P kJL^ and the area 
of the bore cross-section is h c £. If the core density is 80 K 1 /in 2 then, 


2 A c 


= 80 and 


< 7 T d l B gC p .'Vf'(Bg) .65d 


'S' • 

2 (80) P 


„ .01275 (B g )d 

p 


For 50 K#./in 2 gap density h c = . 638-p- 
For 35 K£/in 2 gap density h c = . 446-^ 
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OUTSIDE DIAMETER OF STATOR (D) 


The outside diameters of the stators are d + 2h s + 21^ 
h s = .05 d in all cases 

h,. = .638 A- for 50 k£ /in 2 
= .446 A- for 35 yJL /in 2 

D = (1.1) d + 1 - p 7 -- for 50 YJL /in 2 gap density 

= (1.1+.319) d = 1.419 d for 4P 
= (1.1+ .212) d = 1.312 d for 6P 
- (1.1 + .16) d = 1.26 d for 8P 

D = 1.1 d + . 892 for 35 /in 2 gap density 

= (1.1 * .223) d = 1.32 d for 4P 

= (1.1+. 147) d ■ 1.25 d for 6P 

= (1.1 + .111) d - 1.21 d for 8P 

YOKE THICKNESS HOMOPOLAR INDUCTOR 

The leakage flux of the homopolar inductor is assumed to be all of the flux 
that is not useful in generating voltage and is assumed to be . 4 the useful 
flux. This is a low estimate that can only be realized with careful design 
and is an estimate that makes the homopolar inductor compare favorably 
with the other machines. 
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any diameter of the end-bell the weight of the end bell must be 

EBwt. = /T(D-t- ty) ty m . 283 -j- flux shoe weight estimated at 
.125 (END-BELL WEIGHT) = EBwt. = 77 (1.125). 283 j~D + tyj t y lbs. 


ROTORS OF OUTSIDE-COIL LUNDELLS 

The length of the rotors for outside-coil Lundells must be about to 
accommodate the flux-collector rings at the rotor ends. When the 
environment and speed permits the rotor can be made approximately 50% 
solid. 

Wt. of rotor = d (. 5) (. 283) 

= . 166 d 3 lbs. 


ROTOR AND FIE UP-COILS FOR TWO, INSIDE-COIL LUNDELLS 
(BECKY ROBINSON) 


To accommodate the field coils and yokes in the rotor ends, the rotor of 
the two, inside, stationary-coil Lundell must be approximately 2 JL long. 

Since £_ - ^ , the rotor length is ^ — d. The rotor is considered 80% 
solid to account for the coils and yokes. 


Rotor wt. 



d (.283) 


. 333 d 3 lbs. 
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GENERATOR STATOR AMPERE LOADING 



If: 


1) The slots occupy 50% of the space in the stator bore area. 

2) The no load flux density in the iron is 100 K t/in . 

Then the air gap density will be 50 Kl^/in 2 at no load. 

In practice, the air gap density may vary from 20 K^/in 2 for 
a permanent-magnet generator to 65 KtL/in 2 for a low-reactance 
wound-pole generator. Solid-pole generators usually operate at 
moderate gap densities because of pole-face losses , 

The ampere-wires per inch of bore periphery A represents 
the current loading of the generator and the output power. If the 
value A is doubled and all else held constant, the power output 
doubles. In a 30 KVA, 8-pole MIL-G-6099A generator for air¬ 
craft, the value A will be about 700 to 800. The tooth width will 
about equal the slot width, the gap density is approximately 
50 K£/in 2 and the tooth no load density is about 100 K^/in^. 


t 

\ 

\ 

\ 


i 
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Yoke area = *7 J (D-j- t y ) ty 


- a fj‘ (D) t y approximately 

IT ®g^p 

Yoke flux -2- Kilolines of useful flux 


Yoke flux = 80 x yoke area 

I 


d 

2 


Yoke area = 


Td 2 B g C p 
4(70) 


^ Dt, = ' ^ d2 % c p 
/l T 280 


d 2 B g (. 65) 


D(280) 

_ d (50). 65 0817 d 

“ 1.42 (280) " * 00I7d < 


4 


0 


For 4 poles and Bg = 50 


/ 

Yoke weight = 7f( D + ty) t y d- x 3 (.283) 


r 


= 77*(1. 42 +. 082) .082 ^-x 3 (.283) 

Z i 


= . 164 d^ lbs. for Bg = 50 and P = 4. 


ROTOR, HOMOPOLAR INDUCTOR 

The rotor of the homopolar inductor is assumed to be 60% solid on the 
ends and 100% solid in the center section between the stators. The diameter 
of the center section or shaft area is .8 d. The rotor length is^per section 
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or per section. 


Rotor weight = 


| (.60) x2 + 


' K ( d x - 8 ) 2 1 .283 lbs. 


= d 3 [l. 2 + . 64 ] = . 204 d 3 lbs. 


YOKE AND END-BELLS OF OUTSIDE-COIL LUNDELL GENERATORS 

The useful flux in the yoke is C P Kilolines and the leakage flux 

is JL as great as the useful flux. The flux density in the yoke, due to the 
3 

useful flux alone is 60 K£/in 3 . 


Area yoke - ‘TTD ty approx. 2 

-r . 65(Bo-) 

t„ = thickness of yoke = -x-2—-— t — 2 . 

y Jf D (60) 


Yoke length = z£ = 3| 


Yoke weight = (D +- t y ) t y -= d (. 283) lbs. 


The end bells of the two, outside coil Lundell generators must carry the 
yoke flux and must have approximately constant cross-section area from 
the outer yoke to the auxiliary air gap. IfTT^EB (thickness) = constant at 
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The limit of ampere loading on the stator is reached when ad¬ 
ditional rotor ampere turns, which would be necessary to overcome 
more armature reaction and armature voltage drop, only serve to 
force more flux across the rotor leakage paths. 



A small (1 KW) conventional, 8000 rpm, 400 cycle permanent 
magnet generator with block-type salient-pole magnets will produce 
about 600 ampere wires/inch at maximum output. The regulation 
is about 40% at such a load and in order to maintain a useable 
voltage level about 250 ampere wires/inch is the usual limit. 

A 10 KVA, 4-pole, wound rotor, salient-pole MIL-G-6099A 
generator may have a rated load output of about 500 to 600 ampere- 
wires per inch. 


A 30 KVA, 2-pole, 400 cps, non-salient pole generator may 
only be able to carry 400 ampere-wires/inch if it is to meet 
MIL-G-6099A performance requirements. 

There is no direct and clear-cut way to arrive at the maximum 
stator loading. For different applications, the load limit may be 
entirely different. In general though, for any of the electromagnetic 
generators we are likely to be interested in the loading A will 
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probably lie between 600 and 900. For very small machines and 
PM generators, this value may be as low as 200 while for large 
synchronous condensers on utility systems (50,000 KVA) the value 
may be 1200. 

The ampere loading of the stator is related to the reactances 
of the machine. 

If wire is wound through an iron or steel machine, it acquires 
inductance. 


N (Ni) (P 
i 

Where permeance 

This indicates that since the output current is fixed in any specific 
application, the inductance of the stator winding will vary directly as 
the square of the conductors/inch used or as A . 
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GROUPING OF FRACTIONAL SLOT WINDINGS 


When the stator is comprised of a winding having fractional slots per pole, 
the grouping of the coils can be determined by the following method. 

Express the ratio of the number of slots to the number of poles as an improper 
fraction reduced to its lowest terms. The denominator will then represent the 
number of poles in a repeatable section and the numerator will then be the 
number of slots in a repeatable section. If the number of slots in a repeat- 
able section is not divisible by the number of phases a balanced polyphase 
winding cannot be obtained. 

The maximum number of parallel paths in the winding is found by dividing the 
number of poles of the machine by the number of poles in a repeatable section. 
This gives the maximum possible number of parallels. 

To determine the grouping lay out a table having as many horizontal divisions 
as there are slots per repeatable section, and as many vertical divisions as 
there are poles per repeatable section. Divide the horizontal divisions into 
m number of phases and using a pitch of y = poles per repeatable section lay 
out the winding. 

As an example consider a 3 phase, 20 pole, 84 slot machine, then 84/20 = 

21/5 = (slots per repeatable section)/(poles per repeatable section). The 
maximum number of parallels possible = 20/5 = 4, and the available number 
of parallels will then be 4 and 2. Lay out the table as follows with y = 5 and 
throw = 1-6. . 


1 


1 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

n 

X 

■ 




X 





X 





X 





X 

H 





X 





X 





X 





X 


E3 




X 



I-■ 


X 





X 





X 



a 



X 





X 





X 





X 




Cl 

— 

IT 





X 

__ 




X 





X 






Phase a Phase c Phase b 


sa-53 



The grouping will thus be 21211-21211 = 21211 = and repeat 3 x 4 = 12 times. 
The coils will then be placed as follows: 


Slot No. 

1 

2 

— 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

Phase of 

+ 

+ 

- 

+ 

+ 

- 

+ 

- 

- 

+ 

- 

- 

+ 

- 

+ 

+ 

- 

+ 

+ 

- 

+ 

Coil 

a 

a 

c 

b 

b 

a 

c 

b 

b 

a 

c 

c 

b 

a 

c 

c 

b 

a 

a 

c 

b 

Grouping 

< 

4 

i 

1 

2 

1 

1 

2 

i 

2 

1 

i 

2 

1 

2 

1 

1 


Another way of accomplishing the same result as that given above is as follows: 


Since there are 5 poles per repeatable section and 21 slots per repeatable sec¬ 
tion there must be 7 slots per phase in each repeatable section with these 7 
slots occupying positions over 5 poles. The coils in a section may be arranged 
in any order but the best arrangement in general is that which gives the largest 
distribution factor. Usually the most symmetrical arrangement of the coils is 
best, or in this case 12121. Laying out the winding of the previous example for 
the 12121 arrangement gives: 


Slot No. 

1 

2 

3 

4 

7 

6 

7 

8 

7 

10 

ii 

12 

13 

14 

15 

. 

16 

17 

18 

19 :20 

21 

Phase of 

+ 

- 

- 

+ 

- 


+ 

- 

+ 

+ 

- 

+ 

+ 

- 

+ 

- 

- 

f 

1 

+ 

Coil 

a 

c 

c 

b 

a 

! a 

c 

b 

a 

a 

e 

b 

b 

a 

c 

b 

b 

a 

C S c 

! b 1 

Grouping 

1 

i 

2 

1 


2 

1 

'll 

< 

2 

i 

i 2 

i 

i 

i i ! 

A 

1 

—.. 

2 

i i 

i. 

2 j 1 


which is the same as that arrived at previously with tfye exception that it is 
displaced several slots. 


THE SLOT STAR 

The two above methods of grouping can be verified by the slot star method. 

In general let the number of slots per phase per pole q be represented as 
q = N/B = a + b/B. where a is an integer. Then the winding repeats itself 
after each B poles and the number of recurrent groups is equal to P/B. Each 
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phase has N slots in B poles. Further, each phase has B - b coil groups with 
a coils and b coil groups with a + 1 coils in B poles. If B/m equals an integer 
there cannot be a perfectly balanced fractional slot winding. 

Apply the above rules to the 3 phase, 20 pole, 84 slot example gives q = 84/ 

(20 x 3) = 7/5 = 1 + 2/5. Thus the winding repeats itself after each 5 poles the 
number of recurrent groups is 20/5 = 4. Each phase has 7 slots in 5 poles 
and each phase had 5-2 = 3 coil groups with 1 coil and 2 coil groups with 
1+1 = 2 coils in 5 poles. 

The slot star of a 2 pole, 3 phase integral slot winding with q equals two is 
shown in Figure 1. The angle between two adjacent slots is cc = 180°/mq = 
30°. 

Two adjacent vectors correspond to two adjacent slots and thus slots 1, 2, 7, 
and 8 belong to phase I, slots 3, 4, 9, and 10 belong to phase n and so on. Vec¬ 
tor 7 with which the second pole starts is shifted 180° with respect to vector 1, 
and the same applies for vectors 2 and 8, 3 and 9, etc. Therefore the bottom 
half of the slot star is the same as the top half except that their vectors are 
shifted by 180°. 

The four coils of phase 1 are shown in Figure 2 with the solid lines represent¬ 
ing the tops of the slots and the dotted lines representing the bottom. The con¬ 
nector C takes care that all the emfs add together and thus the slot star of this 
example is completely represented by only half of the circle. 

The above type of star will then apply to the integral slot windings in general. 

In the fractional slot winding, B poles make the recurrent group and just as 
for the integral slot winding the slot star of B poles is represented by half of 
a circle. 

Figure 3 shows the slot star of the 20 pole, 3 phase, 84 slot example. There 
are 5 poles in B with 7x3 = 21 slots per recurrent group, oc = 180/ (3 x 1.4) 

O 
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= 42 and 6/7°. Thus the angles which correspond to the slots are: 

Slot No. 1 2 3 4 5 6 

oc g - 0 42 6/7 85 5/7 128 4/7 171 3/7 214 2/7 = 34 2/7 etc. 

Since the largest distribution factor for the fundamental is obtained when the 
slots belonging to each group are closest together, the first seven slots in the 
slot star will be assigned to phase a. The second seven to phase c, and the 
last seven to phase b. Thus slots 1, 18, 14, 10, 6, 2, and 19 are phase a; 

15, 11, 7, 3, 20, 16, and 12 are phase c; and 8, 4, 21, 17, 13, 9, and 5 are 
phase b. 

The grouping then becomes 


Slot No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

Phase of 

+ 

+ 

- 

+ 

+ 

- 

+ 

- 

- 

+ 

- 

- 

+ 

- 

+ 

+ 

- 

+ 

+ 

- 

+ 

Coil 

a 

a 

c 

b 

b 

a 

c 

b 

b 

a 

c 

c 

b 

a 

c 

c 

b 

a 

a 

c 

b 

! 

Grouping 

2 

; 

1 

2 

i 

1 

2 

i 

2 

1 

i 

2 

1 

2 

1 

I 

1 


This is identical with the grouping of the first method and therefore verifies 
its accuracy. 

In the integral slot winding the beginnings of the phases are displaced by 120 
and 240 electrical degrees. Also in the fractional slot windings the distances 
between the beginnings of the phases can be made 120 and 240°. This will be 
the case when the beginnings are placed in slots 1, 1 + N, and 1 + 2N when B 
is an even number, and in slots 1, 1 + 2N, and 1 + (1 + m) N, when B is an odd 
number. This arrangement of the beginnings of the phases will place them far 
apart from each other mechanically while it is often desirable to have them 
near to each other. In order to place the beginnings of the phases near to each 
other the beginnings can be placed approximately 120 and 240 electrical de¬ 
grees apart and the windings still will be balanced. This is due to the fact that 
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in the fractional slot winding the emfs of the consecutive coil groups are not 
in phase and the sequence of the geometric addition of the single emfs is of no 
influence on the resultant phase emf. So in the example considered the be¬ 
ginnings of the phases can be placed in slot 1 for phase a, slot 4 for phase b, 
and slot 7 for phase c. The angles between the beginnings are then 128 and 
4/7°, and 257 and 1/7°. 


DISTRIBUTION FACTOR 

The voltages induced in the separate coils of a distributed winding are not in 
exact phase and their resultant is therefore less than would be produced in a 
concentrated winding having the same number of turns. The ratio of the volt¬ 
ages produced by distributed and concentrated windings having the same num¬ 
ber of turns is called the distribution factor. In the case of integral slots per 
phase per pole it can be derived as follows: 

0 = electrical angle per phase group 

6 - q° t s 


sin 


«s _ e/2 _ e 
~2 r _ 2r 


r = 


oc 


2 sin 


T~ 


0 qoc s 

sin g = sin —j— 


E 

2r 


qoc 


E = 2r sin 


qoc 


qoc. 


e sin- 


= 21 


>2 sin ^sy 
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Thus K d = 


i oc 
sin s 


. qo = 


~T 


, oc 

q sin s 


~ 2 ~ 


Since a displacement of oc g between the slots is noc g for the nth harmonic, 


qnoc 


sin 


s 


K. for the nth harmonic is K, 
dn dn 


noc 


q sm 
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FRACTIONAL SLOT DISTRIBUTION FACTOR 

Refer to the slot star shown in Figure 3 of the section titled "Grouping of 
Fractional Slot Windings" and it will be noted that to determine for a frac¬ 
tional slot winding it is necessary to distinguish between the angle between 
two slots oc g and the angle between two vectors °c m - This latter angle is the 
magnetic field angle between the slots of the recurrent group and this angle 
determines the phase difference between the vectors. The magnetic field 
angle is oc m = 180°/Nm where N is the same as in "Grouping of Fractional 
Slot Windings". 


It can be seen from the slot star that the fractional slot winding thus behaves 
like a winding with N slots per phase per pole shifted with respect to each 
other by the magnetic field angle oc . Therefore the distribution factor is 


Noc 


sin 


m 


K. = 


oc 


N sin 


m 

T 


The general effect of distributing a winding is to smooth out the wave form by 
diminishing the amplitude of the harmonics with respect to the fundamental. 
The distribution of the armature copper loss is also improved. 
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The distribution factor of the three phase winding is greater than that of the 
two phase winding and for this reason the three phase winding is used where 
there is a free choice of the number of phases. The distribution factor of the 
single phase winding is much smaller than that of either 3 or 2 phase windings 
because the winding is distributed over a larger arc. 


In general only 2/3 of the slots per pole are used in winding single phase ma¬ 
chines. The reason for this can be best shown by an example. Let the number 

of slots per pole equal 9, and if all slots were wound then cc = 180°/9 = 20°, 

s 

and 

sin (9 x 1 

K d = 9 x sin 10° = TT5M = * 640 

If only 6 slots are wound 


sin (6 x 10) 
d 6 x sin 10 


.832 


9 x 640 

Thus, for 9 slots the number of effective turns is only = 1.15 times 

the number of effective turns obtained when using 6 slots, and therefore by 
using 50% more copper with its additional 50% more loss, only 15% more volt¬ 
age has been obtained. 


SKEW FACTOR 

It can be noted in the table of distribution factors that some harmonics have 
the same distribution factor as the fundamental. These harmonics are called 
the slot harmonics and their orders are 


n 



where K is an integer 


The slot harmonics which correspond to K = 1 (slot harmonics of the first 
order) are among the most troublesome harmonics in AC machines. Their 
influence, as well as the influence of other harmonics of higher order, can 
be reduced by skewing and for this reason stator or rotor slots are sometimes 


SA-61 



skewed. The skewing also reduces the flux variation in the fringing of the flux 
at the pole tips due to the slots entering and leaving the polar region. Such a 
flux variation oftentimes contributes to noise. 

Skewing has the same effect as the distribution of a winding over a larger zone 
because it reduces the interlinkages between the field and stator windings. 

This distribution factor due to skewing is called the skew factor. 

For slots that have a large number of slots per coil group the path of the vec¬ 
tors being added approaches the arc of a circle, and when this happens the 
distribution factor can be expressed as the ratio of the chord AE to the arc 
AE and 


qoc g 
sm -j- 


1/2 chord AE 
F 


qoc 

s 

chord AE = 2R sin —^— 
arc AE = Rqoc 

s 

qoc qoc 

2R sin x - sin — 
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d Rqoc s qoc g 
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If Z is the arc which the coil group occupies per pole and t is the pole pitch 
then 
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Let t g ^ be the slot skew in inches for a length equal to the core length of the 
machine and the skew factor K , is 
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The influence of skewing is nigligible for the fundamental, small for the har¬ 
monics of low order, and very considerable for the harmonics of higher order. 
A skew which is equal to one stator slot pitch makes the influence of the dan¬ 
gerous slot harmonics almost negligible. 


PITCH FACTOR 

Fractional pitch windings decrease the length of end connections, reduce slot 
reactance, and provide a means for improving the wave form. They can be 
used to eliminate any one harmonic from the voltage wave as well as to reduce 
other harmonics. However, they require a few more turns or a greater flux 
for the same voltage than a fuH pitch winding. 

Since the two sides of a coil of a fractional pitch winding do not lie under the 
centers of adjacent poles at the same instant the voltages induced in them are 
not in phase when considered around the coil. The voltage produced is there¬ 
fore less than that which would be produced in a full pitch winding. The voltage 
generated in any single turn is the vector difference of the voltages generated 
in the two inductors which form the active sides of the turn. If the throw in 
slots is y, then (y/mq) x 180° will be the angle of phase difference between 

the turns and K is then derived as follows: 

P 

AB = BC = E; since AB = BC, angle C = angle A, and the bisector of (y/mq) x 

o T 

180 will be perpendicular to AC and will bisect AC into two equal parts E ; 

thus 

Kp = 2E /2E = E /E = sin £(y/mq) x (180°/2)J = sin x 90°j 
Since the displacement for any harmonic such as the nth is n times the phase 
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displacement of the fundamental, for any harmonic n is 


K 


= sin(J* x 90°) 

\mq ) 


Any harmonic can be eliminated by choosing a pitch that makes the pitch fac¬ 
tor zero for that harmonic. Thus to eliminate the nth harmonic it is only neces¬ 
sary to select a pitch such that n(y/mq) x 90° equals 180°, 360°, 540°, etc. 

(any multiple of 180°). Eliminating any one harmonic also reduces other har¬ 
monics and the fundamental by different amounts. A pitch of 5/6 will give 
minimum fifth and seventh harmonics and should theoretically give a minimum 
additional rotor surface loss under load conditions. 
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Reactances 


Machine reactances are sometimes expressed in ohms but are more 
often given in per cent or in per unit. The per unit system is discussed 
in Appendix A of Electric Machinery by Fitzgerald and Kingsley, pages 
543-546. In our treatment of generators, the per unit valuesof the 
various reactances are obtained by dividing the ohmic value of the 
reactance by the base ohmic value, which is the phase voltage divided 
by the phase current or Eph/Iph. 

When the generator values are expressed in P. U. (or per cent which 
is 100 x P. U.), all generators can be compared by the same standards. 
A 50,000 KW, 60 cycle generator and a 1 KW, 1000 cycle generator 
look pretty much alike as far as their reactances and per unit values 
are concerned. 

Per unit values can be easily used in system studies and the values 
are easy to convert to a new base when used in a system involving 
several machines. 

Nearly all of the formulae in this study are arranged to give per cent 
values. To convert to ohmic values, multiply the per cent value by 

Eph 1 
Iph X 100 
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Any rotating electrical machine has only one zero-sequence reactance 
and only one negative-sequence reactance but it has several positive- 
sequence reactances. These positive sequence reactances depend upon 
the angular position of the rotor and upon whether the positive-sequence 
currents are steady or are suddenly applied. 
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SYNCHRONOUS REACTANCE 



X 


d 


X 




When a generator operates on a steady state symmetrical short circuit condi¬ 
tion Its terminal voltage Is zero and its saturation is negligible. Since there 
U no terminal voltage the net stator linkage must be zero, and thus the stator 
linkage due to its current acting alone must be exactly equal and opposite in 
direction to the stator linkage due to the field current acting alone. If the field 
current were acting alone with the stator open-circuited a certain terminal 
voltage would exist, and since the armature linkage has the same value as the 
field linkage, the voltage which must be appUed to produce the stator current 
must be exactly the same as the terminal voltage induced under open circuit 
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Therefore the unsaturated synchronous impedance is the ratio of the phase vol¬ 
tage on open circuit resulting from a certain field current to the steady state 
short circuit stator current resulting from the same field current. Further, 
since the value of the effective resistance is usually very small in comparison 
with the reactance it can be neglected, and the impedance and reactance are 
then equal. Thus 



E (open circuit) 
I (short qircuit) 


ohms 


and since the voltage drop due to the stator current is due to the fictituous re¬ 
actance X a( j and the actual reactance X-t 



xt 


Direct Axis Synchronous Reactance 
Definition AIEE TEST CODE 503, Para. 1.820 

"The direct-axis synchronous reactance is the ratio of the fundamental 
component of reactive armature voltage due to the fundamental direct- 
axis component of armature current, to this component of current under 
steady-state conditions and at rated frequency. 

Quadrature Axis Synchronous Reactance Xq 
Definition AIEE TEST CODE 583, Para. 1.830 

"The quadrature-axis synchronous reactance is the ratio of the fundamental 
component of reactive armature voltage, due to the fundamental quadrature- 
axis component of armature current, to this component of current under 
steady-state conditions and at rated frequency. 
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FLUX CIRCUIT FOR STEADY STATE CONDITIONS 



The above schematics show the circuits for steady-state fluxes in the direct 
axis and the quadrature axis. The two conditions shown would be obtained 
when testing for the synchronous reactances, X d , and Xq. Under steady 
state loading conditions, the stator flux wave is stationary with respect 
to the rotor poles, since it is traveling at synchronous speed. Therefore, 
the damper circuit and the field windings have no effect on the flux. 
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REACTANCE OF ARMATURE REACTION 


Since the shape of the resultant MMF wave is independent of the direction of 
armature reaction the effect produced by the MMF's of the field and stator 
windings can be found by treating the two forces as if they each acted alone, 
and thus the forces may be replaced by the voltages they would cause if acting 
separately. If this substitution is made the voltage due to armature reaction 
may be considered as being a voltage drop due to a fictituous reactance X ad , 
and this is called the reactance of armature reaction. It is not an actual re¬ 
actance but under steady operating conditions may be considered as such in 
order to simplify the methods of calculation. It is in phase with the voltage 
drop due to leakage reactance I^X^, and the sum of these two reactances is 
the synchronous reactance X^. 


If the effect of saturation is neglected the saturation curve becomes a straight 
line, and then any change in flux with its corresponding change in voltage, 
produced by any change in MMF is proportional to the change in MMF. Thus, 
if an unsaturated condition is assumed, the reactance of armature reaction in 
per unit value can be seen to be the ratio of the MMF of armature reaction to 
the MMF required by the field to force the flux across the air gap, or 

F 

X ad = (percent) 
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FLUX CIRCUIT FOR TRANSIENT CONDITIONS 



The two schematics above show the flux circuits for transient conditions 
in the salient-pole, wound-pole synchronous generator. These conditions 
obtain when tests are made for x' d and X' q or when the stator armature 
reaction flux wave is changing with respect to the poles. 


X^ Transient Reactance 


When armature current is abruptly applied to a synchronous 
generator it may cause the sudden appearance of an opposing mmf 
opposite each field pole tending to establish flux through the pole 
core. 

In a laminated salient pole rotor with wound poles, the transient 
flux is opposed by currents in the field winding. These initial 
currents decay until only the internal impedance X^ opposes the 
current flow. 

The permeance of the leakage path for the transient flux 
determines the armature transient current level. 
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TRANSIENT REACTANCE IN THE STATIONARY- 


COIL BRUSHLESS GENERATORS 


In the solid-rotor or stationary-coil, a-c generators, 
the field-coil is separated from the main air-gap by 
two auxiliary air-gaps. The transient flux is prevented 
from passing through the field coil and is shunted 
through the field leakage paths. These field leakage 
paths are large and the resulting transient reactance 
is large. 

The following sketches show the path of the transient 
flux in a two, inside, stationary-coil, Lundell generator, 
or Becky-Robinson generator. This sketch is generally 
representative of the condition existing in any of the 
stationary-coil generators. 

During the first cycle or two of the transient, the flux 
will be opposed by eddy currents in the rotor surface. 
These eddy currents reduce the transient reactance 
by about 15% during these first one or two cycles. 
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UNSATURATED TRANSIENT REACTANCE — (X* ) 

du' 


When current is suddenly applied to a generator stator (a step load) there is 
no flux linkage with the field winding during the first instant, and during that 

first instant all flux due to the armature mmf is forced through the field coil 
leakage paths. 



TRANSIENT FLUX IN THE NORTH POLE OF 
A BECKY-ROBINSON LUNDELL, A-C GENERATOR 
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TRANSIENT FLUX IN THE SOUTH POLE OF A BECKY— 
ROBINSON, LUNDELL A-C GENERATOR 



\0 j3 O O/Of 0\0,0,0 O 


\ \ + 
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FLUX CIRCUIT FOR SUBTRANSIENT CONDITIONS 



Xd Subtransient Reactance 


If armature current is abruptly applied to a synchronous 
generator which has a damper circuit on the rotor poles, or 
which generator has a solid pole face, currents will be set up 
in the damper circuit or pole surface. These currents oppose 
the change of flux linkages in the rotor circuit and the permeance 
of the path of the transient leakage flux determines the level of 
the transient current in the armature. 

The subtransient current decays quickly leaving the transient 
current, which decays less quickly until at last only the steady 
state current limited by X^ remains. 

The subtransient reactance of solid-rotor machines is about 
15% lower than the transient reactance of the same machine. 
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X 2 Negative-Sequence Reactance 


Negative sequence currents applied to the stator produce a mmf 
rating backwards at synchronous speed. With respect to the rotor 
which is rotating forward at synchronous speed, the negative se¬ 
quence mmf and flux waves are rotating backwards at twice 
synchronous speed. 

Currents of twice rated frequency are induced in all rotor 
circuits keeping the flux linkages of those circuits at almost con¬ 
stant zero value. The flux due to the negative sequence armature 
current is forced into paths of low permeance which do not link 
any rotor circuits. These paths are the same as for the sub¬ 
transient flux. 

Since the negative-sequence flux wave alternately meets per¬ 
meances of the two rotor axes corresponding to subtransient re¬ 
actance X^' and X^' it is usually calculated as 


V ,f V' ’ 

X,* - *q 

1 2 

the arithmetical mean. 


Xq Zero-Sequence Reactance 


Zero-sequence currents in the armature produce’slot-leakage, 
end-leakage, and differential-leakage fluxes, but not the same 
leakages as are produced by positive-sequence current. The 
leakage flux produced is small and Xq is the smallest of all the 
reactances; lower than X^. 

The value of Xq varies through a wide range and depends 
upon the coil pitch. It is least for 2/3 pitch. 
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LEAKAGE REACTANCE 


In addition to the demagnetizing action of armature reaction there also exists 
a voltage drop due to the leakage fluxes. The lines of leakage flux go partly 
across the slot from one wall to the other (slot leakage), partly from tooth top 
to tooth top, (tooth tip), partly from phase belt to phase belt (zig zag), and 
partly in the end windings (end leakage). Each of the leakage fluxes is directly 
proportional to the current which produces it because its reluctance is princi¬ 
pally in air. 


Consider first the slot leakage flux. This flux includes all the flux which links 
the portion of the conductors that are embedded in the iron, but which does 
not enter the air gap. Assume an integral slot, full pitched winding, and as¬ 
sume that the flux passes directly across the slot. The effect of the notches 
at the wedge will be neglected and the current density in the conductors will 
be assumed constant. Also, for simplification, the distance between the top 
and bottom coil will be assumed negligible. 


The flux produced by the element in the 
bottom of the slot will surround this ele¬ 
ment. It passes across the slot above 
the element and returns in the iron be¬ 
low it. The flux will thus encircle part 
of the coil plus all of the slot above the 
coil. For the part above the coil the flux 
will proportion itself equally over the 
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3.19n a 
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distance h, /0 + h 0 , and since the inductance L equals-the L per slot 

11 4-010 

per inch of core length for this part is: 


3.19n g (hj /2 + h 2 ) 
L = -T- 


4b 10 
s 


For the inductance of the part of the coil itself consider the element dx enclosing 
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part of the conductor. The L per slot per inch of core length for this part will 

,2 


O 4Q N 

be equal to K ~ and since N = 

i rfi X 


10 


‘l/2 


n g b 
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and the total inductance of the bottom coil side is 


l b = 


3 ’ I9 H ,h , V| 3 * 19n s | 1/2 h 2 


In the same manner the inductance of the top coil side will be that due to its 
own current plus that due to the flux in the part h 2 above the coil and thus: 


, 3 - 19 ”s/S/2 \ 

T " 3 j 



Since a full pitched winding has been assumed, the currents of both coil sides 
will be in phase in all slots and the mutual inductance in the top coil due to 
the current in the bottom coil will be equal to that in the bottom coil due to 
the current in the top. Thus L< M in the top coil due to i in the bottom coil will 
be 



3.19n 2 
_s 



current in the top coil enclosed by the flux from 


the bottom. 
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At the distance x the current enclosed is 


3.19n 


l l/2 
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For the portion above the slot 
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as in the previous cases. 


Therefore, since total inductance is L T -h the total inductance 

per slot per inch of core length is 

r% 
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Since the portion of the slot between the coils has been neglected it can now be 
included by making hj ^ in the above derivation equal to 

h. 


T 


and the formula for L will then be: 


3. 19n / h h 


The tooth tip and zigzag leakage components have been derived by Kilgore with 
flux plotting methods and these have been determined to be 
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Adding these to the slot leakage portion gives 


3.19n 


.35b, 


L = 


10 


W 


3b 


s 


«V 


The above derivation has been based upon a full pitched integral slot winding. 
When the winding is chorded, however, as is usually the case in polyphase 
machines, the currents in the two coil sides in a certain number of slots will 
be out of phase and the coefficient of self inductance will be smaller than for 
the full pitch case. For the slots in which both the bottom and top layer belong 
to the same phase the conditions are the same as for full pitch. Both coil 
sides carry currents of the same phase and the phase angle between the cur¬ 
rents is zero. In the slots in which the bottom layer and top layer belong to 
different phases, the phase angle between the currents is not zero and the 
mutual inductance between the two layers will be reduced. The reduction fac¬ 
tor is equal to cos ft . 

Two reduction factors have to be used: K for the slot part in which the 

xco 

conductors lie, and K xt for the part above the conductors. Thus L becomes 

a w, r (m (M b t -»»t' 

io*~ “°l 3 v * \*>J i( v v 

The factors K xcQ and will vary according to the percent pitch and when 
they are derived they are found to be independent of the number of slots per 
phase per pole. Hence, since they depend solely upon the percent pitch, they 
have been combined into a single factor K which is reasonably accurate for 
most machines. It is determined as follows: 


=ifJy + 
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for three phase machines 


K = - 2 - 
x mq 


for two phase machines 
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Therefore the leakage inductance per slot finally becomes 


. 319n s „ l b 2 h l b t - 35b t 

L = ~W * » S + 3b s + 16t s g + ‘s 


and the per unit reactance per phase will be 
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Many attempts have been made to determine accurate formulas for the end 
winding leakage reactance, but the one that appears to be most promising for 
small machines is the one proposed by E. C. Barnes in AIEE Volume 70-1951. 
Expressed in per cent notation this formula will be: 
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Where ~Ti~ is P r °P° rti °nal to the pole pitch and is taken from Graph No. 1 
and K e is proportional to the ratio of the calculated L fi to the L £ of Graph 
No. 1. 
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POTIER REACTANCE 


The terminal voltage of an alternator under load differs from its open circuit 
voltage at the same field excitation. This difference is due to a voltage drop 
through the armature caused by leakage reactance, armature effective resis¬ 
tance, and armature reaction. The relative importance of the three factors 
depends upon the power factor of the load. With a reactive load at zero power 
factor the decrease in the terminal voltage is due almost entirely to the arma¬ 
ture reaction and the armature leakage reactance. Under this condition the 
effective resistance drop is in quadrature with the terminal voltage, and since 
it is small in magnitude, it has little influence on the change in the terminal 
voltage caused by a change in load. Likewise, the resultant field F NL is al¬ 
most exactly equal to the algebraic difference between F NL and F DM , and the 
terminal voltage E is nearly equal to the algebraic difference between E and 

B 

IphX. Under these conditions the armature reaction subtracts almost directly 
from the impressed field and the armature leakage reactance drop subtracts 
almost directly from the generated voltage. It follows from this that if an open 
circuit characteristic OB, and a curve CD are plotted as in Figure 14, show¬ 
ing the variation in the terminal voltage with excitation for the condition of 
constant stator current at a reactive power factor of zero, the two curves are 
so related that any two points, as E and F, which correspond to the same de¬ 
gree of saturation and consequently to the same generated voltage, are dis¬ 
placed from each other horizontally by an amount equal to the armature reac¬ 
tion and vertically by an amount equal to the leakage reactance drop. 

GF represents the armature reaction in equivalent field amperes and GE repre¬ 
sents the leakage reactance drop in volts. The armature leakage reactance per 

pip _. 

phase for a Y connected generator is thus X = / 73 I and this reactance is 

called the Potier reactance. The triangle EJF is known as the Potier Triangle. 
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If the leakage reactance is constant and the increase in field current necessary 
to balance a given number of ampere turns of armature reaction is independ¬ 
ent of the saturation of the magnetic circuit, Potier triangles drawn between 
the open circuit saturation curve and the zero power factor curve at points 
corresponding to different degrees of saturation would be identical. Under 
these conditions the zero power factor curve would have the same shape as 
the open circuit saturation curve but would be displaced from the open circuit 
curve by a distance equal to the length of the hypotenuse of the Potier triangle. 
Actually, the field pole leakage increases somewhat with an increase in satu¬ 
ration, and for this reason the number of field ampere turns which are neces¬ 
sary to balance a fixed number of ampere turns of armature reaction is not 
quite constant. In spite of this change in field pole leakage, the curves have 
nearly enough the same shape for practical purposes and are assumed to have 
the same shape when determining generator performance. 


In order to make use of the Potier method to construct a zero power factor 

curve it is necessary to locate two points on the Potier triangle. Figure 14 

shows how the method is applied. A triangle OCL is constructed with OC equal 

to the amperes corresponding to short circuit ampere turns (F = X , F ) 

v sc d g' 

The altitude LT is made equal to IX^ where Potier’s reactance is calculated 
by a partly empirical method and 


Xp = xA 


F + F 
TR CR 


^TR + F CR + F t + F C_ 


FS 


Ars 


FS 


2d 

Pg. 


A 


rs 


As the point L of the triangle is moved along the no load saturation curve, 
point C traces the zero power factor saturation curve. 
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FIGURE. 14 
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TRANSIENT AND SUBTRANSIENT REACTANCES AND TIME CONSTANTS 


Figure 9 shows the short circuit current in one phase of a three phase alterna¬ 
tor that has had all three phases shorted when operating at no load. The line ef 
drawn midway between the two sides of the envelope is called the direct cur¬ 
rent component and has different magnitudes in the three phases. The envelope 
of the alternating component is redrawn in Figure 10 with its axis horizontal, 
Le -’ with tbe direct current component eliminated. The alternating components 
are the same in all three phases. When an alternating voltage is short cir¬ 
cuited through an inductance and a resistance in series, the short circuit cur¬ 
rent will consist of two components. These components are the unidirectional 
or DC component which decreases logarithmically to zero, and an alternating 
component which is fixed by the voltage, the resistance, and the reactance of 
the circuit, and is constant when the resistance and reactance are constant. 

In the short circuited generator the reactance is not constant until the condi¬ 
tion equivalent to synchronous reactance is reached. When a sudden change 
occurs in the stator current the armature reactance is no longer constant, 
and under this condition voltages are induced by the changing armature reac¬ 
tion in the field winding, and in any other closed windings on the field structure. 
As the stator current builds up the voltages of the field structure will build up 
simultaneously, tending to maintain constant the total number of ampere turns 

acting on the magnetic circuit. These voltages cause transient currents in 
these parts. 

In a polyphase generator the alternating components of the short circuit stator 
currents produce armature reaction which is fixed in direction with respect 
to the poles but decreases from an initial value to a final value fixed by the 
steady state short circuit currents. To balance this increase in armature re¬ 
action there must be an increase in the field current. This increase in field 
current is in the same direction as the initial field current since the armature 
reaction caused by lagging currents is demagnetizing. The change in field 
current decreases and becomes zero, when steady state conditions are reached. 



The D. C. components in the stator currents produce a resultant MMF 
which is fixed with respect to the stator but has fundamental frequency 
with respect to the field. To balance this the field current must contain 
an alternating component of fundamental frequency. This alternating 
component in the field current produces an MMF of fundamental fre¬ 
quency in the air gap which is fixed in direction with respect to the field 
poles. This MMF can be resolved into two oppositely rotating components, 
each of which rotates at synchronous speed with respect to the stator. 

The component which rotates in the direction of the rotor balances the 
the stator MMF caused by the D. C. components in the stator currents. 

The oppositely rotating component rotates at double synchronous speed 
with respect to the stator winding and must be balanced by second har¬ 
monic components in the transient stator currents. All except the 
fundamental alternating components decrease to zero when steady state 
conditions have been reached and the fundamental alternating components 
become the steady state short circuit currents. 

Consider the elementary generator shown in Figure 11. The stator wind¬ 
ing is open and the rotor winding is excited by a D. C. current of magni¬ 
tude If. At the time t = 0, when the axes of both windings are perpen¬ 
dicular to each other, the stator winding is suddenly short circuited. As 
stated previously, the total flux interlinked with each winding under this 
condition will remain constant. Thus, the total flux interlinked with the 
field winding at t = 0 will consist of two parts; one part going through 
the path of the main flux, and the other part going through the leakage 
path of the rotor. The total flux interlinked with the armature at t = 0 


is zero. 
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FIGURE 12 


FIGURE 13 






During the time t the rotor moves through an angle a = wt and this produces a 
current i in the stator winding and also forces a current i { to flow in the field 
winding in order to sustain the field flux interlinkage. The transient currents 
i and i, are determined by the angle a as well as the leakage fluxes of both 

cL f 

windings and they become a maximum when a = V/2, a quarter period after the 
short circuit occurred. 


It can be shown that the maximum transient stator current is determined by the 
equivalent circuit of Figure 12, and the reactance that corresponds to this cir¬ 
cuit is the direct axis unsaturated transient reactance, X du - The field reac¬ 
tance, X„, is given in Kilgore's paper as 
r 

= X ^C?, A F 
F Tr M 

and therefore the unsaturated transient reactance becomes 


The saturated transient reactance has been determined empirically and is ap¬ 
proximately 88% of the unsaturated value. 


In the determination of the transient reactance only the field winding of the 
rotor was considered. If there is a damper winding in the poles, or if eddy 
currents are possible in the rotor iron, subtransient currents similar to the 
transient currents of the field winding are induced in these circuits. These 
circuits will support the field winding for a few cycles and therefore they have 
to be considered as acting in parallel with the field winding. The equivalent 
circuit for this case is shown in Figure 13. X Dd is the leakage reactance of 
the damper winding and eddy current circuits together in the direct axis, and 
per Kilgore’s paper 



= X /Dd 
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wh ere x Dd = [g + <f d + h r2 ] 

S is a depth of penetration factor and varies as . It is equal to 1. 2 at 60 
cycles and (]/^) (l.2) at 400 cycles. The subtransient reactance is there¬ 


fore 


X d ■ + X Dd 


The rate of decrease of the transient and subtransient currents will be deter¬ 
mined by the time constants of the windings involved. The damper winding 
and eddy current circuits have much larger ratios of resistance to leakage 
reactance than the field winding and therefore their influence will be much 
shorter in duration than that of the field winding. As a matter of fact, the 
damper winding and eddy current circuits influence the currents only during 
tlie first few cycles. The field winding determines the decrease of the ampli¬ 
tudes for :i much longer time. The change of the amplitudes during the short 

circuit period is such that the amplitudes are determined first by the subtran- 
« * 

sient reactance then by the transient reactance X^, and finally by the 
synchronous reactance X .. 


The total self inductance of the field winding is 


L 


f 





where is the ratio of the interlinkages of the field with its own flux to the 
maximum interlinkages that would be produced with a uniform gap and a con¬ 
centrated field winding. 
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TIME CONSTANTS 


The time constant is the time in seconds required for the particular component 
to decay to 36.8% of its initial value. The time constant T^ is the time con¬ 
stant of the field winding with the armature circuit open and with negligible 
external resistance and inductance in the field circuit. Therefore, the open 
circuit time constant is 



The armature time constant is the time constant of the DC component and is 



where r 

a 


stator I 2 R (KW) 
rated KVA 


The transient time constant T^ is the time constant of the transient reactance 
component of the alternating wave and with good approximation is 



The subtransient time constant T d is the time constant of the subtransient 

reactance component ol the alternating wave and is approximately . 005 second 

lor 400 cycle machines, (from tests of 60 cycle machines t" 035 second or 
2.1 cycles). 
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Resistances 


Positive Sequence Resistance rj 

The positive sequence resistance of an a-c machine is its a-c 
armature resistance r which is greater than the d-c armature 
resistance r a because of non-uniform current density and local 
iron loss, m is usually disregarded in system studies except 
in the case of a 3 phase short circuit near the generator terminals. 


Zero-Sequence Resistance r f 


This resistance is small and is usually neglected unless a 
lesistor is used in the neutral of a wye connected machine. 

It is the sum oi the d-c resistance of the three places (in 
a 3 phase machine). If a neutral resistor is used, add 3 x the 
resistor also. 


gat i 


o Sequence Resistance v 0 


The negative sequence power input to the a 
1/2 of the copper losses of the rotor c ircuits, 
rotor copper losses are supplied mechanically. 


r:nature supplies 
The other 1/2 of 


The negative sequence resistance depends upon the re¬ 
sistance? of the damper circuit If there is one. It is equal to 
3/2 the rotor damper circuit resistance plus the stator 
resistance 

r 2 = R s + ^ R r 
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GENERATOR VOLTAGE AND OUTPUT EQUATIONS 


To begin a generator design when there is no required envelope or other 

size limitation, time can be saved by using an output formula that relates 

rotor size to output kva. The output equation is given here and then devel- 

-8 

oped from the relation E = N £v B 10” 

c max 

output kva = d r -g (g£M)A Bg 

90 x 10' 

where d f = rotor diameter which is for practical purposes the same as the 

stator I. D. (inches) 

= stator length in inches 

A = ampere wires/inch of bore periphery 

2 

Bg = gap density in kilo lines/in 


All electromagnetic generators convert mechanical energy to electrical 
energy through a change of flux linking the conductors of the output winding. 


The voltage generated in a single conductor is - 

10 ' 


E =-£ v B 


max 

e/coil = 2 Emax sin wt Kp 


„ , 6?rms Coils in Series „ 

Eph = - C5ir x —— xK d 

_erms. t Total Slots K d 
Coil Phases x Parallels 


2^vB max 10^ x x N 

me s 


f .r 
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(N = Conductors/Slot) 
s 

0 T = Bg7Td r ^ 

B max Fundamental = Cj Bg 
^ “max Fund = C, 0< T X 


n-mi K N Q K , 10 
E phase = C A <? T ^ x - g 5 d 


-8 


& 


me 


p rh RPM XT ^ v _ 1ft -8 
C L 60 N K d 10 

if2m 

(Ne = total effective conductors in the stator = 

^ _ 60 Eph (2 m _ 60 E LL ^ m 1()8 

RPM Ne X C, K , RPM Ne C, K , E T T 
id 1 d LL 

Eph 


Q N K 
s p 


C 


) 


C w = Winding Constant = 

Ifi m 

60 x 10 5 x E T T 

d = 

V T RPM C Ne 
w 

C for 3 phase = . 39 C, 

W 1 

C for 2 phase (90°) = . 319 C, 

W 1 


C for 3 phase = . 225 C, 
w 1 
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The derived voltage equation can be found in the AIEE Paper No. 50-201 
dated August, 1950, "Design Calculations for A. C. Generators” by David 
Ginsberg, 



Where = Hypothetical total flux in the air gap of the machine (in 

kilolines). This is the maximum density over the pole 
times the area of the stator bore. 

N = Total effective conductors in the machine, 
e 

RPM = Revolutions per minute. 

C w = Winding Constant 


LL 
m E . 


c i K d 


-— - . 39 Ct for a 3 phase machine 

n 1 


Where m = No. of phases 

K, = Distribution Factor = .955 
d 

C = Ratio — ax ^ n . ¥ — Fundamenta l fj ux wave -phis va i ue 

1 Actual Maximum 

is 1.0 for a sine wave and a sine wave is assumed in the 
formula. 


Then 


KVA 


i l e ll 


fa 
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A'TTd 

N 

e 


since A = 


N I u 
e ph 

IT d 


ampere wires per inch of stator 
bore periphery 


B (7Td 1 )RPM(. 39C,)N |[3 A 7Td 

KVA = -i---. 4 e _ 

60 x 10 8 N 

e 

B 9.85d 2 l RPM .39 C, A (3 
60 x 10 8 


B (d 2 l )(RPM)A 
KVA = -I-£___ 

90 x 10 7 


The basic equation. 


If a pole embrace is used that gives other than Cj = 1.0 (or if the maximum 
of the fundamental flux wave is not equal to the maximum of the actual flux 
wave), the KVA formula should be multiplied by Cj. 

Note. Liwschitz-Garik and Whipple in ’’Electric Machinery” Vol. 2, 

appendix 8, pp 553-557 derived an equivalent formula for either 
motors or generators. 
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DIRECT AXIS COMPONENT CF ARMATURE REACTION 


(¥> 


X -o 








'( 1 r)’ r 


m* 1 


When the air gap under the pole is constant the field produced by the 
direct axis armature reaction MMF wave A will have the approximate shape B in¬ 
sofar as its reaction on the field pole is concerned. Since this MMF produces 
very little demagnetizing effect in the interpolar region it is sufficient to 
consider only the fundamental of the effective part of the MMF curve (the part 
which lies u n de r the pole). The amplitude of this fundamental is given by the 


equation 


, r irr 

A,~ rr ] f(*) s/ * * 

o 


If the ratio of the pole arc to the pole pitch is designated as °c then 
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Although this derivation has been based on the assumption that the air 
gap under the pole is constant, the formula is much more accurate than might 
appear# The formula has been checked by flux plots over a wide range of pole 
shapes from k to 88 poles and found to be reasonable accurate 
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QUADRATURE AXIS COMPONENT OF ARMATURE REACTION 


(/- 
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X = 0 


For all practical purposes the armature MMF in the quadrature axis will 
hare the shape indicated insofar as its reaction on the field is concerned. If 
the assunption ( justified by practice ) is made that the magnitude of the con- 
stant part over the interpolar space is equal to 1/8 , the fundamental of 

this wave will then be found by the equation 
A, Tf cos x Jx 

f( x ) " g from X ~ O to X ~ (/ - oc) 7 ^ 

f(x) - f^ M cos X from ( /- <x)ff fo (/ ocJslL 

f ( x ) ~ Q f rorr> ( f+ 0(.J ^ to 7f ~* { / ~ 

t( x ) ~ f~E,M Co s X from 7T + (/ — CL.) fx. to 7T t* // y- <X.J 

g fpv, from 7T + (/+<*•)% o 2 fT 
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£OS X d V ^ Tr 
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VECTOR DIAGRAM AND CALCULATION OF LOAD 
EXCITATION OF A SALIENT POLE GENERATOR 



Th© vector diagram for a salient pole generator is as shown above. OA is 
the terminal voltage E ph and OB is the phase current I phi drawn at the 
proper power factor angle 6 . AC Is the effective resistance drop l ph r e 

and CD is the leakage reactance drop l ph 1 £ . The stator current is di¬ 

vided into a direct axis component 1^ and a quadrature axis component j , j 
and I , X , will thus be the voltage induced in the stator by the direct 

4 *-4 

axis flux and 1^ X „ will be the voltage induced in the stator by 

the cross flux ^ . The voltage I ^ X will be in quadrature with 

and it is represented by DE. Likewise, the voltage I j X will 
be in quadrature with and it is represented by EF. The voltage CF is 

the voltage that would exist at no load with excitation corresponding to rated 
load, neglecting saturation, and this voltage is represented by the symbol 


LA-lot 



Consider the triangle 01^ B and 

am ^“ n * 

ces & z ‘ Is ; X? ^ S?h cc,s 

The triangles 01^ B and DEG are similar because E0 is perpendicular to 
B, GD is perpendicular to BO, and DE is perpendicular to 01 f . Thus 

angle 0 is equal to angle D. Also, triangles DBG and HFO are similar because 
DE is parallel to HF, EG is concentric with FG, and GD is concentric with GH, 
and therefore angle D and angle H are equal. From a similar triangle relation¬ 
ship the angle A of triangle AM3 is also equal to . 

From triangle DHN 


S/N 
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£>0 
D H 



DH = 

Jj 
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Tj A »s 
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Cjs 

K F 

= CH 
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• j Tph Xj 

5 > /■ 
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os €r 

~ E/>h 

c o s € 


/I AT; 

■ LK 

J/ot) ■ 

r e s 




Thus the nominal voltage e^- is derived as 

OLf LX * JXP 

*4 » £f=h ccs e + !?>*>(’* & * X* S/V ft) 


In the per unit system of notation B^ and I^ are unity at rated 
load conditions and if resistance is neglected the nominal voltage equation 
becomes 


<9f 


€ * XJ 
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From triangle DBG 


Cos 


V =" 
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DE 
D<£ 

-4? ^ _ 


4 >4 


Z) cs 

-*f 4f 




CO "*> lf/'U _ r „ 

C Cj = CD + D<5 sr Tfh y- - 
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Therefore, if resistance is again neglected the angles 0 and p can 
be determined by 

Tan /!/ =■ ^ ^ -S/*^ _ \ f * S‘* £ 

op cot a & 

Ihen load is applied to the generator the flux wave shape will be distort 
ed and the wave form under load will be less effective in generating Yoltage 
than the no load wave form. The reason for this is that the fundamental of 
the wave form under load is relatively less than that of the no load form for 
an equal total flux per pole under each condition. Therefore the flux per pol 
under load conditions will be greater than that required at no load, and the 
flux per pole under load is 

fipt s &>[ € s -- 53 fij 

The leakage flux under load will also increase and 

JZ^X - a£ y rf V £ *■ 0 + cos 8) 4 + Fc. 

L ^ '7'fy r 7fr 

The total flux per pole in the rotor thus becomes 

fiftt ~ fipt + fift 

and the density of the pole under load is 



The total ampere turns per pole under load can thus be calculated as 



SINGLE PHASE CALCULATIONS 

When two phases of a normal three phase generator are used for single 

phase operation the machine reactances are first calculated on the basis of a 

three phase machine and the single phase reactances then become* 

;ja-io8 



*</('+) 

Kf(^) 

0 $) 

The negative sequence component of the single phase armature reaction 
wave will induce large currents in the damper winding and this will cause high 
damper losses under load conditions. Until further experience is accumulated 
it is not felt advisable to exceed current densities of 10,000 arapa/sq. in. 
in the damper bars under normal load conditions when the current is calculated 
in the following manner. 

tfes. Smq, AT.= B atf (3 pfmte )] 
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SECTION TA 










THE EFFECT OF INCREASING THE AIR-GAP 


The losses and reactances of a 30 KVA, 320 CPS, 4800 rpm aircraft 
generator were calculated for various air-gap lengths. The stator 
and its output winding were both held unchanged for the range of air- 
gap lengths. 

As would be expected under these conditions, the field excitation 
requirement soon limits the output and, for this particular machine, 
the thermal limit for continuous full load operation would be reached 
when the air-gap was increased from .035” to .055”, an increase of 
57%. 


TA-l 





TO '? INCH 46 0860 

















ampere jwires 


















SALIENT POLE SYNCHRONOUS DESIGN SHEET 


STATOR 


ROTOR 


Core Length_ 

DBS x 2 _ 

s lot • i± l '" 

. .122 x 

Sise Slots 


n. 

9.25 


3.0.. 


Type Wdg. 


83.3% 


Throw _ 

Skew h Dlst. Fact. 

Chord Fact._ 

Cond. Per Slot_ 


ser. Y 
1 - 11 


Cond. Per 8 lot__ 

Total Eff. Cond-l8%.8_ 

.075 x .162 SGHF 


Cond. Sise_ 

Cond. Area_ 

ns 11 


wdg . Cons t. 
Total Fltpi^ 

Gap Area_ 

Oap Dens i ty 
Pole Const. 


.0119*+ 


SATURATION 


Air Gap A.T._ 

Stator A.T.__ 

Pole A.T- 

Ho. Damper Bars- No Load A.T.— 

Bar Site dla> _ Bated AT 

Bar P<trh »*+03 h .038 b __ .040 Overload A.T._ 

Shor t Ci r c . A . 


0322 


9 


Wt. Of Copper 


% Load 


LOSStS-IHICIINCY 


% Load 


Amos/In 


3600 7650 


1/2 Mean Turn_ 

Rea . Per Ph . At _ jj 

Eddy Fact. Top- 

Eddy Fact. Bottom 


Dema g . Fact. Cm_ 

Amp. Cond. Per In. 

React. Factor_ 

Cond. Perm.- 

End Perm.- 

Leakage React.- 

Air Gap Perm.- 

React . Of Arm. Xad 


6.99 

300*F To59l 


1.026 


RIACT..TIMI CONST. 


Unsat. Trans._ C 

Sat. Trans.- \ 

Subtrans. Xd. Je2*3Q Xq—1 

Neg. Seauence-- 

Zero Sequence---~—~ 

300° F 

Open Circ. Time Con.- 

Arm. Time Con- 

Trans. Time Con.- 

Subtrans. Time Con.— 


L Losses 


89.2 


Stator Vatts/In 
Rotor Watts/In. 


Volts _ 83 _Amps. 


.Type Cooling 

5_Ph. 320 C vc le. _J+800_RPM 









































P ole f face ;Loss ^urvesj J*or A Wide Speed Range 


Solid-Pole! 40 KfA Genjeratoif Operating at 
8 ,00 0 RPM |and wjth a W «L. <|ap Dehsity iof 
28 K,6/in^j !jj 


-f- 



JPull Load T 

j_ ^ _A 

sr | I 

4 — 1 — 4 - 


! ^ 

—U -4- 


No toad 


! 

; 

_J 

L_ 

1 

j 

! 

i 

! 

! 

i 

! 

i 





















-FACE LdSSES 



TA-3 














INCH 46 0860 



TA-4 



46 0860 










